Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,587,531 G→A L195F (CTC→TTC)  ugpQ ← glycerophosphodiester phosphodiesterase UgpQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,5310GA98.4% 203.6 / NA 64L195F (CTC→TTC) ugpQglycerophosphodiester phosphodiesterase UgpQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (41/22);  minor base T (1/0);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCG  >  NC_000913/3587285‑3587757
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    
aaTCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCTGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAAc                                                                                                                                                                                                                                  <  1:331958/251‑1 (MQ=255)
    agcagcATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTAt                                                                                                                                                                                                                              <  1:53002/251‑1 (MQ=255)
         caTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTg                                                                                                                                                                                                                         <  1:204405/251‑1 (MQ=255)
               ctgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCCTCGGTACAGATGCAATCCACCACCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCCCTCGCGCTTTATTGAATAACTTATGATTGGGATg                                                                                                                                                                                                                    >  2:373745/1‑250 (MQ=255)
                  ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAAt                                                                                                                                                                                                                 >  1:281749/1‑250 (MQ=255)
                       gaggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGGATAGAGAc                                                                                                                                                                                                           >  1:257391/1‑251 (MQ=255)
                        aggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCACTCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACg                                                                                                                                                                                                          >  1:270165/1‑251 (MQ=255)
                         ggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCACTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGc                                                                                                                                                                                                         >  2:424737/1‑251 (MQ=255)
                           agaAACGACGACATGCTGAATCAGCCCAGCACGACGCAGCCTGACCTTCCGTTGAAAACTATTGGGCCGTCAAGTTCCGCCCAATCACGTCAATCGCATCGGCACAGCTGCACTCCACACCCCCGCGCAGCAACCCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGACACCGAAGTCCGGCGTCTTTCCACTGCATCACTCGCGCATTCTTGATTAACttgtgctgtagataacaagcaacgc                                                                                                                                                                                                         >  2:241725/1‑226 (MQ=255)
                                          cTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTc                                                                                                                                                                                        >  1:14534/1‑251 (MQ=255)
                                                     cAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCa                                                                                                                                                                             <  1:373745/251‑1 (MQ=255)
                                                        cACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAgtc                                                                                                                                                                          <  1:175607/251‑1 (MQ=255)
                                                        aaCGCCGCAGCCGTACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAAGCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGAGTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAgtc                                                                                                                                                                          <  2:543005/250‑1 (MQ=255)
                                                            ccgcAGCTTGACATTCCGTTGAAAACTAGTGGGCCGTAAAGTTCTGGCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGtcgtcg                                                                                                                                                                      <  1:315695/251‑1 (MQ=255)
                                                                   cTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACt                                                                                                                                                               >  1:413020/1‑251 (MQ=255)
                                                                        aTTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGCTACGCAGCCCAGCCGCGCGGTCAGTTAGCGCCAGTCATCGAGcccataatcc                                                                                                                                                          >  2:478302/1‑243 (MQ=255)
                                                                                    aCTATTGGGCCGTAAAGTACGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑cg                                                                                                                                            <  1:316218/251‑1 (MQ=255)
                                                                                       aTTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGGCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑cgc                                                                                                                                           <  1:458528/249‑1 (MQ=255)
                                                                                         tGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAAACAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACCAAC‑‑CGCGCgg                                                                                                                                       >  1:261741/1‑251 (MQ=255)
                                                                                           ggCCGTAAAGTTCGGACCAATCCCGTCTAACGTAGCGGTCCAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTAGTGACGGTATAAACCAGAATCCGCAGTCCGGAGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTAGACGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCa                                                                                                                                     <  2:382215/251‑1 (MQ=255)
                                                                                           ggCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCa                                                                                                                                     <  1:424737/251‑1 (MQ=255)
                                                                                                 aaagggcGGACCAATCACGTCAATCGTATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCGTCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCTCAGTCCGGCGTCTTTCAACTGCATCACCCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTGCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTcc                                                                                                                                <  2:382270/244‑1 (MQ=255)
                                                                                                 aaaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGCATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCg                                                                                                                               >  1:34316/1‑251 (MQ=255)
                                                                                                 aaaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTcc                                                                                                                                <  1:488677/250‑1 (MQ=255)
                                                                                                 aaaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCACAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CTCGCGGCAGTTccc                                                                                                                               >  2:135446/1‑250 (MQ=255)
                                                                                                       cGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTgccg                                                                                                                         >  1:425015/1‑251 (MQ=255)
                                                                                                          aCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTgccgccc                                                                                                                      >  1:103608/1‑250 (MQ=255)
                                                                                                                 aCGTCAATCGCATCGGTAAAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAACCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCCCTCATCCCGCACCACAC‑‑CGCACGGCAGTCCCGGTGCCGCCtgttgttc                                                                                                               >  2:384509/1‑249 (MQ=255)
                                                                                                                  cGTCAATCGCGACGATACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGGCCGACGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGCTGAATAGAGACGCCGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCATAAGGACCGATGCCGCCTGTTGTGc                                                                                                               >  2:492561/1‑250 (MQ=255)
                                                                                                                      aaTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCtt                                                                                                          >  2:67305/1‑251 (MQ=255)
                                                                                                                      aaTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCtt                                                                                                          >  1:259186/1‑251 (MQ=255)
                                                                                                                          gCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGCATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTa                                                                                                       >  1:324436/1‑250 (MQ=255)
                                                                                                                           cATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTaa                                                                                                      >  1:348292/1‑250 (MQ=255)
                                                                                                                           cATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAATCGTCGCGCCACTCATACAGCAACAAAA‑‑CGCGCGGCAGTTCCGTTGCCGCCTGATGTGCAGCTTCTaa                                                                                                      >  2:487876/1‑250 (MQ=255)
                                                                                                                              cGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCa                                                                                                  <  1:190705/251‑1 (MQ=255)
                                                                                                                                gTACAGATTCGATCCACCCCCCAGCGCAGCAACGCTGCTACGCGCTGGGGTATATTGACGGAATAAACCAGAATCCGCAGTCAGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGCAAAACTTCTGATTGAGAAGAATAGAAACGCAGCACAGCCGCGCGGTCAGTTCGCGCCAGGAGTCGCGCAACTCATCCAG‑‑ACAAACAACGCGCGGCAGGTCCTGTGCCACGTGTTGGGCAGCTTAGAAAGCAt                                                                                                 >  2:185168/1‑250 (MQ=255)
                                                                                                                                 tACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGcc                                                                                                                                                                              >  1:413788/1‑174 (MQ=255)
                                                                                                                                 gaCAGATGAAATCCACACCCCAGCGCAGCAACTCTGCGGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCTGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGAGTGAGATGAATAGAGACGCAGTCCATCCGCGCGGTCAGTTCGCGcc                                                                                                                                                                              <  2:413788/173‑1 (MQ=255)
                                                                                                                                  aCAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGGCGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGATTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAA‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATcc                                                                                               >  1:237650/1‑249 (MQ=255)
                                                                                                                                          aaTCCACACCCCCGCGCAGCAACTCTGCTGCGCGCTGGAGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGAGTCTTTCAACTGCATCACTCGCGCTGTCTTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCAAGCCGCTCGGACAGTTCGCGCAAGTCGTCGCGCCACTAATCCAGCAACAAAC‑‑CGCGCGGCAGTGCCGATGCCGCCTGTTAGACAGCTTCTAAAGCATCAAtccataa                                                                                       >  2:464158/1‑245 (MQ=255)
                                                                                                                                           aTCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCCCTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGTACAGTTCGCGCCAGTCGTCGCGCCCCTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCCGTTGTGCAGCTTCTAAAGCATCACTCTCAAACg                                                                                     >  2:285759/1‑251 (MQ=255)
                                                                                                                                                 cccccAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAAGCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGGCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCATTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGAca                                                                               <  1:519113/250‑1 (MQ=255)
                                                                                                                                                   cccAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTCTTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGGCAGTTCACGCCAGTCGTCGCGCCACTCATCCAGCACCAACC‑‑CGCGCGGCAGTTACCGTGCCGCCTGTTGTGCAACATCTAAAGCCTCAATCTCAAACGATGActac                                                                             >  2:61051/1‑248 (MQ=255)
                                                                                                                                                      aacgcAGCAACTCGGCTGGGCGCTGGGGTTTATTGACGTTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGATCTTTATTGAATAACTTATGATTGAGAGGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTACCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACagcagc                                                                          <  2:224281/249‑1 (MQ=255)
                                                                                                                                                                 tctGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCAt                                                               >  1:28300/1‑251 (MQ=255)
                                                                                                                                                                      tgcGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCCGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTAATTGAATAACTTATGATTGAGATGAATAGAGACGCAACCCAGACGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACCCATCCAGCAACAAAC‑‑CGCCCGGCCGATTCGGAGCCGCCTCTTGTGACGCTTATAAAGCATTCATCTCAACAGATGACATACGCGGCGGCGTCATaacg                                                           >  2:43691/1‑247 (MQ=255)
                                                                                                                                                                        cgcgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGCTTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCACCAAAC‑‑CGAGCGGCAGTTCCGGTGCCGCCTGTTGGGCAGCTTCTAAAGCATCAATCTAAACCGATGACAGCATCGGCGGCGTCATAACGGcc                                                        >  2:472904/1‑251 (MQ=255)
                                                                                                                                                                           gcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCgcgc                                                                                                                                                                               >  1:213492/1‑131 (MQ=255)
                                                                                                                                                                           gcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCgcgc                                                                                                                                                                               <  2:213492/131‑1 (MQ=255)
                                                                                                                                                                                  tttATTGACGGTATAAACCAGACTCCGCAGTCCGGCGTCTTTCAACTGCATCACTAGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCACCAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATAAATCTCAAACGATGACAGCAGCGGCGGCGTCATAACGGCCCACAGCTcgc                                              >  1:563706/1‑251 (MQ=255)
                                                                                                                                                                                     aTTGACGGTATAATCCAGAATCCGCAGTCCGGCGGCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATTAATAGATACGcagcccagcc                                                                                                                                                                                                <  2:345507/104‑1 (MQ=255)
                                                                                                                                                                                     aTTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGcagcccagcc                                                                                                                                                                                                >  1:345507/1‑104 (MQ=255)
                                                                                                                                                                                     aTTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTg                                           <  1:451250/251‑1 (MQ=255)
                                                                                                                                                                                           ggTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCa                                     <  1:195781/251‑1 (MQ=255)
                                                                                                                                                                                            gTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTAGCGTgccgcc                                      >  1:339581/1‑249 (MQ=255)
                                                                                                                                                                                                aaaCCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCACCAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCCGCGGCGGCGTCATACCGGCCCACAGCTCCCGTGCCGCCAGCGcc                                >  2:429887/1‑251 (MQ=255)
                                                                                                                                                                                                        aTCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCAtttt                        <  1:67305/251‑1 (MQ=255)
                                                                                                                                                                                                              aGTCCGGCGTCTTTCAACTGCATAACTCGAGCTTTATTGAATAACCTATGCTTGAGATGAATAGAGACCCAGCCCAGCCGCGCGGGCAGTTCGCGCCAGTCGTCGCGCCACTCATCCCGCAACAAAC‑‑CGCGAGGCAGTTCCGGTGCCGCCTGTTGTGCCGCTTCTAAAGCCTCAATCTCAAAAGAGTACAGCAGCGGCGGCGCCATACCGGCCCACAGCTCGCGCTCCGCACGCGACCCCATTTTGCCCgg                  >  2:505910/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATACATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCAAAGCTCGAGTGCAGCCAGCGCCACCATTTTGCCCGTTAAt              >  1:370290/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggCGTCTTTCAACTGCATCACTCGCGATTTATTGAATACCTTATGATTGCGATGAATAGAGACGCAGCCAAGCCGCGCGTTCAGTTCGCGCCAGTCGTCGCGCCACTCATACACCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGCTGACAGCAGCGGCGGCGTCTTACCGGCCCAAAGCTCACGTACAGCCGCCCCCACCATTTTCCCCCTTattg             >  2:373080/1‑248 (MQ=255)
                                                                                                                                                                                                                     cGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGc           >  1:218681/1‑251 (MQ=255)
                                                                                                                                                                                                                           tCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCCTTTTGCCCGTTAATGGCccggtg     >  1:217920/1‑251 (MQ=255)
                                                                                                                                                                                                                            cAACTGAATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGAGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCACCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCCGAGCCCCCCATTTGCCCGtaacaggccaggggc    >  2:232568/1‑237 (MQ=255)
                                                                                                                                                                                                                              aCTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggtgccg  >  1:477651/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    
AATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAAC‑‑CGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCG  >  NC_000913/3587285‑3587757

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: