Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 997,725 T→C 5.9% intergenic (‑127/+105) 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,7250TC 5.9% 298.2 / 6.3 103intergenic (‑127/+105)PP_0860/PP_0861sulfite reductase flavoprotein subunit/outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (48/47);  new base C (4/2);  total (53/50)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

Only 100 of 108 total aligned reads displayed.

AATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCG  >  NC_002947/997576‑997873
                                                                                                                                                     |                                                                                                                                                    
aaTTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCct                                                                                                                                                      >  2:1524996/1‑150 (MQ=255)
  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)
                tcttcACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGGGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGcc                                                                                                                                      >  2:1689750/1‑150 (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)
                      ccacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCt                                                                                                                                >  2:556015/1‑150 (MQ=255)
                      ccacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTATTCGTGTTGTCTGTCCGCCCCTCTCGCGGCTGTAGCCGCCCCTg                                                                                                                               >  2:419392/1‑150 (MQ=255)
                       cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTg                                                                                                                               >  1:1525265/1‑150 (MQ=255)
                       cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTg                                                                                                                               >  1:609384/1‑150 (MQ=255)
                       cacGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTg                                                                                                                               <  2:1525265/150‑1 (MQ=255)
                          gTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCAc                                                                                                                            <  1:1524996/150‑1 (MQ=255)
                             tctcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAATGTGCTTCTGCACCgg                                                                                                                         <  2:959121/150‑1 (MQ=255)
                               tcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAg                                                                                                                                        >  2:1428622/1‑133 (MQ=255)
                               tcACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCCGGCCCTTTCGCGGCTGAAGCCGCTCCTGcgttgccg                                                                                                                       >  1:1565908/1‑143 (MQ=38)
                                 aCCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCAcc                                                                                                                           <  2:1601945/144‑1 (MQ=255)
                                 aCCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCAcc                                                                                                                           >  1:1601945/1‑144 (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)
                                            ctgcCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCAGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGc                                                                                                          <  2:1303629/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)
                                                                      aTGCCTTGATACACAAAAGCCCCTCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCCTTTCGCGGCTGAAGCCGCTCCTGCGTTGCCGGTGCAGGAGCGGCTTTAGCCGCGAAGAGGCCCGCACAga                                                                                <  1:1273772/150‑2 (MQ=9)
                                                                          cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGc                                                                                         >  2:767872/1‑137 (MQ=255)
                                                                          cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGc                                                                                         <  1:767872/137‑1 (MQ=255)
                                                                          cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCa                                                                            >  1:902288/1‑150 (MQ=255)
                                                                          cTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCa                                                                            >  1:1462527/1‑150 (MQ=255)
                                                                           ttGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAg                                                                           <  2:1906106/150‑1 (MQ=255)
                                                                               tACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCt                                                                       >  2:1285262/1‑150 (MQ=255)
                                                                                    aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAgg                                                                                           >  1:1801574/1‑125 (MQ=255)
                                                                                    aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAgg                                                                                           <  2:1801574/125‑1 (MQ=255)
                                                                                    aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAgg                                                                                >  1:1530064/1‑136 (MQ=255)
                                                                                    aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAgg                                                                                <  2:1530064/136‑1 (MQ=255)
                                                                                    aaaaGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAAc                                                                  <  1:1494451/150‑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                                                   >  2:329304/1‑150 (MQ=255)
                                                                                                   gggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCaacccaaccc                                                   >  1:1539833/1‑150 (MQ=255)
                                                                                                    ggAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGTCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCt                                                  >  2:1219632/1‑150 (MQ=255)
                                                                                                     gAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCgg                                                                                                                                              <  2:755196/57‑1 (MQ=255)
                                                                                                     gAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCgg                                                                                                                                              >  1:755196/1‑57 (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                                      <  1:556015/150‑1 (MQ=255)
                                                                                                                    ggAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCAc                                  >  2:694449/1‑150 (MQ=255)
                                                                                                                    ggAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCAc                                  >  2:1103622/1‑150 (MQ=255)
                                                                                                                     gAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACg                                 <  2:1149516/150‑1 (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)
                                                                                                                                       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                 <  1:1113309/137‑1 (MQ=255)
                                                                                                                                                  gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgc    <  2:1670970/150‑1 (MQ=255)
                                                                                                                                                  gCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgc    >  2:1589258/1‑150 (MQ=255)
                                                                                                                                                  gCCTCTTCGCGGCTAAAGCCGCTCCTGCAACGGCAACGCAGGGGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTg                 <  1:1113277/137‑1 (MQ=255)
                                                                                                                                                    ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           >  2:781937/1‑125 (MQ=255)
                                                                                                                                                    ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCtggtg                           <  1:781937/125‑1 (MQ=255)
                                                                                                                                                    ctctTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGgccg  >  2:1872060/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTGGAACAGCGTCTTCTTCACCACGTTTCTCACCTTGAGTTGCTGCCGCGACCCGCGCGGCGTGCATTATGCCTTGATACACAAAAGCCCCGCTCACGGGAATGAACGGGGCCGGGAGGGCTGCTTTTCGTGTTGTCTGTGCGGGCCTCTTCGCGGCTAAAGCCGCTCCTGCACCGGCAACGCAGGAGCGGCTTCAGCCGCGAAAGGGCCGGCACAGGCGCAGCCCTTCAACAGATCAACCCAACCCTGACATCAGAAGTGCACGCTGGTGGTCAACAGTGCCGTGCGCCCGGCCG  >  NC_002947/997576‑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: