Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 997,726 C→T 5.8% intergenic (‑128/+104) PP_0860 ← / ← PP_0861 sulfite reductase flavoprotein subunit/outer membrane ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947997,7260CT 5.8% 304.7 / 6.5 103intergenic (‑128/+104)PP_0860/PP_0861sulfite reductase flavoprotein subunit/outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (48/48);  new base T (4/2);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 107 total aligned reads displayed.

TTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCG  >  NC_002947/997578‑997873
                                                                                                                                                    |                                                                                                                                                   
ttGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCctct                                                                                                                                                    >  1:1149516/1‑150 (MQ=255)
 tGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCAGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACTTTAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTGTGTGCGAGCCTCtt                                                                                                                                                   <  2:860083/150‑1 (MQ=255)
    aaCAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTcgc                                                                                                                                                >  1:959121/1‑150 (MQ=255)
      cAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCgg                                                                                                                                              <  1:1515975/150‑1 (MQ=255)
      cAGCGTCTTCTTCACCACGGTTCTCACCTGGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCGCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCgg                                                                                                                                              <  2:1856810/150‑1 (MQ=255)
              tcttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGCCTGTGCCGGCCCTTTCGCGGc                                                                                                                                             >  2:1166739/1‑143 (MQ=255)
              tcttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGCCTGTGCCGGCCCTTTCGCGGc                                                                                                                                             <  1:1166739/143‑1 (MQ=255)
              tcttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGCCTGTGCCGGCCCTTTCGCGGCTGAAGcc                                                                                                                                      >  1:1285630/1‑150 (MQ=255)
               cttcACCACGTTTCTCACCTTGCGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCAGGGAGGCCTGCTTTGCGTGTTGTCTGTGCGGGCCTATTCGCGGCTAAAGccc                                                                                                                                     >  2:269572/1‑149 (MQ=255)
                ttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGc                                                                                                                                    <  1:1612728/150‑1 (MQ=255)
                 tcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCt                                                                                                                                   >  1:1670970/1‑150 (MQ=255)
                 tcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCt                                                                                                                                   >  1:996046/1‑150 (MQ=255)
                   accacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTcc                                                                                                                                 <  2:637618/150‑1 (MQ=255)
                    ccacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTATTCGTGTTGTCTGTCCGCCCCTCTCGCGGCTGTAGCCGCCCCTg                                                                                                                               >  2:419392/1‑150 (MQ=255)
                     cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTg                                                                                                                               <  2:1525265/150‑1 (MQ=255)
                     cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTg                                                                                                                               >  1:609384/1‑150 (MQ=255)
                     cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTg                                                                                                                               >  1:1525265/1‑150 (MQ=255)
                        gTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCAc                                                                                                                            <  1:1524996/150‑1 (MQ=255)
                           tctcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAATGTGCTTCTGCACCgg                                                                                                                         <  2:959121/150‑1 (MQ=255)
                             tcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAg                                                                                                                                        <  1:1428622/133‑1 (MQ=255)
                             tcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAg                                                                                                                                        >  2:1428622/1‑133 (MQ=255)
                             tcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCCGGCCCTTTCGCGGCTGAAGCCGCTCCTGcgttgccg                                                                                                                       >  1:1565908/1‑143 (MQ=38)
                               aCCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCAcc                                                                                                                           >  1:1601945/1‑144 (MQ=255)
                               aCCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCAcc                                                                                                                           <  2:1601945/144‑1 (MQ=255)
                                ccTTGAGTTGCTGCCGAGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGAAGGGCTGCTTCACGGATCGCCGCTGCAACGTGCAGCGCGTCTAAAACCGCTCCTGCACCGGCAACg                                                                                                                    >  2:909332/1‑150 (MQ=255)
                                  ttGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCa                                                                                                                  >  1:1541118/1‑150 (MQ=255)
                                       ttGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGc                                                                                                             >  2:1312372/1‑150 (MQ=255)
                                        tgctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCg                                                                                                            <  1:269572/150‑1 (MQ=255)
                                          ctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCggg                                                                                                          >  2:327603/1‑149 (MQ=255)
                                          ctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGc                                                                                                          <  2:1304787/150‑1 (MQ=255)
                                          ctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGc                                                                                                          <  1:1104083/150‑1 (MQ=255)
                                          ctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGc                                                                                                          <  1:310521/150‑1 (MQ=255)
                                               gcgACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAg                                                                                                     <  1:1835712/150‑1 (MQ=255)
                                                    ccGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCCACGAAGGAGCGGCTTCAg                                                                                                     >  1:99619/1‑145 (MQ=255)
                                                     cgcgcgGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGAGGCGGGAGGGCTGCTTTTCGTGTTGCCTGTGCCGGCCCTTTCGCGGCTGAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGa                                                                                               >  1:348298/1‑150 (MQ=17)
                                                       cgcgGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGaaa                                                                                             <  2:1741110/150‑1 (MQ=255)
                                                        gcggcgTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAg                                                                                            <  1:440326/150‑1 (MQ=255)
                                                         cggcgTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGAGAAAgg                                                                                           <  2:1541118/150‑1 (MQ=255)
                                                            cgTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGcc                                                                                        >  1:1361715/1‑150 (MQ=255)
                                                             gTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTGAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGaagaggccc                                                                                       <  1:1346957/150‑8 (MQ=18)
                                                             gTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCg                                                                                       <  2:1603812/150‑1 (MQ=255)
                                                                  ttATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGcaca                                                                                  >  2:1494451/1‑150 (MQ=255)
                                                                    aTGCCTTGATACACAAAAGCCCCTCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCCTTTCGCGGCTGAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCCGCACAga                                                                                <  1:1273772/150‑2 (MQ=9)
                                                                        cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGc                                                                                         <  1:767872/137‑1 (MQ=255)
                                                                        cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGc                                                                                         >  2:767872/1‑137 (MQ=255)
                                                                        cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCa                                                                            >  1:902288/1‑150 (MQ=255)
                                                                        cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCa                                                                            >  1:1462527/1‑150 (MQ=255)
                                                                         ttGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAg                                                                           <  2:1906106/150‑1 (MQ=255)
                                                                         ttGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAg                                                                           >  2:1762687/1‑150 (MQ=255)
                                                                             tACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCt                                                                       >  2:1285262/1‑150 (MQ=255)
                                                                                  aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAgg                                                                                           <  2:1801574/125‑1 (MQ=255)
                                                                                  aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAgg                                                                                >  1:1530064/1‑136 (MQ=255)
                                                                                  aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAgg                                                                                <  2:1530064/136‑1 (MQ=255)
                                                                                   aaaGCCCCGGGCGCGGGAAGGAGCGGGGCCGGGGGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCGCGTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACa                                                                 <  2:807369/150‑1 (MQ=255)
                                                                                       ccccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCcgc                                                                                                 <  2:560996/114‑1 (MQ=255)
                                                                                       ccccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCcgc                                                                                                 >  1:560996/1‑114 (MQ=255)
                                                                                         ccGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaa                                                           >  1:550041/1‑150 (MQ=255)
                                                                                          cgcgCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaac                                                          <  2:474303/150‑1 (MQ=255)
                                                                                          cGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCNTCAGCCNCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaac                                                          <  1:1312372/150‑1 (MQ=255)
                                                                                               aCGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATcaacccaac                                                     >  1:277521/1‑150 (MQ=255)
                                                                                                 gggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCaacccaaccc                                                   >  1:1539833/1‑150 (MQ=255)
                                                                                                 gggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCaacccaaccc                                                   >  2:329304/1‑150 (MQ=255)
                                                                                                  ggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGTCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCt                                                  >  2:1219632/1‑150 (MQ=255)
                                                                                                   gAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCgg                                                                                                                                              <  2:755196/57‑1 (MQ=255)
                                                                                                    aaTGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGa                                                <  2:1462527/150‑1 (MQ=255)
                                                                                                        aaaggggCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGAGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCTGCGAAAGGGCCGGCACAGGCGCAGCCCTTGAACAGATCAACCCAACCCTTACagc                                            >  2:84295/4‑148 (MQ=255)
                                                                                                         aCGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCa                                           >  2:1227933/1‑150 (MQ=255)
                                                                                                            gggCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGTCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGaa                                        <  1:1219632/150‑1 (MQ=255)
                                                                                                            gggCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGaa                                        <  2:550041/150‑1 (MQ=255)
                                                                                                            gggCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGaa                                        <  1:329304/150‑1 (MQ=255)
                                                                                                              gCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGt                                      <  2:609384/150‑1 (MQ=255)
                                                                                                              gCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGt                                      <  1:556015/150‑1 (MQ=255)
                                                                                                                  ggAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCAc                                  >  2:1103622/1‑150 (MQ=255)
                                                                                                                  ggAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCAc                                  >  2:694449/1‑150 (MQ=255)
                                                                                                                   gAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACg                                 <  2:1149516/150‑1 (MQ=255)
                                                                                                                      ggCTGCTTTTCGTGTTGTCTGTCCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtg                              >  1:551905/1‑150 (MQ=255)
                                                                                                                       gctgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCGCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCGGGCACAGGCGCAGCCCATCAACAGATCAACCCAACCCTGACATCATGATTGCACGCtgg                             <  2:1249615/150‑1 (MQ=255)
                                                                                                                        ctgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGCGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggt                            >  1:686349/1‑150 (MQ=255)
                                                                                                                         tgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           <  1:1103622/150‑1 (MQ=255)
                                                                                                                         tgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           <  2:996046/150‑1 (MQ=255)
                                                                                                                         tgctTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           >  1:1672942/1‑150 (MQ=255)
                                                                                                                            ttttCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTggtggta                        >  2:53435/1‑149 (MQ=255)
                                                                                                                                 gtgtTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAg                   >  2:1708265/1‑150 (MQ=255)
                                                                                                                                  tgttgtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGt                  >  1:1663463/1‑150 (MQ=255)
                                                                                                                                     tgtCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGcc               <  2:1361715/150‑1 (MQ=255)
                                                                                                                                       tCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGt             >  2:714075/1‑150 (MQ=255)
                                                                                                                                         tgtgCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGc                                                                                                    >  2:1549688/1‑61 (MQ=255)
                                                                                                                                         tgtgCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGc                                                                                                    <  1:1549688/61‑1 (MQ=255)
                                                                                                                                         tgtgCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTCCAGCACCGAAAGGCATGGTACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTgc           <  1:1689750/150‑1 (MQ=255)
                                                                                                                                          gtgCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTgcg          >  1:1881782/1‑150 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCGCGCTGGTGGTCAACAGGGCGGTGCGCCCGgc    >  1:1248064/1‑150 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTg                 >  2:1113309/1‑137 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTg                 >  2:1113277/1‑137 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTg                 <  1:1113309/137‑1 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgc    >  2:1589258/1‑150 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgc    <  2:1670970/150‑1 (MQ=255)
                                                                                                                                                gCCTCTTCGCGGCTAAAGCCGCTCCTGCAACGGCAACGCAGGGGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTg                 <  1:1113277/137‑1 (MQ=255)
                                                                                                                                                  ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           <  1:781937/125‑1 (MQ=255)
                                                                                                                                                  ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           >  2:781937/1‑125 (MQ=255)
                                                                                                                                                  ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccg  >  2:1872060/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCG  >  NC_002947/997578‑997873

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