Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,653,953 T→A 9.2% intergenic (‑241/+55) 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,9530TA 9.2% 743.5 / 36.2 248intergenic (‑241/+55)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (122/99);  new base A (19/6);  total (143/105)
Fisher's exact test for biased strand distribution p-value = 5.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

AAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGG  >  NC_002947/1653804‑1654102
                                                                                                                                                     |                                                                                                                                                      
aaGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCt                                                                                                                                                        <  2:5392067/150‑1 (MQ=255)
aaGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCt                                                                                                                                                        >  1:8254636/1‑150 (MQ=255)
     tttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGACTGCGCCGGTTGCCCTACCag                                                                                                                                                   >  1:5352583/1‑149 (MQ=255)
     tttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAt                                                                                                                                                   >  1:4376071/1‑150 (MQ=255)
      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)
                ggcggGGAGGATATCGGGGAAGTGATAAGAACAAAGACAACTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                        <  1:7973162/150‑1 (MQ=21)
                ggcggGGAGGATATAGGGGAAGTTATAAGAAAATAGACAAATTGTCGCAGAGGATTGTTGTTAAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                        <  2:1344585/150‑1 (MQ=255)
                ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                        <  1:7708142/150‑1 (MQ=255)
                ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                        <  2:1580898/150‑1 (MQ=255)
                ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                        <  1:3312425/150‑1 (MQ=255)
                ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                        <  2:655272/150‑1 (MQ=25)
                 gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGAGTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCggg                                                                                                                                       <  2:2918312/150‑1 (MQ=21)
                   ggggAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTg                                                                                                                                     >  1:4325176/1‑150 (MQ=255)
                    gggAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGa                                                                                                                                   >  2:8551776/1‑150 (MQ=255)
                    gggAGGATAGAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGa                                                                                                                                   <  1:2959902/150‑1 (MQ=255)
                           tataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCGACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGtt                                                                                                                            >  2:3032260/1‑150 (MQ=255)
                            ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                           <  1:6371372/150‑1 (MQ=255)
                             taGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTcc                                                                                                                          <  1:6728014/150‑1 (MQ=255)
                               ggggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTa                                                                                                                        >  2:206128/1‑150 (MQ=255)
                                gggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAAAGCGCCGCTTGACCTACCAACGTATTCGCGGGTGAACCCGTTCCTac                                                                                                                       >  1:5075038/1‑150 (MQ=255)
                                 ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTAGTTGAACGGGCAACAGTCGTCATAAATCGAGGCAGGTGTAGGAGCGGGTTGACCCGCCAATGAGACCGTGGCCCGACCAAAGTATTAGCGTGTGAAACCCtacataca                                                                                                                      >  1:2323690/1‑150 (MQ=255)
                                     gTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGa                                                                                                                  >  2:5666299/1‑150 (MQ=255)
                                          aacaaaaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtc                                                                                                             <  1:3506790/147‑1 (MQ=255)
                                            gAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgt                                                                                                           <  2:181860/150‑1 (MQ=255)
                                             aaaaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgtc                                                                                                          >  1:7167310/1‑150 (MQ=255)
                                               aaaaGACAAATTGTCGCAGAGTATTGTTGTTGAACGGGAAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCTTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGc                                                                                                        <  2:713797/150‑1 (MQ=255)
                                               aaaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGc                                                                                                        <  1:8812246/150‑1 (MQ=255)
                                                aaaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCt                                                                                                       >  2:1919117/1‑150 (MQ=255)
                                                    aCAAATTGGCGCAGAGGATTGTTGTTGAACGGGCAACTGTCGTCATTAATCGAGGCAGGTTTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCACTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGTTCTTCGTCGCTtgtg                                                                                                   <  1:582986/150‑1 (MQ=39)
                                                      aaaTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGt                                                                                                 >  2:614253/1‑150 (MQ=255)
                                                       acttGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGtt                                                                                                <  1:6667702/148‑1 (MQ=255)
                                                              gCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAgg                                                                                         >  1:1985392/1‑150 (MQ=255)
                                                               cAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTGACCCGCGAATGGGCCGGTTGGCCTACCATCGTATTAGCGGGTGAAACCGTTGCTACAACGATCTTCGTGGCTTGTGGTTGAGAAggt                                                                                        >  1:8321222/1‑150 (MQ=255)
                                                                agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                       <  2:3340881/150‑1 (MQ=255)
                                                                agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                       <  1:1549775/150‑1 (MQ=255)
                                                                agagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtg                                                                                       <  2:5116759/150‑1 (MQ=255)
                                                                 gagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtgg                                                                                      <  1:5219592/150‑1 (MQ=255)
                                                                  agcaTTGTTCTTGAACGGGCAACAGTCGTCATTAAACGAGGCAGGTGTAGGAGCGGTTTCACCCGCGAAAGAGCCCGTTGCCCTACCATCGAATACGCCGGTGAACCAGATCATACACCGAAATACGTCGCTTGTGGATCAGAAggagta                                                                                     >  2:3401008/4‑146 (MQ=255)
                                                                    gATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTAc                                                                                   >  2:4769686/1‑150 (MQ=255)
                                                                        gttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCaa                                                                               <  2:4400874/150‑1 (MQ=255)
                                                                         ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAc                                                                              >  2:78046/1‑150 (MQ=255)
                                                                          tgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAcc                                                                             <  1:6555661/150‑1 (MQ=255)
                                                                           gttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAccc                                                                            <  2:6921338/150‑1 (MQ=255)
                                                                              gAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGAc                                                                         >  2:6729065/1‑150 (MQ=255)
                                                                               aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCCCGAAAGCGCCGGGTGCCCTACCAAAGTATTCGCGGGTGAAACCGTTCCTAAACCGATCATCGTCGCATGTTGTTGAGAATGTGGTAAGCAccccccct                                                                        >  1:6896007/1‑146 (MQ=255)
                                                                               aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTAACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACt                                                                        >  1:2155106/1‑150 (MQ=255)
                                                                               aaCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGGTAACCCCCGAAAGCGCCCGTTGCCCGACCAACGGAATCGCGGGTGTAACCCGTCCTACACCCCAAGTCGTCGCGGGTGGTTGAGAAGGTGGTTAGCAAaccccca                                                                        >  1:1502386/1‑146 (MQ=255)
                                                                                aCGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACtt                                                                       <  1:4880920/150‑1 (MQ=255)
                                                                                  gggCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACGTcc                                                                     >  2:8253888/1‑150 (MQ=255)
                                                                                    gCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCAc                                                                   >  2:2723500/1‑150 (MQ=255)
                                                                                        cAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTa                                                               <  1:2734012/150‑1 (MQ=255)
                                                                                             gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGccc                                                          <  1:2362329/150‑1 (MQ=255)
                                                                                              tcATTACTCGCCGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                                         <  1:4593472/150‑1 (MQ=255)
                                                                                              tcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                                         <  2:5555079/150‑1 (MQ=255)
                                                                                              tcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                                         >  2:8706153/1‑150 (MQ=255)
                                                                                               cATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCg                                                        <  2:7946792/150‑1 (MQ=255)
                                                                                                aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGCATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                       >  2:4085387/1‑150 (MQ=255)
                                                                                                   aaTCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGGTCAGAAGGTGGTACGCAACCCGACTTCCACGCGACGCCCCGgcgc                                                    >  2:3857413/1‑150 (MQ=255)
                                                                                                    aTCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGcc                                                   >  2:4257446/1‑150 (MQ=255)
                                                                                                         gatcggTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGc                                              >  2:2472288/5‑150 (MQ=17)
                                                                                                            aGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGat                                           >  2:829781/1‑150 (MQ=255)
                                                                                                               tgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATc                                        >  2:3994586/1‑150 (MQ=255)
                                                                                                                 tAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATcgc                                      >  2:6816893/1‑150 (MQ=255)
                                                                                                                   ggAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCgcgc                                    >  1:1080013/1‑150 (MQ=18)
                                                                                                                    gaacGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                   >  2:5118065/4‑150 (MQ=21)
                                                                                                                    gAGCGGGTTCACCCGCGCATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                   <  1:32580/150‑1 (MQ=255)
                                                                                                                    gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                   >  2:2652730/1‑150 (MQ=255)
                                                                                                                    gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                   <  2:5889223/150‑1 (MQ=255)
                                                                                                                    gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                   <  1:3704755/150‑1 (MQ=255)
                                                                                                                        gggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGAt                               <  1:1298272/150‑1 (MQ=255)
                                                                                                                         ggTTCACCCGCGACTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                              <  1:1410574/150‑1 (MQ=255)
                                                                                                                         ggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTTGTACGCCACCCGACTTCCAAGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATc                              <  2:2323690/150‑1 (MQ=255)
                                                                                                                           ttCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGa                            >  1:8465548/1‑150 (MQ=255)
                                                                                                                             cACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAAc                          >  1:4492685/1‑150 (MQ=255)
                                                                                                                              aCCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACt                         >  2:3742616/1‑150 (MQ=255)
                                                                                                                              aCCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACt                         <  2:6648313/150‑1 (MQ=255)
                                                                                                                               cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg                        >  1:4817863/1‑150 (MQ=255)
                                                                                                                                ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                       >  1:3084507/1‑150 (MQ=255)
                                                                                                                                ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                       >  1:3865504/1‑150 (MQ=255)
                                                                                                                                ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                       <  1:4345156/150‑1 (MQ=255)
                                                                                                                                ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                       >  1:4955204/1‑150 (MQ=255)
                                                                                                                                 cgcgAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGc                      >  1:7879019/1‑150 (MQ=255)
                                                                                                                                    gAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATa                   <  2:7999069/150‑1 (MQ=255)
                                                                                                                                    gAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATa                   <  2:6233482/150‑1 (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)
                                                                                                                                            cGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACgg           >  1:6480034/1‑150 (MQ=255)
                                                                                                                                                  cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggt     <  1:8551776/150‑1 (MQ=255)
                                                                                                                                                  cccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggt     <  2:1320348/150‑1 (MQ=255)
                                                                                                                                                   ccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtc    >  1:5186738/1‑150 (MQ=255)
                                                                                                                                                   ccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtc    >  2:6359123/1‑150 (MQ=255)
                                                                                                                                                   ccTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtc    >  2:3080699/1‑150 (MQ=255)
                                                                                                                                                    cTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcg   >  1:5548269/1‑150 (MQ=255)
                                                                                                                                                    cTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGGTCAGAAGGTGGTACGCAACCCGACTTCCCCGCTACGCCCCGGCGCCGGTGCGATATCGCGCCGGATCGAACTGGCATAGCGCACGGTCTggtcg   >  2:3304969/1‑150 (MQ=255)
                                                                                                                                                     tACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcgg  >  1:2236722/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGG  >  NC_002947/1653804‑1654102

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: