Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,073,173 T→A 5.3% intergenic (+160/+67) PP_1846 → / ← PP_1847 group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,073,1730TA 5.3% 404.1 / 8.0 130intergenic (+160/+67)PP_1846/PP_1847group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (59/61);  new base A (2/4);  total (61/69)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 132 total aligned reads displayed.

CGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACA  >  NC_002947/2073025‑2073321
                                                                                                                                                    |                                                                                                                                                    
cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTc                                                                                                                                                     <  1:1377236/150‑1 (MQ=255)
cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTc                                                                                                                                                     <  1:5113672/150‑1 (MQ=255)
cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTc                                                                                                                                                     >  1:5127112/1‑150 (MQ=255)
cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTc                                                                                                                                                     <  2:1513014/150‑1 (MQ=255)
 ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGTATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCg                                                                                                                                                    <  1:4896737/150‑1 (MQ=255)
   aaCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGcc                                                                                                                                                  >  1:4481310/1‑150 (MQ=255)
    aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                                 <  1:4474635/150‑1 (MQ=255)
    aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                                 <  1:4407280/150‑1 (MQ=255)
    aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                                 <  2:1275242/150‑1 (MQ=255)
    aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                                 <  2:1001044/150‑1 (MQ=255)
      cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCgt                                                                                                                                               >  1:4245169/1‑150 (MQ=255)
      cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCgt                                                                                                                                               <  1:5227438/150‑1 (MQ=255)
         aCGCACGGTGGGGTGCGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGtt                                                                                                                                            <  1:4709544/150‑1 (MQ=255)
            cACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCGAAGCCCgtgccatc                                                                                                                                         <  2:4128599/150‑6 (MQ=255)
             aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTcgcg                                                                                                                                        <  2:2155292/150‑1 (MQ=255)
             aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTcgcg                                                                                                                                        >  2:1961098/1‑150 (MQ=255)
             aCGGTGGGGTGAGAGGTCGGCGGGTGTGACGCCGCCGCCGCCTCGATTGTCCGGGCCATGGGGGGTGCTGGGCAGCCCTTCGCGACACAAGCGGGCTGGCCCGCGAAGCCAACACCGGGGGGGATGGCACGGGCTTCGCCCGTGTTcgcg                                                                                                                                        <  1:562757/150‑1 (MQ=21)
              cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCACCACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                       <  2:4273420/150‑1 (MQ=255)
              cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                       <  2:3282442/150‑1 (MQ=255)
              cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                       <  2:1098024/150‑1 (MQ=255)
              cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                       <  2:2947465/150‑1 (MQ=255)
                  ggtgTGAGAGGACGGCGGGGGGGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCaa                                                                                                                                   <  1:3445889/150‑1 (MQ=255)
                    tgtgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGc                                                                                                                                 <  1:3217046/150‑1 (MQ=255)
                        agagGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCCAGCCCGc                                                                                                                             >  1:5243359/1‑150 (MQ=255)
                         gagGACGGCGGGTGTGAATCCGCCGCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCt                                                                                                                            >  1:160086/1‑150 (MQ=255)
                          agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTc                                                                                                                           >  1:2754167/1‑150 (MQ=255)
                          agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGAGTTCGCGGGCAAGCCCGCTc                                                                                                                           >  2:1165584/1‑150 (MQ=255)
                            gACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTccc                                                                                                                         <  1:5446410/150‑1 (MQ=255)
                               ggcggGTGTTAATCCGCCTCCTACTCTTTTTTCCGGGCCTTTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCAGGGGCTTCGCCCGTTTTCGCGGGCAAGCCCGCTCCcaca                                                                                                                      <  1:4501360/150‑1 (MQ=25)
                               ggcggGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCcaca                                                                                                                      >  1:2453096/1‑150 (MQ=255)
                                    gtgtGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTcg                                                                                                                 >  1:4132313/1‑150 (MQ=255)
                                    gtgtGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTcg                                                                                                                 >  1:4132429/1‑150 (MQ=255)
                                          aTCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAg                                                                                                           >  2:2364318/1‑150 (MQ=255)
                                          aTCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAg                                                                                                           >  2:5230821/1‑150 (MQ=255)
                                              gCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGcc                                                                                                       >  2:356374/1‑150 (MQ=255)
                                               cctcctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCg                                                                                                      >  1:2741695/1‑150 (MQ=255)
                                                   ctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCaga                                                                                                  >  1:4159828/1‑150 (MQ=255)
                                                   ctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCaga                                                                                                  >  2:4171939/1‑150 (MQ=255)
                                                   ctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTGCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCaga                                                                                                  <  2:5127112/150‑1 (MQ=255)
                                                     aCTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGc                                                                                                <  2:2662275/150‑1 (MQ=255)
                                                     aCTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGc                                                                                                >  1:5334191/1‑150 (MQ=255)
                                                      cTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                                                                                               >  1:4323656/1‑150 (MQ=255)
                                                       tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCgg                                                                                              >  2:3458479/1‑150 (MQ=255)
                                                       tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCgg                                                                                              >  1:3042671/1‑150 (MQ=255)
                                                                    ccATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCCATCTa                                                                                 >  1:5249215/1‑150 (MQ=255)
                                                                    ccATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTa                                                                                 >  1:5249042/1‑150 (MQ=255)
                                                                     atttGGGGGTGCTGTGCAGCCAATCGCTACACATGGGGGGTTGCCCGCGAAGCCAACACCGTGGTGTTTGGCAGGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAg                                                                                <  1:3550505/148‑1 (MQ=25)
                                                                       ttGGGGCTGCTGTGCCGCCCATCGCGACACAAGCGGGCTTGCCGGCGAAGCCAACACCGTGGTGGATGGAACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGaa                                                                              <  2:3659735/150‑1 (MQ=255)
                                                                       ttGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGaa                                                                              <  2:2234936/150‑1 (MQ=255)
                                                                       tgggggCTGGTGTGCAGCCCTTCGCGACACAAGGGGGGTGGCCCGCGACGCCAACACCGGGGTGGTTGGCACGGGCTTCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGaa                                                                              <  2:2878240/148‑1 (MQ=37)
                                                                            gctgctGTGCAGCCCATCGCGACACACGGGGGCTTGCCCGCGAAGCCAACACGGTGGTGGTTGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGa                                                                         <  2:3473120/150‑1 (MQ=255)
                                                                            gctgctGTGCAGCCCATCGCGAAACCCGCGGGCTTGCCCGCGACTCCACCACCTTGGTGGATGGAACGGGCTTCGCCGTTGTGCGCGGTCACGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGa                                                                         <  2:738281/150‑1 (MQ=255)
                                                                            gctgctGTGCAGCCCATCGCCACACACGGGGGCTGCCCCGCTAAGCCAACACCTGGGTGGTGGGCACGGGCTCCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAACGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGa                                                                         <  2:1542792/150‑1 (MQ=25)
                                                                            gctgcgGTGCCGCCCACCGCGACACACGGGGGCTTGCCCGCGAAGCCAACCCCGGGGTGTATGGCAGGGGCTCCGCCCGTGTCCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGa                                                                         <  2:1543086/150‑1 (MQ=25)
                                                                                gtttgCCGCCCACCGCGAAACAAGGGGGCTGGCCCGCCACCCAAACCCGGTGTGGGTTGGCCGGGGCTGCGCCGGTGTTGGGGGGCACGCCCGCTCCCACCGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGa                                                                     <  1:3254611/147‑1 (MQ=255)
                                                                                 tgtgCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaa                                                                    >  2:3811956/1‑150 (MQ=255)
                                                                                 tgtgCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaa                                                                    >  2:3184358/1‑150 (MQ=255)
                                                                                  gtgCAGCCCATTGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCca                                                                                                                        >  1:5371392/1‑97 (MQ=255)
                                                                                        cccATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaa                                                             <  2:4085443/150‑1 (MQ=255)
                                                                                         ccATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaag                                                            >  1:5356768/1‑150 (MQ=255)
                                                                                           aTCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaagcc                                                          >  1:3502955/1‑150 (MQ=255)
                                                                                                   acaAGCGGGCTGGCCCGCTAAGCCAACACCTTGGTGCATGGCACGGGATTCGCCCTTTTTCGCGGCCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTa                                                  <  2:666695/150‑1 (MQ=255)
                                                                                                   acaAGCGGCCTTGCCCGCCACCCAAACACCGTGGTGTATGCAACGGGCTTCCCCCGTGTACGCGGGCAACCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTa                                                  <  2:666776/150‑1 (MQ=255)
                                                                                                     aaGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                       >  1:1136730/1‑63 (MQ=37)
                                                                                                     aaGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                                       <  2:1136730/63‑1 (MQ=37)
                                                                                                     aaGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCACTCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTAcg                                                >  1:4272978/1‑150 (MQ=255)
                                                                                                       gCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTAcgcg                                              <  2:2453096/150‑1 (MQ=255)
                                                                                                             ttGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGccc                                        >  1:2671722/1‑150 (MQ=255)
                                                                                                                 ccGCGACCCCACCACGGGGGGGGATGGCCCGGGCTTCGCCCGTGTTGGCGGGCACCCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTg                                    <  1:1659558/150‑1 (MQ=255)
                                                                                                                 ccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTtgtgtcgcgatgggctgcacagcagccccaatggccctgaaaatcgagtaggaggcggattcacacccgccgtcctctcacaccaccg                                    <  2:2076656‑M2/150‑89 (MQ=255)
                                                                                                                 ccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTtgtgtcgcgatgggctgcacagcagccccaatggccctgaaaatcgagtaggaggcggattcacacccgccgtcctctcacaccaccg                                    <  1:3624669‑M2/150‑89 (MQ=255)
                                                                                                                 ccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTg                                    >  2:4564437/1‑150 (MQ=255)
                                                                                                                 ccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTg                                    >  1:2965689/1‑150 (MQ=255)
                                                                                                                 ccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGGTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCGAGAAGCTGTTCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTg                                    >  1:372637/1‑150 (MQ=255)
                                                                                                                   gcgAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGa                                  >  1:2870999/1‑150 (MQ=255)
                                                                                                                    cgAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGaa                                 <  2:3042671/150‑1 (MQ=255)
                                                                                                                       aGCCAACCCCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAg                              <  2:5356768/150‑1 (MQ=255)
                                                                                                                       aGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTtgtgtcgcgatgggctgcacagcagccccaatggccctgaaaatcgagtaggaggcggattcacacccgccgtcctctcacaccaccgtgcgta                              <  1:3140871‑M2/150‑95 (MQ=255)
                                                                                                                       aGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAg                              >  2:3435787/1‑150 (MQ=255)
                                                                                                                          cAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                           <  2:1277884/150‑1 (MQ=255)
                                                                                                                          cAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                           <  2:5439687/150‑1 (MQ=255)
                                                                                                                          cAACACCGTGGTGGATAGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                           <  2:3595652/150‑1 (MQ=255)
                                                                                                                            acacCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAAGCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGACCAGGTCGt                         >  2:4410223/1‑150 (MQ=255)
                                                                                                                             cacCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGtt                        >  2:2643813/1‑150 (MQ=255)
                                                                                                                                 gtggtgGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATcc                    >  1:2464318/1‑150 (MQ=255)
                                                                                                                                  tggtggTTGGCAGGGGCTTCGCCGGTGTTCGGGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATccc                   <  1:3521876/150‑1 (MQ=255)
                                                                                                                                     tggATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCa                <  2:2754167/150‑1 (MQ=255)
                                                                                                                                         tGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTa            >  2:4755380/1‑150 (MQ=255)
                                                                                                                                           gCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAAt          >  2:1983782/1‑150 (MQ=255)
                                                                                                                                           gCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAAt          >  2:4384793/1‑150 (MQ=255)
                                                                                                                                           gCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAAt          >  2:4513531/1‑150 (MQ=255)
                                                                                                                                           acACGGGCGAAGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                           >  2:760839/2‑133 (MQ=255)
                                                                                                                                           acACGGGCGAAGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                           >  2:760926/2‑133 (MQ=255)
                                                                                                                                           acACGGGCGAAGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                           <  1:760839/132‑1 (MQ=255)
                                                                                                                                            cACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTg                                           <  2:354367/116‑1 (MQ=255)
                                                                                                                                            cACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTCATc         >  2:1157393/1‑150 (MQ=255)
                                                                                                                                                 gCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGACCAGGTAGTTATCCCCCAGGTAATCCCCca    <  2:3772262/150‑1 (MQ=255)
                                                                                                                                                 gCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCca    <  2:3772257/150‑1 (MQ=255)
                                                                                                                                                   ttCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCcaca  >  1:4760364/1‑150 (MQ=255)
                                                                                                                                                   ttCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCcaca  >  1:2555176/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCACA  >  NC_002947/2073025‑2073321

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: