Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,073,861 C→G 6.9% L214F (TTG→TTC pstB‑II ← phosphate ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,073,8610CG 6.9% 559.5 / 22.0 171L214F (TTG→TTCpstB‑IIphosphate ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (92/67);  new base G (6/6);  total (98/73)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

AACAGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTCGT  >  NC_002947/6073712‑6074009
                                                                                                                                                     |                                                                                                                                                    
aaCAGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTc                                                                                                                                                      >  2:5392252/1‑150 (MQ=255)
aaCAGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTc                                                                                                                                                      >  1:3774007/1‑150 (MQ=255)
  cAGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCaa                                                                                                                                                    >  2:1202046/1‑150 (MQ=255)
       tatcggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCg                                                                                                                                               <  1:3203433/150‑1 (MQ=255)
       tatcggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCg                                                                                                                                               <  1:3203198/150‑1 (MQ=255)
          cggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTGAAGGTCgaag                                                                                                                                            <  2:3250038/150‑3 (MQ=255)
          cggtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCCTCGACCTTGAAGGTCgaag                                                                                                                                            <  2:3105365/150‑3 (MQ=255)
            gtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATc                                                                                                                                          <  2:4295313/150‑1 (MQ=255)
            gtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATc                                                                                                                                          >  2:3158105/1‑150 (MQ=255)
            gtatcCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTACTTCGACCTTCAAGGTCGAGATc                                                                                                                                          <  1:1303548/150‑1 (MQ=255)
                  ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCaa                                                                                                                                                    >  1:1338846/1‑134 (MQ=255)
                  ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCaa                                                                                                                                                    <  2:1338955/134‑1 (MQ=255)
                  ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTcc                                                                                                                                    >  2:886641/1‑150 (MQ=255)
                  ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTcc                                                                                                                                    >  2:4337716/1‑150 (MQ=255)
                  ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTcc                                                                                                                                    <  1:1709573/150‑1 (MQ=255)
                  ccGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTcc                                                                                                                                    >  1:2665849/1‑150 (MQ=255)
                    gAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAg                                                                                                                                  <  1:2005474/150‑1 (MQ=255)
                     aaTTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAgc                                                                                                                                 <  1:1584930/150‑1 (MQ=255)
                        tCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGcc                                                                                                                              <  2:1360594/150‑1 (MQ=255)
                          gACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGa                                                                                                                            >  1:2404817/1‑150 (MQ=255)
                          gACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGa                                                                                                                            >  2:4265509/1‑150 (MQ=255)
                             cAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCa                                                                                                                         <  2:835320/150‑1 (MQ=255)
                             cAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCa                                                                                                                         >  2:5441191/1‑150 (MQ=255)
                               gTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACg                                                                                                                       <  1:554137/150‑1 (MQ=255)
                                     ccATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGt                                                                                                                 <  1:886641/150‑1 (MQ=255)
                                     ccATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGt                                                                                                                 >  1:4195929/1‑150 (MQ=255)
                                       aTGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCa                                                                                                               <  2:309697/150‑1 (MQ=255)
                                        tGTACATGTACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGTTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGTCCTTCCAGGTCGAGATCGGGTCCATCGCCGAGCACGGTTCGTCaa                                                                                                              <  1:1566739/150‑1 (MQ=255)
                                        tGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCaa                                                                                                              >  1:3877349/1‑150 (MQ=255)
                                        tGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCaa                                                                                                              >  1:4276787/1‑150 (MQ=255)
                                          tACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagc                                                                                                            <  2:2932942/150‑1 (MQ=255)
                                           aCATGAACGCGGTGTAATCCGAGACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagca                                                                                                           >  1:522170/1‑150 (MQ=255)
                                            cATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagcag                                                                                                          >  2:3827471/1‑150 (MQ=255)
                                              tGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagcagca                                                                                                        <  1:4085315/150‑1 (MQ=255)
                                               gAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAagcagcag                                                                                                       <  2:4742944/150‑1 (MQ=255)
                                                        gtAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGCGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCgg                                                                                              <  1:2705473/150‑1 (MQ=255)
                                                         tAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTc                                                                                                                                                <  2:3021489/99‑1 (MQ=255)
                                                         tAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTc                                                                                                                                                >  1:3021489/1‑99 (MQ=255)
                                                           aTCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTc                                                                                           >  1:4696498/1‑150 (MQ=255)
                                                           aTCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTc                                                                                           >  1:4859941/1‑150 (MQ=255)
                                                              cGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCAc                                                                                        <  1:3344226/150‑1 (MQ=255)
                                                               gAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACg                                                                                       >  2:4478728/1‑150 (MQ=255)
                                                                    acGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCgg                                                                                                                                        <  2:2510477/96‑1 (MQ=255)
                                                                    acGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCgg                                                                                                                                        >  1:2510477/1‑96 (MQ=255)
                                                                      gcgcCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATgg                                                                                >  2:4654203/1‑150 (MQ=255)
                                                                            cTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCggg                                                                          <  2:1902299/150‑1 (MQ=255)
                                                                             ttGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGc                                                                         >  2:2031277/1‑150 (MQ=255)
                                                                             ttGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGACAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGc                                                                         >  2:4534442/1‑150 (MQ=255)
                                                                                 tgcATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATg                                                                     >  2:873695/1‑150 (MQ=255)
                                                                                   cATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGAc                                                                   >  1:3353084/1‑150 (MQ=255)
                                                                                   cATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGAc                                                                   >  1:3352723/1‑150 (MQ=255)
                                                                                   cATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGAc                                                                   <  2:3352723/150‑1 (MQ=255)
                                                                                   cATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGAc                                                                   <  2:3353084/150‑1 (MQ=255)
                                                                                     tgttgtGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCa                                                                 <  2:4702819/150‑1 (MQ=255)
                                                                                      gttgtGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAg                                                                >  2:3144301/1‑150 (MQ=255)
                                                                                       ttgtGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGa                                                               <  2:2404817/150‑1 (MQ=255)
                                                                                          tgGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgc                                                            >  1:2616054/1‑150 (MQ=255)
                                                                                             gTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACgctgc                                                         >  1:4920169/1‑150 (MQ=255)
                                                                                              tCACGATGACGATGGTGTACTTGGATTTCAATGCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGAccgaaga                                                        <  2:3417162/150‑7 (MQ=255)
                                                                                                  gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGTTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGctgctgctg                                                    >  2:136646/1‑150 (MQ=255)
                                                                                                  gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGctgctgctg                                                    >  1:3059720/1‑150 (MQ=255)
                                                                                                  gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGctgctgctg                                                    >  2:4799118/1‑150 (MQ=255)
                                                                                                  gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGctgctgctg                                                    <  1:3364830/150‑1 (MQ=255)
                                                                                                  gatgacgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGctgctgctg                                                    >  1:4709934/1‑150 (MQ=255)
                                                                                                       cgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCAc                                               <  1:981195/150‑1 (MQ=255)
                                                                                                       cgatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCAc                                               <  2:4108346/150‑1 (MQ=255)
                                                                                                        gatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCAcc                                              >  2:2217865/1‑150 (MQ=255)
                                                                                                        gatgGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGACAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCAcc                                              >  2:2860959/1‑150 (MQ=255)
                                                                                                            gtgtACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGAc                                          <  1:4946994/150‑1 (MQ=255)
                                                                                                                aCTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGc                                      >  1:1427744/1‑150 (MQ=255)
                                                                                                                aCTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGc                                      >  2:2929668/1‑150 (MQ=255)
                                                                                                                  ttGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGccc                                    >  2:117756/1‑150 (MQ=255)
                                                                                                                     gATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAgc                                 <  1:3154538/150‑1 (MQ=255)
                                                                                                                       tttCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAgcgc                               >  1:3621249/1‑150 (MQ=255)
                                                                                                                           aaTTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGAc                           >  1:571712/1‑150 (MQ=255)
                                                                                                                            aTTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACt                          >  1:4388615/1‑150 (MQ=255)
                                                                                                                             ttCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTc                         <  1:1202046/150‑1 (MQ=255)
                                                                                                                             ttCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTc                         <  1:4092187/150‑1 (MQ=255)
                                                                                                                                aTAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTg                      >  1:46195/1‑150 (MQ=255)
                                                                                                                                aTAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTg                      >  1:368757/1‑150 (MQ=255)
                                                                                                                                    aTCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAgg                  >  2:1290727/1‑150 (MQ=255)
                                                                                                                                    aTCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAgg                  <  1:3972147/150‑1 (MQ=255)
                                                                                                                                    aTCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCAACCGGACAGGCCCACCGCCGCCTCGTGAAgg                  >  2:1525446/1‑150 (MQ=255)
                                                                                                                                        gTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGt              >  2:5347212/1‑150 (MQ=255)
                                                                                                                                        cgatctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGt              >  2:4696811/6‑150 (MQ=255)
                                                                                                                                         ttcttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTc             <  2:4063933/150‑1 (MQ=255)
                                                                                                                                         ttcttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGAGCAGGCGGTc             >  2:4581079/1‑150 (MQ=255)
                                                                                                                                         gatctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCg                                             >  2:223336/5‑118 (MQ=255)
                                                                                                                                         gatctcGACCTTGAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTc             >  2:1814300/5‑150 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGa                                                                       >  1:1606634/1‑89 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGa                                                                       <  2:1606634/89‑1 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGc                 >  1:2819251/1‑143 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt          >  2:4167955/1‑150 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt          >  2:4936249/1‑150 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt          >  2:1681339/1‑150 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt          >  1:1401714/1‑150 (MQ=255)
                                                                                                                                            ttcGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCtt          >  2:822056/1‑150 (MQ=255)
                                                                                                                                               gACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGAc       >  2:2722777/1‑150 (MQ=255)
                                                                                                                                                aCCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGAcc      >  1:3981658/1‑150 (MQ=255)
                                                                                                                                                   ttCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTCg   >  1:4507789/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACAGGGTATCGGTATCCCCGAATTCGACCAGTTTGCCCATGTACATGAACGCGGTGTAATCCGAAACACGCGCCGCTTGCTGCATGTTGTGGGTCACGATGACGATGGTGTACTTGGATTTCAATTCATAGATCAGTTCTTCGACCTTCAAGGTCGAGATCGGGTCCAGCGCCGAGCACGGTTCGTCAAGCAGCAGCACCTCCGGCTCCACGGCAATGGTGCGGGCGATGACCAGACGCTGCTGCTGGCCACCGGACAGGCCCAGCGCCGACTCGTGCAGGCGGTCCTTGACCTCGT  >  NC_002947/6073712‑6074009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: