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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,649,0090TG 5.7% 421.0 / 11.2 142intergenic (‑246/+50)oprB‑II/PP_1446carbohydrate‑selective porin/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (55/71);  new base G (5/4);  total (67/75)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 145 total aligned reads displayed.

TTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAGG  >  NC_002947/1648860‑1649158
                                                                                                                                                     |                                                                                                                                                     
tttCAGCTTGCGGCGGGGAGGATATAGGGGAAGGTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGTTTCACCCGCGAATGCGCCGGTTGCCCTACCAt                                                                                                                                                       <  1:5161712/150‑1 (MQ=255)
 ttCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATc                                                                                                                                                      <  2:4994427/150‑1 (MQ=255)
      cTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                                                 <  2:165890/150‑1 (MQ=255)
      cTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTAtt                                                                                                                                                 <  2:165772/150‑1 (MQ=255)
           ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                            <  2:3705268/150‑1 (MQ=255)
           ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                            <  1:153633/150‑1 (MQ=255)
           ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCgg                                                                                                                                            <  2:4359134/150‑1 (MQ=255)
           ggcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAACGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCgg                                                                                                                                            <  2:2995660/150‑1 (MQ=25)
            gcggGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCggg                                                                                                                                           <  2:1418335/150‑1 (MQ=25)
                 gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                                      <  2:2743446/150‑1 (MQ=255)
                 gaggaTATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAAc                                                                                                                                      >  1:5214424/1‑150 (MQ=255)
                     atatTGGTGAAGTTTTAATAAAAAATAAAAATTGTCGCAGAGGATTGTTGTTTAACGGGCAACTTTCGTCATTAATCGAGGCAGTTGTAGGTGCGGTTTCACCCGCGAATACGTTGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCgc                                                                                                                                  <  1:3825113/150‑2 (MQ=12)
                     atatAGGGGAAGTTATCAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGt                                                                                                                                  <  1:2966700/150‑1 (MQ=255)
                       ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                                <  1:5355374/150‑1 (MQ=255)
                       ataGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTc                                                                                                                                <  1:2994195/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                                                                                                                                >  2:4472818/1‑145 (MQ=255)
                            ggAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTaca                                                                                                                           <  2:1281576/150‑1 (MQ=255)
                             gAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTacac                                                                                                                          <  1:1964902/150‑1 (MQ=255)
                              aaGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACAcc                                                                                                                         >  1:4371978/1‑150 (MQ=255)
                                 ttATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCTAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGAt                                                                                                                      <  1:1015155/150‑1 (MQ=255)
                                 ttATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGAt                                                                                                                      <  1:2689450/150‑1 (MQ=255)
                                    taAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCtt                                                                                                                   >  2:3665117/1‑150 (MQ=255)
                                      aGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCGGTTGCCCTACCAACGTATTAGCGGGTGAACCCCGTCCAACACCCATCttag                                                                                                                 >  1:1134650/1‑148 (MQ=255)
                                       gAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTCCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTtcgt                                                                                                                >  2:2771992/1‑150 (MQ=255)
                                            aaGACCAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCCGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                                           >  2:3269482/1‑150 (MQ=255)
                                            aaGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCtt                                                                                                           <  1:4244974/150‑1 (MQ=255)
                                             aGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTtg                                                                                                          <  2:4285843/150‑1 (MQ=255)
                                               aCAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTtgtg                                                                                                        <  2:1803253/150‑1 (MQ=255)
                                               aCAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCCGTTGCCCTACCAACGTATTAGCGGGTGAACCCCTTCCTACACCGATCTTCGTCGCTtgtg                                                                                                        >  1:4700125/1‑150 (MQ=255)
                                                 aaaTTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCAACCGCGAATGCGCCCGTTGCCCTACCAACGTATTCGCGGGTGAACCCGGTCCTACACCGATAGTCGTCGCGTGTGGt                                                                                                      >  1:3419599/1‑150 (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)
                                                            tagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtgg                                                                                           <  1:5261427/149‑1 (MQ=255)
                                                            gagGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAggtgg                                                                                           <  1:3475867/150‑1 (MQ=255)
                                                              ggATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTa                                                                                         <  1:1627540/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)
                                                                    ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATc                                                                                                                     <  1:2240969/116‑1 (MQ=255)
                                                                    ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATc                                                                                                                     >  2:2240969/1‑116 (MQ=255)
                                                                    ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTc                                                                                                    <  1:4798347/133‑1 (MQ=255)
                                                                    ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTc                                                                                                    >  2:4798347/1‑133 (MQ=255)
                                                                    ttgttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAc                                                                                   >  2:3891934/1‑150 (MQ=255)
                                                                      gttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAccc                                                                                 >  2:3177593/1‑150 (MQ=255)
                                                                      gttgAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAAccc                                                                                 <  2:27196/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)
                                                                                      tcgtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGc                                                                 >  2:1263597/1‑150 (MQ=255)
                                                                                        gtcATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCCCGAATGCGCCCGTTTGCCTACCATCGTATTCGCGGGTGAACCCGGTCCAACACCGGTCTTCGTCGCGTGTGGTTAAGAAGGTGGTAAGCAACCCGCCATCCACGCTACGccc                                                               >  1:2874394/1‑150 (MQ=255)
                                                                                           aTTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCgg                                                            >  2:153681/1‑150 (MQ=255)
                                                                                            ttAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGCATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCGTGTGGTTCAGAAGGTGGTACGCCACCCGACTTCCACGCTACGCCCCGgc                                                           >  1:2668139/1‑150 (MQ=255)
                                                                                                  gAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCCTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGt                                                     <  2:2328852/150‑1 (MQ=255)
                                                                                                           gtAGGAGCGGGTTCACCCGCGAAATCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATcg                                            >  2:2436726/1‑150 (MQ=255)
                                                                                                            tAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATagc                                           >  1:2251122/1‑150 (MQ=255)
                                                                                                              ggAACGGGTTAACCCGCGAATACGATGGTAGGGCAACCGGCGCATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCgcgc                                         >  1:1054804/1‑150 (MQ=21)
                                                                                                               gAGCGGGTTCACCCGCGAATTCCCCGGATGCCCTACCATCGTATTCGCCGGTGAACCCGTTCCTACACCGATCTTAGTCGCTTGTGGTTCAGAAGGTGGTAAGCAACCCGACAACCACGCTACGCCCCGGCGCCCGTGCGATATCGCGCa                                        >  2:1681457/1‑150 (MQ=255)
                                                                                                               gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGc                                                   <  1:1786564/139‑1 (MQ=255)
                                                                                                               gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                        <  2:3610747/150‑1 (MQ=255)
                                                                                                               gAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCa                                        <  2:2576431/150‑1 (MQ=255)
                                                                                                               gAGCGGGTTCACCCGCGAATGCGCCCGGTGCCCGACCCTCGTAATCGCCGGTGAACCCCTTCCAACACCGATCGTCGTCGCGTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCGACGCCCCGGCGCCCGTGCGATATCgcgcc                                        >  1:4082255/1‑149 (MQ=255)
                                                                                                                 gCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAgg                                      >  1:3307938/1‑150 (MQ=255)
                                                                                                                    ggTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGcccc                                                              >  2:4931224/1‑123 (MQ=255)
                                                                                                                     gTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCg                                  >  2:1599800/1‑150 (MQ=255)
                                                                                                                      ttCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGa                                 <  1:709409/150‑1 (MQ=255)
                                                                                                                      ttCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGa                                 >  1:2394322/1‑150 (MQ=255)
                                                                                                                          cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg                             >  1:245150/1‑150 (MQ=255)
                                                                                                                          cccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTg                             <  2:362541/150‑1 (MQ=255)
                                                                                                                           ccGCGACTGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGCCTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            <  1:693918/150‑1 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            <  1:712832/150‑1 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:1787322/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:1858324/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  1:1194303/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTgg                            >  2:2143458/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)
                                                                                                                                   gcgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAgcgc                    >  1:4975443/1‑150 (MQ=255)
                                                                                                                                   gcgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAgcgc                    >  2:2085285/1‑150 (MQ=255)
                                                                                                                                    cgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCa                   >  1:3323291/1‑150 (MQ=255)
                                                                                                                                    cgcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCa                   >  2:2928876/1‑150 (MQ=255)
                                                                                                                                     gcCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCAc                  >  1:1762911/1‑150 (MQ=255)
                                                                                                                                         gTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTc              <  2:2639950/150‑1 (MQ=255)
                                                                                                                                          ttGCCCTACCATAGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCt             <  1:4321154/150‑1 (MQ=255)
                                                                                                                                           tGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTg            <  2:3419599/150‑1 (MQ=255)
                                                                                                                                               cTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcg        <  2:2355463/150‑1 (MQ=255)
                                                                                                                                                tACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCg                                                       <  2:3317785/102‑1 (MQ=255)
                                                                                                                                                tACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcgg       >  2:2642472/1‑150 (MQ=255)
                                                                                                                                                 aCCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGTGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTggtcggt      >  1:1256639/1‑150 (MQ=255)
                                                                                                                                                    aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAg   <  2:349786/150‑1 (MQ=255)
                                                                                                                                                    aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAg   <  1:3177593/150‑1 (MQ=255)
                                                                                                                                                    aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAg   <  1:2801706/150‑1 (MQ=255)
                                                                                                                                                    aTCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAg   <  2:1134650/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCAGCTTGCGGCGGGGAGGATATAGGGGAAGTTATAAGAAAAAAGACAAATTGTCGCAGAGGATTGTTGTTGAACGGGCAACAGTCGTCATTAATCGAGGCAGGTGTAGGAGCGGGTTCACCCGCGAATGCGCCGGTTGCCCTACCATCGTATTCGCGGGTGAACCCGTTCCTACACCGATCTTCGTCGCTTGTGGTTCAGAAGGTGGTACGCAACCCGACTTCCACGCTACGCCCCGGCGCCGGTGCGATATCGCGCAGGATCGAACTGGCATAGCGCACGGTCTGGTCGGTCAGG  >  NC_002947/1648860‑1649158

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: