Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,653,951 C→G 6.7% intergenic (‑239/+57) 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,9510CG 6.7% 790.9 / 34.9 239intergenic (‑239/+57)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base C (124/97);  new base G (10/6);  total (134/105)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 256 total aligned reads displayed.

TGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTC  >  NC_002947/1653802‑1654100
                                                                                                                                                     |                                                                                                                                                     
tGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAACAAGACACATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGcc                                                                                                                                                       <  2:1272992/150‑1 (MQ=255)
tGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGcc                                                                                                                                                       <  1:1470062/150‑1 (MQ=255)
tGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGcc                                                                                                                                                       <  1:4871115/150‑1 (MQ=255)
  aaGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCt                                                                                                                                                     <  2:5392589/150‑1 (MQ=255)
  aaGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCt                                                                                                                                                     >  1:8254636/1‑150 (MQ=255)
   aGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTTCCCTa                                                                                                                                                    >  1:6648313/1‑150 (MQ=255)
       tttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAt                                                                                                                                                >  1:4376071/1‑150 (MQ=255)
       tttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAAGGCGCCGGTTGCCCTACCAt                                                                                                                                                >  1:5352473/1‑150 (MQ=255)
        ttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGACGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATc                                                                                                                                               <  2:3265843/150‑1 (MQ=255)
        ttCAGCTTGCGGCGGGGAGGATAGAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATc                                                                                                                                               <  2:4484813/150‑1 (MQ=255)
        ttCAGCGTGCGGCGGGGAGGAGAGAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTTTTAAACGGTCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCTAATGCGCCGGTTGCCCTACCATc                                                                                                                                               <  2:7758310/150‑1 (MQ=255)
             cTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                                          >  1:489537/1‑150 (MQ=255)
                 cggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTcgcg                                                                                                                                      <  1:6153761/150‑1 (MQ=255)
                  ggcggGGAGGATATCGGGGAAGTGATAAGAACAAAGACAACTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                     <  1:7973162/150‑1 (MQ=21)
                  ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                     <  1:7708142/150‑1 (MQ=255)
                  ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                     <  2:7977837/150‑1 (MQ=255)
                   gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAGTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCggg                                                                                                                                    <  2:2918312/150‑1 (MQ=21)
                     ggggAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTg                                                                                                                                  >  1:4325176/1‑150 (MQ=255)
                      gggAGGATAGAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGa                                                                                                                                 <  1:2959902/150‑1 (MQ=255)
                       ggaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCCCCGGTTGCCCAACCAACGTAATAGCGGGTGaa                                                                                                                                >  1:5210624/1‑150 (MQ=255)
                       ggaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGaa                                                                                                                                <  2:64048/150‑1 (MQ=21)
                             tataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCCGTTTCCCTACCAACGTATTCGCGGGTGAAacccat                                                                                                                          >  1:412601/1‑147 (MQ=255)
                              atcggggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                         <  1:5010368/147‑1 (MQ=255)
                              ataGGGTAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                         <  2:6689373/150‑1 (MQ=255)
                              ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                         <  1:6371372/150‑1 (MQ=255)
                               taGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTcc                                                                                                                        <  1:6728014/150‑1 (MQ=255)
                                 ggggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCAACCAAAGGATTCGCGGGTGAACCCCTTCCTa                                                                                                                      >  1:6711003/1‑150 (MQ=255)
                                  gggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAAAGCGCCGCTTGACCTACCAACGTATTCGCGGGTGAACCCGTTCCTac                                                                                                                     >  1:5075038/1‑150 (MQ=255)
                                   ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTAGTTGAACGGGCAACAGTCGTCATAAATCGAGGCAGGTGTAGGAGCGGGTTGACCCGCCAATGAGACCGTGGCCCGACCAAAGTATTAGCGTGTGAAACCCtacataca                                                                                                                    >  1:2323690/1‑150 (MQ=255)
                                            aaGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtc                                                                                                           <  2:6336502/150‑1 (MQ=255)
                                            aaGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtc                                                                                                           <  1:4353785/150‑1 (MQ=255)
                                              gAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgt                                                                                                         >  2:5999876/1‑150 (MQ=255)
                                                 aaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGc                                                                                                      <  2:713757/150‑1 (MQ=255)
                                                 aaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGc                                                                                                      <  1:8812246/150‑1 (MQ=255)
                                                      aCAAATTGGCGCAGAGGATTGTTGTTGAACGGGCAACTGTCGTCATTAATCGAGGCAGGTTTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCACTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGTTCTTCGTCGCTtgtg                                                                                                 <  1:582986/150‑1 (MQ=39)
                                                        aaaTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTAACCCGCGAATGCGCCGGTTGGCCTACCATCGTATTAGCGGGTGTACCCGTTGCTACACCGATCTTAGTGGCTTGTGGt                                                                                               >  1:192901/1‑150 (MQ=255)
                                                             gTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGa                                                                                          <  2:5676949/150‑1 (MQ=255)
                                                                 cAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTGACCCGCGAATGGGCCGGTTGGCCTACCATCGTATTAGCGGGTGAAACCGTTGCTACAACGATCTTCGTGGCTTGTGGTTGAGAAggt                                                                                      >  1:8321222/1‑150 (MQ=255)
                                                                  agagGATTGTTGTTGAACGGGCAACGGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                     <  2:5232352/150‑1 (MQ=255)
                                                                  agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                     <  2:5597073/150‑1 (MQ=255)
                                                                  agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                     <  2:1830446/150‑1 (MQ=255)
                                                                    agGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCAACACCGATCTTCGTCGCTTGTGGTTCAGAAggtggt                                                                                   >  2:3400521/1‑150 (MQ=255)
                                                                           ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCAACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAc                                                                            >  2:1537140/1‑150 (MQ=255)
                                                                            tgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAcc                                                                           <  1:6555661/150‑1 (MQ=255)
                                                                                gAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGAACTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGAc                                                                       >  1:2820489/1‑150 (MQ=255)
                                                                                 aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTAGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACt                                                                      >  1:7264115/1‑150 (MQ=255)
                                                                                 aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCCCGAAAGCGCCGGGTGCCCTACCAAAGTATTCGCGGGTGAAACCGTTCCTAAACCGATCATCGTCGCATGTTGTTGAGAATGTGGTAAGCAccccccct                                                                      >  1:6896007/1‑146 (MQ=255)
                                                                                  aCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATTCGCCCGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACtt                                                                     >  1:5499367/1‑150 (MQ=255)
                                                                                    gggCAACAGTCGTCATTCATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTcc                                                                   >  2:220178/1‑150 (MQ=255)
                                                                                    gggCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGCACGCAACCCGACTTcc                                                                   >  2:2613702/1‑150 (MQ=255)
                                                                                      gCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATTCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTAGTCGCTTGTTGTTCAGAAGGTGGTAAGCAAACCGACTTCCAc                                                                 >  2:6177779/1‑150 (MQ=39)
                                                                                      gCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCAc                                                                 >  1:2038021/1‑150 (MQ=255)
                                                                                          cAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTa                                                             <  1:2734012/150‑1 (MQ=255)
                                                                                               gtcATTAATCGAGGCAGTTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGccc                                                        <  1:2155483/150‑1 (MQ=255)
                                                                                               gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGccc                                                        <  1:2362329/150‑1 (MQ=255)
                                                                                                 cATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCg                                                      <  2:7946792/150‑1 (MQ=255)
                                                                                                  aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGTAACCCGACTTCCACGCTACGCCCCgg                                                     >  2:1513547/1‑150 (MQ=255)
                                                                                                  aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                     >  2:3890609/1‑150 (MQ=255)
                                                                                                  aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                     >  2:6581110/1‑150 (MQ=255)
                                                                                                  aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTAACATTGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACATCCACGCTACGCCCCgg                                                     >  1:3064839/1‑150 (MQ=255)
                                                                                                      aTCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGcc                                                 >  2:6888742/1‑150 (MQ=255)
                                                                                                      aTCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGcc                                                 >  2:4257446/1‑150 (MQ=255)
                                                                                                       tCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATTCGCCGGTTGCCCTACCATCGTATTAGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGGGTGTGGTTGAGAAGGTGGTACGCAACCCGACTTCCACGCTAAGCCCCGGCGCCg                                                >  1:3421956/1‑150 (MQ=39)
                                                                                                       tCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCg                                                >  1:7948354/1‑150 (MQ=255)
                                                                                                           gatcggTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTAGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGc                                            >  2:2471264/5‑150 (MQ=17)
                                                                                                            gCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCGTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCg                                           >  2:2700091/1‑150 (MQ=255)
                                                                                                              aGGTGTAGGAGCGGGTTCACCCGCGACTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCACCGGCGCCGGGGCgta                                         <  2:1326275/150‑3 (MQ=255)
                                                                                                              aGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCgta                                         <  1:8651420/150‑3 (MQ=255)
                                                                                                                 tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGTAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATc                                      <  1:1513547/150‑1 (MQ=255)
                                                                                                                 tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATc                                      <  2:8111300/150‑1 (MQ=255)
                                                                                                                 tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATc                                      >  2:3994586/1‑150 (MQ=255)
                                                                                                                     ggAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCgcgc                                  >  1:1080013/1‑150 (MQ=18)
                                                                                                                      gaacGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                 >  2:5117651/4‑150 (MQ=21)
                                                                                                                      gAGCGGGTTCACCCGCGCATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                 <  1:32580/150‑1 (MQ=255)
                                                                                                                      gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                 <  1:3704755/150‑1 (MQ=255)
                                                                                                                      gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                 <  1:32634/150‑1 (MQ=255)
                                                                                                                      gAGCGGGTTCACCCGCGAATGCACCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                 <  2:876587/150‑1 (MQ=255)
                                                                                                                      gAGCGGGTTCACCCGCGAATGCACCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                 <  2:3437678/150‑1 (MQ=255)
                                                                                                                        gCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAgg                               >  2:2995908/1‑150 (MQ=255)
                                                                                                                          gggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGAt                             <  1:1298272/150‑1 (MQ=255)
                                                                                                                           ggTTCACCCGCGACTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                            <  1:1410574/150‑1 (MQ=255)
                                                                                                                           ggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGAGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                            >  2:6726510/1‑150 (MQ=255)
                                                                                                                            gTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCg                           >  2:6790907/1‑150 (MQ=255)
                                                                                                                             ttCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGa                          >  1:8465548/1‑150 (MQ=255)
                                                                                                                               cACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAAc                        >  1:4492685/1‑150 (MQ=255)
                                                                                                                               cACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGAt                             <  2:8395869/145‑1 (MQ=25)
                                                                                                                                aCCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACt                       >  2:8598764/1‑150 (MQ=255)
                                                                                                                                aCCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACt                       >  2:3742616/1‑150 (MQ=255)
                                                                                                                                 cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg                      >  1:4817863/1‑150 (MQ=255)
                                                                                                                                  ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                     >  1:3865504/1‑150 (MQ=255)
                                                                                                                                   cgcgAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGc                    >  1:4604801/1‑150 (MQ=255)
                                                                                                                                   cgcgAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCGTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGCTCGAACTGGc                    >  2:1812157/1‑150 (MQ=255)
                                                                                                                                           cgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCa            >  2:2054333/1‑150 (MQ=255)
                                                                                                                                            gcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCat           >  1:1434480/1‑149 (MQ=255)
                                                                                                                                              cGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACgg         <  2:5940099/150‑1 (MQ=255)
                                                                                                                                              cGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACgg         >  1:6480034/1‑150 (MQ=255)
                                                                                                                                               ggTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGt        <  2:8254636/150‑1 (MQ=255)
                                                                                                                                                    cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggt   >  2:5113887/1‑150 (MQ=255)
                                                                                                                                                    cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggt   >  2:4113147/1‑150 (MQ=255)
                                                                                                                                                     ccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtc  >  2:3080699/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTC  >  NC_002947/1653802‑1654100

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: