Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,633 G→A W580* (TGG→TGA cstA → pyruvate transporter CstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,6330GA98.2% 194.1 / ‑6.1 55W580* (TGG→TGAcstApyruvate transporter CstA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/24);  minor base C (1/0);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCG  >  NC_000913/631383‑631859
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
tGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGGATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGGGTTGTTCAAGATGAAACGTCAACGTTACGCCTgg                                                                                                                                                                                                                                    >  2:38678/1‑251 (MQ=255)
             gctgGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTgg                                                                                                                                                                                                                          <  1:56712/248‑1 (MQ=255)
             gctgGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTAc                                                                                                                                                                                                                       >  1:117764/1‑251 (MQ=255)
               tgGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCAGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCATATGCTGGCAGGGTTGGCGCGGATGCTCTGTGCCGGGGGGTTGTTCAAGATGAACCTTCAACATTACGCCTGAGTGGCGCTGGtgcc                                                                                                                                                                                                                      >  2:341608/1‑247 (MQ=255)
                    gtggtgTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGc                                                                                                                                                                                                                <  1:365101/251‑1 (MQ=255)
                           ctcCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGCTCCGCTGGGCGGCATTAACACCCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGgctg                                                                                                                                                                                                         >  1:284247/1‑251 (MQ=255)
                           ctcCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGAGGCTCTGGGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGgctg                                                                                                                                                                                                         >  1:151664/1‑251 (MQ=255)
                                       aGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTAcc                                                                                                                                                                                             >  2:99599/1‑251 (MQ=255)
                                            aCCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGAc                                                                                                                                                                                        <  1:369986/251‑1 (MQ=255)
                                             ccGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGAcc                                                                                                                                                                                       <  2:240937/251‑1 (MQ=255)
                                                 ttCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTAGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGc                                                                                                                                                                                     >  1:203411/1‑249 (MQ=255)
                                                   cACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCCGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCACGGCTGCTGATTTGTACCCTGACCGCAGGc                                                                                                                                                                                 >  2:205959/1‑251 (MQ=255)
                                                       gcctgcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTg                                                                                                                                                                               <  2:253126/249‑1 (MQ=255)
                                                         ctgcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCCACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGccg                                                                                                                                                                           >  2:134965/1‑249 (MQ=255)
                                                           gcTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCGCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGaa                                                                                                                                                                         <  2:166974/251‑1 (MQ=255)
                                                                         aaCAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCgg                                                                                                                                                           <  1:100887/251‑1 (MQ=255)
                                                                                      gTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCgg                                                                                                                                              <  2:273504/251‑1 (MQ=255)
                                                                                         cTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCtt                                                                                                                                           >  1:215959/1‑251 (MQ=255)
                                                                                           ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTcc                                                                                                                                         <  1:99599/251‑1 (MQ=255)
                                                                                           ggCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTcc                                                                                                                                         >  1:246459/1‑251 (MQ=255)
                                                                                            gCGTGTGGCTACTTCCACCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCt                                                                                                                                        <  2:117764/251‑1 (MQ=255)
                                                                                                  ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCAt                                                                                                                                  <  1:121319/251‑1 (MQ=255)
                                                                                                  ggCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCAt                                                                                                                                  <  2:184877/251‑1 (MQ=255)
                                                                                                      aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCt                                                                                                                              >  2:321802/1‑251 (MQ=255)
                                                                                                      aCTTCCTCCATCAGGGCGTGGTCGATCCGCTGCGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGACTGGCTGCTGATTTGTACCCTGACAGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCt                                                                                                                              >  1:365071/1‑251 (MQ=255)
                                                                                                         tcctccATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCta                                                                                                                             <  1:134965/249‑1 (MQ=255)
                                                                                                                gcAGGGCGGGGTCGTTCGGCTGGGCGGCAGTAACACTCTGTGGCCGCGGTTCGGTATTTCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGTTTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTAAGCCCGGATGTTAAAGTCGGCGTCCGGGCCAttcagaataatttcc                                                                                                                    <  2:39275/250‑14 (MQ=255)
                                                                                                                          gTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGACTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGAt                                                                                                          <  2:267210/251‑1 (MQ=255)
                                                                                                                          gTCGATCCGCTGGGCGACATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGAt                                                                                                          >  1:364291/1‑251 (MQ=255)
                                                                                                                            cGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATc                                                                                                         <  2:189338/250‑1 (MQ=255)
                                                                                                                            cGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCg                                                                                                        >  1:400520/1‑251 (MQ=255)
                                                                                                                                ccGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGAc                                                                                                      >  2:287772/1‑249 (MQ=255)
                                                                                                                                 cGCTGGGGGGCCTTAACCCTCTGTGGCCGCTGTTTGGGATTGCCAACCAGATGCTGGCAGGGTTGGCGCTGATGCTCTGTGCCGTGGGGTTGTTCACGATGAAACGTCCACGTTACGCCTGCGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCCGAAAGCGGTTAGCCCAGATGCGAAAGTCGGCTTCCTGGCCATTGCTACTAAGGTCCAGGCAATGATAGACAgg                                                                                                   >  2:456993/1‑250 (MQ=255)
                                                                                                                                 cGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGc                                                                                                   >  1:18110/1‑251 (MQ=255)
                                                                                                                                            aTTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTcc                                                                                        <  1:60452/251‑1 (MQ=255)
                                                                                                                                                   ctctGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGc                                                                                   >  1:43271/1‑249 (MQ=255)
                                                                                                                                                         ggCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATAc                                                                            <  1:361487/250‑1 (MQ=255)
                                                                                                                                                         ggCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACt                                                                           >  2:185729/1‑251 (MQ=255)
                                                                                                                                                                 tttGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGCCAGCGGCAATATTCCGTCGCAGTATACTGAGTcac                                                                    >  2:398872/1‑250 (MQ=255)
                                                                                                                                                                  ttGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTcaca                                                                   <  1:23313/250‑1 (MQ=255)
                                                                                                                                                                     gTATTGCCATCCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCTTGGTGTTGTTCAGGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGTCCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTg                                                               <  2:284247/251‑1 (MQ=255)
                                                                                                                                                                                 aGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATAATGAGTCACAGCTGGCGCAACTGgtg                                                   >  2:256406/1‑251 (MQ=255)
                                                                                                                                                                                   atgttGGCAGGGATGGCGCTGATGCGCTGTGCCGTGGTGTTGTTCACGAGGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCCACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCCGGCAGAAAGCGTTTAACCCGGATGCGAAAGACGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATCTTCCGGCGCAGTATACGGAGTCACAGCTGGCGCAACTGGNGtt                                                 >  2:409582/5‑251 (MQ=255)
                                                                                                                                                                                       tGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTcaac                                             >  1:144585/1‑251 (MQ=255)
                                                                                                                                                                                          cAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCaacaac                                          >  1:379682/1‑251 (MQ=255)
                                                                                                                                                                                                  gcgatGATGCTCTGTGCCGTGGTGGTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGCACCAACGGCCTGGTTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAACCGTTTAGCCCGGATGAGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGAATACTGGGTCACAGCTGGCGCAACTGGTGTTCAAAAACCggtcgg                                   >  2:72387/5‑245 (MQ=255)
                                                                                                                                                                                                   cgcTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGAt                                 <  1:205959/251‑1 (MQ=255)
                                                                                                                                                                                                       gATGCGCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCg                             <  1:219585/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACNAtc                  >  2:415334/1‑251 (MQ=255)
                                                                                                                                                                                                                       gtgtTGTTCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATcttctt             >  2:341262/1‑251 (MQ=255)
                                                                                                                                                                                                                            gttCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtgg        >  1:233875/1‑251 (MQ=255)
                                                                                                                                                                                                                             ttCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtggt       >  1:70375/1‑251 (MQ=255)
                                                                                                                                                                                                                              tCAAGATGAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTAtggtg      >  2:169396/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gATNAAACGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCg  <  2:421089/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaCGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCa                                                                                               >  2:254440/1‑154 (MQ=255)
                                                                                                                                                                                                                                      aaaCGTCAACGTTACGCCTGAGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCa                                                                                               <  1:254440/154‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
TGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATTCACTGCCTGCTAACCTGCTGGCAACAGCGCTGTGCGTACTGGCGTGGGGCTACTTCCTCCATCAGGGCGTGGTCGATCCGCTGGGCGGCATTAACACTCTGTGGCCGCTGTTTGGTATTGCCAACCAGATGCTGGCAGGGATGGCGCTGATGCTCTGTGCCGTGGTGTTGTTCAAGATGAAACGTCAACGTTACGCCTGGGTGGCGCTGGTACCAACGGCCTGGCTGCTGATTTGTACCCTGACCGCAGGCTGGCAGAAAGCGTTTAGCCCGGATGCGAAAGTCGGCTTCCTGGCCATTGCTAATAAGTTCCAGGCAATGATCGACAGCGGCAATATTCCGTCGCAGTATACTGAGTCACAGCTGGCGCAACTGGTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCG  >  NC_000913/631383‑631859

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