Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,653,958 T→G 6.3% intergenic (‑246/+50) oprB‑II ← / ← PP_1446 carbohydrate‑selective porin/TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,653,9580TG 6.3% 795.8 / 32.6 256intergenic (‑246/+50)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (123/105);  new base G (10/6);  total (145/111)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 264 total aligned reads displayed.

TTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAGG  >  NC_002947/1653809‑1654107
                                                                                                                                                     |                                                                                                                                                     
 ttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATc                                                                                                                                                      <  2:1521295/150‑1 (MQ=255)
    aGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTa                                                                                                                                                   >  1:4017237/1‑150 (MQ=255)
        tGCGGCGGTGAGTATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTcg                                                                                                                                               <  2:5968347/150‑1 (MQ=255)
          cggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTcgcg                                                                                                                                             <  1:6153761/150‑1 (MQ=255)
           ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                            <  1:3312425/150‑1 (MQ=255)
           ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                            <  1:6554304/150‑1 (MQ=255)
            gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAGTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCggg                                                                                                                                           <  2:2918312/150‑1 (MQ=21)
             cggGGAGGATATAGGGGAAGTTATAAGACATAAGTCAAATTGTCGCAGAGGATTGTTGTTGACCGGGCAACAGTCTTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGt                                                                                                                                          <  1:7752098/150‑1 (MQ=255)
              ggggAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTg                                                                                                                                         >  1:4325176/1‑150 (MQ=255)
               gggAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGa                                                                                                                                        >  2:8551776/1‑150 (MQ=255)
                ggaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCCCCGGTTGCCCAACCAACGTAATAGCGGGTGaa                                                                                                                                       >  1:5210624/1‑150 (MQ=255)
                ggaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGaa                                                                                                                                       <  2:64048/150‑1 (MQ=21)
                 gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                                      <  1:2752695/150‑1 (MQ=255)
                      tataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCGACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGtt                                                                                                                                 >  2:3032260/1‑150 (MQ=255)
                       ataGGGTAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                                <  2:6689373/150‑1 (MQ=255)
                       ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                                <  1:6371372/150‑1 (MQ=255)
                       agaGGGGAAGTTATAAGAAAAAAGACAAATTGACGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                                <  2:7368447/148‑1 (MQ=255)
                        taGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTcc                                                                                                                               <  1:6728014/150‑1 (MQ=255)
                        taGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTcc                                                                                                                               >  2:577242/1‑150 (MQ=255)
                           gggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTac                                                                                                                            >  1:5075033/1‑150 (MQ=255)
                            ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTaca                                                                                                                           <  1:2316240/150‑1 (MQ=255)
                            ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTAGTTGAACGGGCAACAGTCGTCATAAATCGAGGCAGGTGTAGGAGCGGGTTGACCCGCCAATGAGACCGTGGCCCGACCAAAGTATTAGCGTGTGAAACCCtacataca                                                                                                                           >  1:2323690/1‑150 (MQ=255)
                                     aaGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtc                                                                                                                  <  1:4354900/150‑1 (MQ=255)
                                     aaGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtc                                                                                                                  <  2:6336502/150‑1 (MQ=255)
                                      aGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCCTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcg                                                                                                                 >  2:210510/1‑150 (MQ=255)
                                       gAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgt                                                                                                                >  2:5230121/1‑150 (MQ=255)
                                           aaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCt                                                                                                            >  2:1919117/1‑150 (MQ=255)
                                                         gCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAgg                                                                                              >  1:1985392/1‑150 (MQ=255)
                                                          cAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggt                                                                                             >  2:1253302/1‑150 (MQ=255)
                                                           agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                            <  2:5597073/150‑1 (MQ=255)
                                                           agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                            <  2:1830446/150‑1 (MQ=255)
                                                           agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                            <  1:1549775/150‑1 (MQ=255)
                                                           agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCAACCCTCGTATTAGCGGGTGAACCCCGTCCAAAACCCAACTTCGTCGCTTGTGGTTCAGAAggtg                                                                                            >  1:6627954/1‑150 (MQ=38)
                                                             agcaTTGTTCTTGAACGGGCAACAGTCGTCATTAAACGAGGCAGGTGTAGGAGCGGTTTCACCCGCGAAAGAGCCCGTTGCCCTACCATCGAATACGCCGGTGAACCAGATCATACACCGAAATACGTCGCTTGTGGATCAGAAggagta                                                                                          >  2:3401008/4‑146 (MQ=255)
                                                             agGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCAACACCGATCTTCGTCGCTTGTGGTTCAGAAggtggt                                                                                          >  2:3400521/1‑150 (MQ=255)
                                                                   gttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCaa                                                                                    <  2:4400874/150‑1 (MQ=255)
                                                                    ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAc                                                                                   >  2:78046/1‑150 (MQ=255)
                                                                      gttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAccc                                                                                 <  1:5566661/150‑1 (MQ=255)
                                                                      gttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAccc                                                                                 <  2:6921338/150‑1 (MQ=255)
                                                                         gAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGAc                                                                              <  2:1790252/150‑1 (MQ=255)
                                                                          aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTAGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACt                                                                             >  1:7264115/1‑150 (MQ=255)
                                                                          aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCCCGAAAGCGCCGGGTGCCCTACCAAAGTATTCGCGGGTGAAACCGTTCCTAAACCGATCATCGTCGCATGTTGTTGAGAATGTGGTAAGCAccccccct                                                                             >  1:6896007/1‑146 (MQ=255)
                                                                          aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTAACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACt                                                                             >  1:2155106/1‑150 (MQ=255)
                                                                           aCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACtt                                                                            <  1:4880920/150‑1 (MQ=255)
                                                                             gggCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTGACCCCCGAATGCGCCGGTTGGCCAACCAACGTATTAGCGGGTGAACCCCTTACTAAAACCATCATCGTAGCTTGTTGTTCAGAAGGTGGTGCGCAACCCCAcatac                                                                          >  1:4463063/1‑146 (MQ=255)
                                                                                        gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCGGGTGCCCTAACAAAGTATTCGCGGGTGGACCCGTTCCAACACCCATCTTAGTCGCGTGTGGTTCAGAAGGGGGTAAGCAACCCCACATACACGCGACGccc                                                               >  1:8279312/1‑150 (MQ=255)
                                                                                         tcATTACTCGCCGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                                              <  1:4593472/150‑1 (MQ=255)
                                                                                         tcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                                              >  2:8706153/1‑150 (MQ=255)
                                                                                          cATTAATCGAGGCAGGTGTAGGAGCGGTTTCACCCGCGAATGCGCCGGTTGCACTACCATCGTATTCGCGGTTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCg                                                             <  1:8476181/150‑1 (MQ=255)
                                                                                          cATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCg                                                             <  2:795384/150‑1 (MQ=255)
                                                                                           aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATTGTATTCGCGGGTGAACCCGTTCCTAAACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCGACGCCCCgg                                                            >  1:3064382/1‑150 (MQ=255)
                                                                                           aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                            >  2:6581110/1‑150 (MQ=255)
                                                                                           aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                            >  2:3890609/1‑150 (MQ=255)
                                                                                                tCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATTCGCCGGTTGCCCTACCATCGTATTAGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGGGTGTGGTTGAGAAGGTGGTACGCAACCCGACTTCCACGCTAAGCCCCGGCGCCg                                                       >  1:3421956/1‑150 (MQ=39)
                                                                                                tCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCg                                                       >  1:7948354/1‑150 (MQ=255)
                                                                                                    gatcggTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGc                                                   >  2:2472288/5‑150 (MQ=17)
                                                                                                       aGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGat                                                >  2:829781/1‑150 (MQ=255)
                                                                                                          tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCTCGCTACGCCCCGGCGCCGGTGCGATATc                                             <  2:600624/150‑1 (MQ=255)
                                                                                                          tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATc                                             <  2:8111300/150‑1 (MQ=255)
                                                                                                          tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATc                                             >  2:3994586/1‑150 (MQ=255)
                                                                                                            tAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATcgc                                           >  2:6816893/1‑150 (MQ=255)
                                                                                                            tAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATcgc                                           >  1:8410283/1‑150 (MQ=18)
                                                                                                              ggAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCgcgc                                         >  1:1080013/1‑150 (MQ=18)
                                                                                                               gaacGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                        >  2:5117651/4‑150 (MQ=21)
                                                                                                               gAGCGGGTTCACCCGCGCATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                        <  1:32580/150‑1 (MQ=255)
                                                                                                               gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                        >  2:2652730/1‑150 (MQ=255)
                                                                                                               gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                        <  1:3704755/150‑1 (MQ=255)
                                                                                                                   gggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGAt                                    <  1:1298272/150‑1 (MQ=255)
                                                                                                                    ggTTCACCCGCGACTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                                   <  1:1410574/150‑1 (MQ=255)
                                                                                                                         aCCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACt                              >  2:3742616/1‑150 (MQ=255)
                                                                                                                          cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg                             >  1:2473405/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:3865504/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:4747184/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:4955094/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:546494/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:5816688/1‑150 (MQ=255)
                                                                                                                            cgcgAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGc                           >  2:2826695/1‑150 (MQ=255)
                                                                                                                             gcgAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCa                          >  2:5097091/1‑150 (MQ=255)
                                                                                                                                  tGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTCCGCTACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGCTCGAACTGGCATcgcg                     >  2:8816888/1‑150 (MQ=255)
                                                                                                                                   gcgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAgcgc                    >  1:2873977/1‑150 (MQ=255)
                                                                                                                                    cgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCa                   <  2:6719898/150‑1 (MQ=255)
                                                                                                                                     gcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCat                  >  1:1434480/1‑149 (MQ=255)
                                                                                                                                         gTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTc              <  2:1441955/150‑1 (MQ=255)
                                                                                                                                             cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggt          <  2:1320348/150‑1 (MQ=255)
                                                                                                                                             cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggt          >  2:4113147/1‑150 (MQ=255)
                                                                                                                                              ccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtc         >  2:3080699/1‑150 (MQ=255)
                                                                                                                                              ccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtc         >  1:5186738/1‑150 (MQ=255)
                                                                                                                                               cTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGGTCAGAAGGTGGTACGCAACCCGACTTCCCCGCTACGCCCCGGCGCCGGTGCGATATCGCGCCGGATCGAACTGGCATAGCGCACGGTCTggtcg        >  2:3304969/1‑150 (MQ=255)
                                                                                                                                                tACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcgg       >  2:675046/1‑150 (MQ=255)
                                                                                                                                                tACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcgg       >  2:4309073/1‑150 (MQ=255)
                                                                                                                                                  ccATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcggtc     <  2:5538993/150‑1 (MQ=255)
                                                                                                                                                   cATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCa    >  2:2365576/1‑150 (MQ=255)
                                                                                                                                                   cATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCa    <  1:1489944/150‑1 (MQ=255)
                                                                                                                                                    aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAg   <  1:3075123/150‑1 (MQ=255)
                                                                                                                                                    aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAg   <  1:6469131/150‑1 (MQ=255)
                                                                                                                                                     tCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAgg  >  2:8001641/1‑150 (MQ=255)
                                                                                                                                                     tCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAgg  >  1:4353530/1‑150 (MQ=255)
                                                                                                                                                     tCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAgg  >  1:8733973/1‑150 (MQ=255)
                                                                                                                                                     tCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAgg  >  2:618122/1‑150 (MQ=255)
                                                                                                                                                     tCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAgg  >  1:6942302/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAGG  >  NC_002947/1653809‑1654107

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