Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,648,991 G→A 5.2% intergenic (‑228/+68) 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,648,9910GA 5.2% 406.8 / 7.6 135intergenic (‑228/+68)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base G (57/69);  new base A (3/4);  total (62/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 140 total aligned reads displayed.

GATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGC  >  NC_002947/1648843‑1649140
                                                                                                                                                    |                                                                                                                                                     
gATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATgc                                                                                                                                                      >  2:5355374/1‑150 (MQ=255)
 aTTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATgcg                                                                                                                                                     >  1:2355463/1‑150 (MQ=255)
  ttGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAAGTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATgcgc                                                                                                                                                    <  2:5087983/150‑1 (MQ=255)
     ttCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCgg                                                                                                                                                 >  1:352458/1‑150 (MQ=255)
       cGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGcc                                                                                                                                                   >  2:1670366/1‑146 (MQ=255)
       cGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGcc                                                                                                                                                   <  1:1670366/146‑1 (MQ=255)
           gAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGcc                                                                                                                                            >  1:1545975/1‑149 (MQ=255)
            aaGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCt                                                                                                                                          >  1:2576431/1‑150 (MQ=255)
            aaGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCt                                                                                                                                          >  1:349786/1‑150 (MQ=255)
                gTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCa                                                                                                                                      <  2:5227406/150‑1 (MQ=255)
                 tttCAGCTTGCGGCGGGGAGGATATAGGGGAAGGTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGTTTCACCCGCGAATGCGCCGGTTGCCCTACCAt                                                                                                                                     <  1:5161712/150‑1 (MQ=255)
                    cAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGt                                                                                                                                  >  2:2943269/1‑150 (MQ=255)
                       cTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                               <  2:165890/150‑1 (MQ=255)
                       cTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                               <  1:4755625/150‑1 (MQ=255)
                            ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          <  1:1837007/150‑1 (MQ=255)
                            ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          <  1:153633/150‑1 (MQ=255)
                            ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          <  2:258691/150‑1 (MQ=255)
                            ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          <  2:3705268/150‑1 (MQ=255)
                            ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          <  1:5088741/150‑1 (MQ=255)
                            ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                          <  2:2995660/150‑1 (MQ=25)
                                  gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                    >  1:5214424/1‑150 (MQ=255)
                                  gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                    <  2:2743672/150‑1 (MQ=255)
                                  gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                    <  2:2743446/150‑1 (MQ=255)
                                  gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAAc                                                                                                                    >  1:2486332/1‑150 (MQ=255)
                                     gaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCg                                                                                                                 >  2:481076/1‑150 (MQ=255)
                                      atatTGGTGAAGTTTTAATAAAAAATAAAAATTGTCGCAGAGGATTGTTGTTTAACGGGCAACTTTCGTCATTAATCGAGGCAGTTGTAGGTGCGGTTTCACCCGCGAATACGTTGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCgc                                                                                                                <  1:3825113/150‑2 (MQ=12)
                                        ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                              <  1:2994195/150‑1 (MQ=255)
                                        ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                              <  1:5355374/150‑1 (MQ=255)
                                        ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                              <  2:2332874/150‑1 (MQ=255)
                                           ggggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTa                                                                                                           <  2:4368456/150‑1 (MQ=255)
                                             ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                              >  2:4473127/1‑145 (MQ=255)
                                             ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                              <  1:4473127/145‑1 (MQ=255)
                                             ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                              <  1:4472818/145‑1 (MQ=255)
                                             ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTaca                                                                                                         <  2:1281576/150‑1 (MQ=255)
                                               aaGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTTAACCCGTTCCTACAcc                                                                                                       <  1:223571/150‑1 (MQ=255)
                                                  ttATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCTAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGAt                                                                                                    <  1:1015155/150‑1 (MQ=255)
                                                  ttATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGAt                                                                                                    <  1:2689450/150‑1 (MQ=255)
                                                       aGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCGGTTGCCCTACCAACGTATTAGCGGGTGAACCCCGTCCAACACCCATCttag                                                                                               >  1:1134650/1‑148 (MQ=255)
                                                             aaGACCAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCCGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                         >  2:3269482/1‑150 (MQ=255)
                                                             aaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                         <  1:4245170/150‑1 (MQ=255)
                                                             aaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                         <  1:4244974/150‑1 (MQ=255)
                                                              aGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTtg                                                                                        <  2:4285843/150‑1 (MQ=255)
                                                                aCAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTtgtg                                                                                      <  2:1803253/150‑1 (MQ=255)
                                                                   acttGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGtt                                                                                   <  2:3623111/148‑1 (MQ=255)
                                                                     ttGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCa                                                                                 <  2:1427403/150‑1 (MQ=255)
                                                                     ttGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCAACGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCa                                                                                 >  2:1832920/1‑150 (MQ=255)
                                                                          gCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAgg                                                                            >  2:848351/1‑150 (MQ=255)
                                                                           cAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggt                                                                           <  2:1012056/150‑1 (MQ=255)
                                                                             gagGATTGTTGTTGAACGGGCCACAGTCGGCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGTATGCGCCGGTTGCCCTACCATAGTATTCCCGGGTGAACCAGTTCCTACACCGATCTTCGTCGCTTGTGGGTCAGAAggtgg                                                                         >  1:2715775/1‑150 (MQ=255)
                                                                             gagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtgg                                                                         <  1:3475867/150‑1 (MQ=255)
                                                                                gATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCGTGTAGATCGGAAgtgcaccc                                                                      >  1:4948190/1‑143 (MQ=38)
                                                                                   tgttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCa                                                                   >  1:1290398/1‑150 (MQ=255)
                                                                                     ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          >  1:3679230/1‑93 (MQ=255)
                                                                                     ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                          <  2:3679230/93‑1 (MQ=255)
                                                                                     ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATc                                                                                                   <  1:2240969/116‑1 (MQ=255)
                                                                                     ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAc                                                                 >  2:3891934/1‑150 (MQ=255)
                                                                                     ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAc                                                                 <  2:398580/150‑1 (MQ=255)
                                                                                       gttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAccc                                                               >  2:3177593/1‑150 (MQ=255)
                                                                                        ttgAACGGTCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGTGAACCCGGTCCAACACCGAACATCGTCGCATGTGGGTCAGAAGGTGGTACGCAAACCGa                                                             >  2:1625526/1‑150 (MQ=255)
                                                                                         tgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGa                                                             <  1:2257276/150‑1 (MQ=255)
                                                                                              gggCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTcc                                                        <  2:1290398/150‑1 (MQ=255)
                                                                                               ggCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCa                                                       <  1:3217751/150‑1 (MQ=255)
                                                                                               ggCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCCGTTTCCCGACCAACGTATTAGCGGGTGAAACCCTTCCAACACCGATATTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCCACTTACa                                                       >  1:4766165/1‑150 (MQ=255)
                                                                                                gCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGCATGCGCCGGTTGCCCTACCAGCGGATTCGCGGGTGAACCCGGTCCGACACCGCTCTTCGTCGCGTGGGGTTCAGAAGGGGGTACGCAACCCCACATCCAc                                                      >  2:693918/1‑150 (MQ=255)
                                                                                                       tcgtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                         <  1:4794026/108‑1 (MQ=255)
                                                                                                       tcgtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                         >  2:4794026/1‑108 (MQ=255)
                                                                                                         gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGccc                                             >  2:2179943/1‑150 (MQ=255)
                                                                                                         gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCCGTTTGCCTACCATCGTATTCGCGGGTGAACCCGGTCCAACACCGGTCTTCGTCGCGTGTGGTTAAGAAGGTGGTAAGCAACCCGCCATCCACGCTACGccc                                             >  1:2874394/1‑150 (MQ=255)
                                                                                                            aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                          >  2:153681/1‑150 (MQ=255)
                                                                                                                  gaagatcggTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGCTCCTACACCTATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCgg                                    >  2:3454348/8‑150 (MQ=14)
                                                                                                                   gAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCCTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGt                                   <  2:2328852/150‑1 (MQ=255)
                                                                                                                    aGGCAGGTGTAGGCGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAGTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTg                                  <  1:1688374/150‑1 (MQ=255)
                                                                                                                          gtgtAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCCTTCCTACACCGATCTTCGTCGCGTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCAACCCCCCGGCGCCGGTGCGatat                            >  1:3072855/1‑150 (MQ=255)
                                                                                                                               ggAACGGGTTAACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCgcgc                       >  1:1054804/1‑150 (MQ=21)
                                                                                                                                gAGCGGGTTCACCCGCGAATTCCCCGGATGCCCTACCATCGTATTCGCCGGTGAACCCGTTCCTACACCGATCTTAGTCGCTTGTGGTTCAGAAGGTGGTAAGCAACCCGACAACCACGCTACGCCCCGGCGCCCGTGCGATATCGCGCa                      >  2:1681457/1‑150 (MQ=255)
                                                                                                                                gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                      <  2:2576431/150‑1 (MQ=255)
                                                                                                                                gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                      <  2:1505273/150‑1 (MQ=255)
                                                                                                                                gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                      <  2:3610747/150‑1 (MQ=255)
                                                                                                                                gAGCGGGTTCACCCGCGAATGCGCCCGGTGCCCGACCCTCGTAATCGCCGGTGAACCCCTTCCAACACCGATCGTCGTCGCGTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCGACGCCCCGGCGCCCGTGCGATATCgcgcc                      >  1:4082255/1‑149 (MQ=255)
                                                                                                                                  gCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAgg                    >  1:3307938/1‑150 (MQ=255)
                                                                                                                                   cGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGa                   <  2:4371978/150‑1 (MQ=255)
                                                                                                                                     ggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                            >  2:4931224/1‑123 (MQ=255)
                                                                                                                                     ggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                            <  1:4931224/123‑1 (MQ=255)
                                                                                                                                      gTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCg                >  2:1599800/1‑150 (MQ=255)
                                                                                                                                       ttCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGa               <  1:709409/150‑1 (MQ=255)
                                                                                                                                       ttCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGa               >  1:2394322/1‑150 (MQ=255)
                                                                                                                                          aCCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCg                                     >  1:665920/1‑125 (MQ=21)
                                                                                                                                           cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg           <  2:362541/150‑1 (MQ=255)
                                                                                                                                           cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg           <  1:2958084/150‑1 (MQ=255)
                                                                                                                                            ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg          <  1:712832/150‑1 (MQ=255)
                                                                                                                                            ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg          >  1:1194303/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg          >  2:2143458/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg          >  1:1787322/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg          >  1:1858324/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATAGCGCAGGATCGAACTgg          >  1:2942354/1‑150 (MQ=255)
                                                                                                                                              gcgAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCa        <  1:3244714/150‑1 (MQ=255)
                                                                                                                                                gAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATa      <  1:1078209/150‑1 (MQ=255)
                                                                                                                                                  aTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGCTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAgc    <  2:2718001/150‑1 (MQ=255)
                                                                                                                                                  aTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGACGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACCCCCCGGCGCCGGTGCGATATCGCGCAGGAACGAACTGGCATAgc    >  2:4619669/1‑150 (MQ=255)
                                                                                                                                                    gcgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAgcgc  >  2:2085285/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATTGTTCGATGAAGCGTTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGC  >  NC_002947/1648843‑1649140

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: