Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,202 A→T 5.1% intergenic (+45/+26) smrB → / ← yfcO putative DNA endonuclease/DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,2020AT 5.1% 395.9 / 5.3 138intergenic (+45/+26)smrB/yfcOputative DNA endonuclease/DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (59/71);  new base T (2/6);  total (61/77)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

CGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTG  >  NC_000913/2449054‑2449351
                                                                                                                                                    |                                                                                                                                                     
cGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGAATGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAt                                                                                                                                                      <  1:890399/150‑1 (MQ=255)
  cATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGc                                                                                                                                                    >  2:48036/1‑150 (MQ=255)
     gtgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcccc                                                                                                                                                 >  1:1748691/1‑150 (MQ=255)
     gggATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGGGAGGCGGCGTGGGTGGGGTTGATTGAAGTCGATGAGTGGCGGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcccc                                                                                                                                                 <  1:980154/148‑1 (MQ=255)
      tgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcct                                                                                                                                                >  2:34781/1‑150 (MQ=255)
      tgATGGCCGTTCATCATGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATAGAAGTCGAAGAGGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcct                                                                                                                                                <  2:1204585/150‑1 (MQ=255)
        aTGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgc                                                                                                                                              <  2:569631/150‑1 (MQ=255)
         tGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgcc                                                                                                                                             <  1:445223/150‑1 (MQ=255)
         tGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgcc                                                                                                                                             <  2:1838215/150‑1 (MQ=255)
         tGGCCTTACATCAGGAACCACAAGACTAGGCCGGTGATGCGGCGTTGTTGGTGGTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATCGCAGGGGCATTCGCCcctgcc                                                                                                                                             <  2:921578/150‑1 (MQ=255)
            ccTTTCATCAGGCACCAAACGAGTATGGCGGTGATGCGGCGTGGTTGGTGTTGATTGAAGTCGAAGAGGGGCTGCCGCCGGAGTGGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTgg                                                                                                                                          <  1:1417098/150‑1 (MQ=255)
            ccTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTgg                                                                                                                                          <  2:572611/150‑1 (MQ=255)
              tttCATCCGGCCCCAACCGAGTCTGGCGGGGCTGCGGCGGTGGTGGGGTGGATTGTAGGCGCTGAGTGGCTGCCTCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  2:76413/150‑1 (MQ=255)
              gtgcatcaGGCACCAAAAGAGTCTGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGACGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGAAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  1:276286/147‑1 (MQ=255)
              attCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  2:1491823/149‑1 (MQ=255)
                tcatcaGGCACCAAACGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCGGTGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTa                                                                                                                                      <  1:1759669/150‑1 (MQ=255)
                  atcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCGAGGTACGGGTAATGGCAGGGGCGCATGCCCCTGCCTGGTGGATc                                                                                                                                    >  2:255319/1‑150 (MQ=255)
                  atcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTTCGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATc                                                                                                                                    >  2:1335290/1‑150 (MQ=255)
                     aGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGGGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAATCCCCAGCCTGGTGTATCTaa                                                                                                                                 >  1:1871336/1‑150 (MQ=255)
                      ggCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGTAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAg                                                                                                                                >  2:1459113/1‑150 (MQ=255)
                      ggCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAg                                                                                                                                <  1:1267532/150‑1 (MQ=255)
                       gCACCAACAGAGTCTGGCGGTGATGCGGCGTGGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGTAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAgg                                                                                                                               <  2:1114667/150‑1 (MQ=255)
                           cAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCCGGGGCGAATGCCCCTGCCGGGTGTATCTAAGGGCTa                                                                                                                           >  1:1458014/1‑150 (MQ=255)
                            aaatGAGGATGGCGGTGATGCGGCGTTGTTGGTGTGGATTGAAGTGGAAGAGTGGCTGCCGGCGGAGTGGCGCTGAGGAGGTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTAc                                                                                                                          <  2:647974/150‑1 (MQ=255)
                            aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTTACTGACTACGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTAc                                                                                                                          >  2:1067067/1‑150 (MQ=255)
                            aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCttc                                                                                                                          >  1:965792/1‑148 (MQ=255)
                               cgagTATGGCGGTGTTGCGGCGTTGTGGGGGTGGATGGAAGTCGACGCGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAAc                                                                                                                       <  2:508139/149‑1 (MQ=255)
                                gagTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACt                                                                                                                      <  1:1800539/150‑1 (MQ=255)
                                 agTATGGCGGTGATGGGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTc                                                                                                                     <  2:480125/150‑1 (MQ=255)
                                 agTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTc                                                                                                                     <  1:667534/150‑1 (MQ=255)
                                  gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCACCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  1:1129442/1‑150 (MQ=255)
                                   tttGGCGGTGATGCGGCGTTGTTGGTGTTGTTTGAACTCGAAGAGGGGCTGCCCCCTGAGTTGCCATGAGGAGTTGCGCGCGTTTCGCCAGGAACGGGAAATCGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAt                                                                                                                   <  1:880715/148‑1 (MQ=255)
                                       gCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                                                                                                               <  2:103462/150‑1 (MQ=255)
                                       gCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                                                                                                               >  2:1808817/1‑150 (MQ=255)
                                       gCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTAGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGATACAACTCATTCAg                                                                                                               >  2:1574157/1‑150 (MQ=255)
                                       gCGGTGATGCGGCGTTGGTGGTGTTGTTTGAAGGCGACGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                                                                                                               <  1:537384/150‑1 (MQ=255)
                                          gTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTAt                                                                                                            <  2:78863/150‑1 (MQ=255)
                                             aTGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTTTATAAGGGCTACAACTCCTTCAGTATGTa                                                                                                         >  1:1447099/1‑150 (MQ=255)
                                             aTGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTa                                                                                                         <  1:1304398/150‑1 (MQ=255)
                                              tGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCGGGTGTATCTATTGGCTACAACTCATTCTGGATggta                                                                                                        >  2:662320/1‑147 (MQ=255)
                                               gcggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCCCGTTGCGCCAGGTAAGGGTAATGGCAGGGGCGAAAGCCCCTGCCAGGTGTATCTAATGGCTACAACTCAATCAGTATGTaaa                                                                                                       >  2:1769935/1‑148 (MQ=255)
                                                cggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTg                                                                                                      <  2:431454/150‑1 (MQ=255)
                                                cggcgTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATATAAGGGCTACAACTCATTCAGTATGTACTg                                                                                                      >  2:1443394/1‑150 (MQ=255)
                                                     ttgttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAAcc                                                                                                 <  2:503300/150‑1 (MQ=255)
                                                        ttgGTGAGGATTGAAGTCGATGAGTGGCTGCCGCCGGAGTGGCCCTGAGGAGTTGCGCGCGGTGCGCCAGGGACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGc                                                                                              <  2:863831/150‑1 (MQ=255)
                                                         tgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCa                                                                                             >  2:1289804/1‑150 (MQ=255)
                                                         tgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCa                                                                                             >  1:1712822/1‑150 (MQ=255)
                                                          ggTGTTGATTGAAGTCGACGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCGGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCaa                                                                                            <  1:1195996/150‑1 (MQ=255)
                                                                gattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTTCGCCAGGTACGGGTAATGGCCGGGGCGAATGCCCCTGCCGGGTGTATCTAAGGGCTACAACTCATTCAGGATGTACAGGCACCGGCAAGTAGaa                                                                                      >  2:439874/1‑150 (MQ=255)
                                                                gattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGaa                                                                                      <  2:1748691/150‑1 (MQ=255)
                                                                gattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGaa                                                                                      <  1:1002596/150‑1 (MQ=255)
                                                                gattgaAGTAGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCCGGGGCGAATGCCCCTGCCTGGTGTATCTTAGGGATAAAACTAATTCAGTTTGAACTGGAACCCGCAATTAGaa                                                                                      >  2:1199023/1‑150 (MQ=255)
                                                                 attgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAAt                                                                                     >  2:1760368/1‑150 (MQ=255)
                                                                    gaAGTCTAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACCACTCATTCAGTATTTACTGGAAACGGCAAGTAGAATCGa                                                                                  >  2:624094/1‑150 (MQ=255)
                                                                    gaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGa                                                                                  >  2:179628/1‑150 (MQ=255)
                                                                    gaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGa                                                                                  >  1:1118772/1‑150 (MQ=255)
                                                                    gaAGTCGAAGAGTGGCAGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGGTTCTAAGGGCTACAACTCATTCAGTATGTACTGGCGCCGGCAAGTAGAATCGa                                                                                  >  1:197781/1‑150 (MQ=255)
                                                                       gTCGAAGAGTGGCTGCCGCCGGAGTTGCACTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTAATTCAGTATGTAGTGGATACGGCAAGTAGAATCGATTg                                                                               >  1:1362925/1‑150 (MQ=255)
                                                                       gTCGAAGAGGGGCTGCCGCCGGAGTTGCCCTGAGGAGCTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGGGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTg                                                                               <  2:1011546/150‑1 (MQ=255)
                                                                        tCGATGTGTGGCTGCCGCCGGAGTGGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTgg                                                                              <  1:1391677/150‑1 (MQ=255)
                                                                         cgtagagGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGTTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGc                                                                             <  1:172264/147‑1 (MQ=255)
                                                                         cGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGc                                                                             >  1:1029890/1‑150 (MQ=255)
                                                                          gAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCt                                                                            <  1:500696/150‑1 (MQ=255)
                                                                           aaGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTTAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTg                                                                           >  2:683551/1‑150 (MQ=255)
                                                                           aaGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTg                                                                           >  1:1480918/1‑150 (MQ=255)
                                                                            agcggggcggcCGGGGGCGTTGCCCTTAGGAGTTGAGCGCGTTGCGCCAGGTCCGGGTCATGGCCGGGGCGATTGCCGCTGCCTGGTGTTTCAAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGc                                                                          <  1:648056/145‑1 (MQ=255)
                                                                             gagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                         <  2:1049163/150‑1 (MQ=255)
                                                                              agTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCAt                                                                        <  2:506956/150‑1 (MQ=255)
                                                                                 ggCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATAct                                                                     >  1:409823/1‑150 (MQ=255)
                                                                                     gccgccGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgctg                                                                 <  2:30636/150‑1 (MQ=255)
                                                                                          gggAGTGGCCCTGAGGTGTGGTGCGCGTGGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATc                                                            <  1:968166/149‑1 (MQ=255)
                                                                                              gTTGCGCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGGAATGGCCGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAaca                                                        <  1:1180059/150‑1 (MQ=255)
                                                                                                   ccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcac                                                   >  1:1228634/1‑150 (MQ=255)
                                                                                                    ctgtcgagTTGCGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAGTGGCTGCATACTGGTGGTATCAAcccgcacg                                                  <  2:1280473/145‑1 (MQ=255)
                                                                                                    cTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacg                                                  >  2:1384848/1‑150 (MQ=255)
                                                                                                        ggagTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCa                                              <  1:864435/150‑1 (MQ=255)
                                                                                                        ggagTGGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCa                                              <  2:43239/150‑1 (MQ=255)
                                                                                                             tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCg                                         >  1:732630/1‑150 (MQ=255)
                                                                                                             ggAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCg                                         <  1:541101/149‑1 (MQ=255)
                                                                                                               tgcggGTTGCGCCAGGTACGGGTAATGGCAGGGGCTAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgtg                                       <  1:1769935/149‑1 (MQ=255)
                                                                                                                   ggTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAAc                                   <  1:1601006/149‑1 (MQ=255)
                                                                                                                       gcgcCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGtt                               <  1:1034578/150‑1 (MQ=255)
                                                                                                                        cgcCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGttt                              <  2:1313228/150‑1 (MQ=255)
                                                                                                                         gcCAGGTACGGGTAATGGCAGGGGCGAATTCCCCTGCCTGGTGTATATAAGGGCTACAACTCATTAAGTATGTACTGGAACCGGCAAGTAGAATAGATTTGCTGCATAATGCTGGTTTAAAAACGCACGACCAAATCGTGTAACAGtttt                             >  2:1590267/1‑149 (MQ=255)
                                                                                                                            aGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTcc                          <  1:1510826/150‑1 (MQ=255)
                                                                                                                              gTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTTCTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATAAACACGCACGACCATAACGTGTAAATGTTTGTCCGt                        >  1:200122/1‑150 (MQ=255)
                                                                                                                                aCGGGTAATGGCAGTGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTta                      <  1:479363/150‑1 (MQ=255)
                                                                                                                                   ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAg                   >  2:1652428/1‑150 (MQ=255)
                                                                                                                                    gTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                  <  2:1529927/150‑1 (MQ=255)
                                                                                                                                    gTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                  <  2:1665444/150‑1 (MQ=255)
                                                                                                                                         ggCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTc             >  2:1830524/1‑150 (MQ=255)
                                                                                                                                         ggCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATGGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATATGGTGTc             >  1:76363/1‑150 (MQ=255)
                                                                                                                                         ggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTc             >  2:645636/1‑150 (MQ=255)
                                                                                                                                           cAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACCCGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGa           >  2:154243/1‑150 (MQ=255)
                                                                                                                                             gtgtcGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGAtt         <  1:290153/146‑1 (MQ=255)
                                                                                                                                               ggCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTg       >  1:415270/1‑150 (MQ=255)
                                                                                                                                                 cGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc     <  1:34781/150‑1 (MQ=255)
                                                                                                                                                 cGAAGGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc     <  2:919174/150‑1 (MQ=255)
                                                                                                                                                    aTGCCCCTGCCTTGTGTATCTAAGGGCTAAAAATCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGATGCATACTGCTGGTATCAACACGCACGACCATATCTTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCTTg  <  2:1529660/150‑1 (MQ=255)
                                                                                                                                                    aTGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTg  <  1:1799890/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTG  >  NC_000913/2449054‑2449351

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