Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,628,194 C→T 5.1% intergenic (+38/‑208) PP_5517 → / → PP_2295 hypothetical protein/antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,628,1940CT 5.1% 1323.4 / 36.6 409intergenic (+38/‑208)PP_5517/PP_2295hypothetical protein/antirestriction protein
Reads supporting (aligned to +/- strand):  ref base C (211/177);  new base T (7/14);  total (218/191)
Fisher's exact test for biased strand distribution p-value = 7.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 411 total aligned reads displayed.

AGTTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTACCACATA  >  NC_002947/2628046‑2628343
                                                                                                                                                    |                                                                                                                                                     
aGTTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCg                                                                                                                                                      <  1:7047626/150‑1 (MQ=255)
 gTTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCgg                                                                                                                                                     >  2:5221078/1‑150 (MQ=255)
  ttCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCggg                                                                                                                                                    >  2:554958/1‑150 (MQ=255)
    cGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGcc                                                                                                                                                  <  2:6884770/150‑1 (MQ=255)
     ggggATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGccc                                                                                                                                                 >  1:1753358/1‑150 (MQ=255)
        gATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATTGGGCCCctt                                                                                                                                              <  1:1013714/150‑1 (MQ=255)
        gATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCttt                                                                                                                                              >  2:1849004/1‑150 (MQ=255)
        gATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCttt                                                                                                                                              >  1:5809438/1‑150 (MQ=255)
         aTGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCtttt                                                                                                                                             <  1:5495251/150‑1 (MQ=255)
           gCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCGATTGGGCCCcttgag                                                                                                                                           <  2:5480280/150‑4 (MQ=255)
           gCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTc                                                                                                                                           <  1:6053385/150‑1 (MQ=255)
             gtgtGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTcgc                                                                                                                                         <  2:8249360/150‑1 (MQ=255)
              tgtgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTcgcg                                                                                                                                        >  2:3219983/1‑150 (MQ=255)
                tgGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTc                                                                                                                                      >  2:4288517/1‑150 (MQ=255)
                 ggAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCt                                                                                                                                     <  1:6812638/150‑1 (MQ=255)
                 ggAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCt                                                                                                                                     <  1:2468951/150‑1 (MQ=255)
                 ggAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCt                                                                                                                                     >  2:3603063/1‑150 (MQ=255)
                   aGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGt                                                                                                                                   >  1:815418/1‑150 (MQ=255)
                     gtgtTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTAc                                                                                                                                 <  2:5809569/150‑1 (MQ=255)
                       gtTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGt                                                                                                                               <  2:696753/150‑1 (MQ=255)
                         tCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCACGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCg                                                                                                                             <  1:2259883/150‑1 (MQ=255)
                         tCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCg                                                                                                                             >  1:7328117/1‑150 (MQ=255)
                         tCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACACAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCg                                                                                                                             >  2:7583380/1‑150 (MQ=255)
                           cAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAt                                                                                                                           >  2:4038242/1‑150 (MQ=255)
                           cAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAt                                                                                                                           >  2:5881978/1‑150 (MQ=255)
                            aaCAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAtt                                                                                                                          >  1:5126392/1‑150 (MQ=255)
                            aaCAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAtt                                                                                                                          >  1:3226891/1‑150 (MQ=255)
                            aaCAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGAGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAtt                                                                                                                          <  1:6877612/150‑1 (MQ=255)
                            aaCAGAATGAGTACGCGGCCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGAtt                                                                                                                          >  1:3923426/1‑150 (MQ=255)
                             aCAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTc                                                                                                                         >  1:5960727/1‑150 (MQ=255)
                                 aaTGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGga                                                                                                                     >  1:1526032/1‑150 (MQ=255)
                                 aaTGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGga                                                                                                                     <  1:170503/150‑1 (MQ=255)
                                      gTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAAtt                                                                                                                >  2:5441632/1‑150 (MQ=255)
                                       tACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTc                                                                                                               >  2:4128032/1‑150 (MQ=255)
                                          gcgGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCa                                                                                                            <  2:2405351/150‑1 (MQ=255)
                                          gcgGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCa                                                                                                            >  1:185206/1‑150 (MQ=255)
                                              tCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAgggg                                                                                                        <  1:6402292/150‑1 (MQ=255)
                                                aCCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCTTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCt                                                                                                      >  2:382379/1‑150 (MQ=255)
                                                  cttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTaa                                                                                                    >  2:2679042/1‑150 (MQ=255)
                                                   ttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAg                                                                                                   >  1:8553763/1‑150 (MQ=255)
                                                   ttccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAg                                                                                                   >  1:6135971/1‑150 (MQ=255)
                                                    tccttTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGc                                                                                                  >  1:4442640/1‑150 (MQ=255)
                                                       tttCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCg                                                                                               >  2:8512662/1‑150 (MQ=255)
                                                          ctcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACa                                                                                            <  2:6858211/150‑1 (MQ=255)
                                                          ctcCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACa                                                                                            <  2:6099395/150‑1 (MQ=255)
                                                              gAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAAt                                                                                        <  2:2315831/150‑1 (MQ=255)
                                                              gAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCACTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAAt                                                                                        >  1:3379782/1‑150 (MQ=255)
                                                                  ggTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                                    <  1:7491173/150‑1 (MQ=255)
                                                                  ggTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                                    <  1:7806182/150‑1 (MQ=255)
                                                                  ggTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                                    <  1:5594710/150‑1 (MQ=255)
                                                                  ggTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgc                                                                                    <  1:4650368/150‑1 (MQ=255)
                                                                   gTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCcgcg                                                                                   >  1:2216655/1‑150 (MQ=255)
                                                                       aGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGa                                                                               >  1:8359026/1‑150 (MQ=255)
                                                                       aGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGa                                                                               >  1:1449405/1‑150 (MQ=255)
                                                                        gCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGaa                                                                              <  1:7707846/150‑1 (MQ=255)
                                                                        gCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGaa                                                                              <  2:3267527/150‑1 (MQ=255)
                                                                         ccTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAg                                                                             >  2:4515552/1‑150 (MQ=255)
                                                                            tCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAgcgc                                                                          >  2:4899256/1‑150 (MQ=255)
                                                                              tacacctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAgcgcgc                                                                        >  1:6568998/1‑150 (MQ=255)
                                                                                 acctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCaa                                                                     >  2:5319682/1‑150 (MQ=255)
                                                                                   ctacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCaaaa                                                                   >  2:2621102/1‑150 (MQ=255)
                                                                                    tacaccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAAt                                                                  <  1:273557/150‑1 (MQ=255)
                                                                                      caccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTa                                                                <  2:6608671/150‑1 (MQ=255)
                                                                                      caccAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTa                                                                >  2:2938030/1‑150 (MQ=255)
                                                                                       accAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAc                                                               >  2:1582539/1‑150 (MQ=255)
                                                                                       accAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTAc                                                               >  2:6377383/1‑150 (MQ=255)
                                                                                            gTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGaaaa                                                          >  1:1209126/1‑150 (MQ=255)
                                                                                             ttCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAac                                                         >  2:211758/1‑150 (MQ=255)
                                                                                                 tGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGACGCGCGCCAAAATTACGAAAACACGt                                                     >  2:8610798/1‑150 (MQ=255)
                                                                                                 tGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGt                                                     >  1:1678171/1‑150 (MQ=255)
                                                                                                     gAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCAtt                                                 >  1:2940330/1‑150 (MQ=255)
                                                                                                              aCGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAAt                                        <  1:5346369/150‑1 (MQ=255)
                                                                                                                gCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCACTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACATCATTTCCACAAATgg                                      >  1:4052610/1‑150 (MQ=255)
                                                                                                                 cAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGa                                     >  2:8035021/1‑150 (MQ=255)
                                                                                                                  aaGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaa                                    <  1:412610/150‑1 (MQ=255)
                                                                                                                   aGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaa                                   <  2:770796/150‑1 (MQ=255)
                                                                                                                   aGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaa                                   <  2:70755/150‑1 (MQ=255)
                                                                                                                     tCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaata                                 >  1:3692779/1‑150 (MQ=255)
                                                                                                                     ctaTACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaata                                 <  2:126562/148‑1 (MQ=255)
                                                                                                                        tACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGaaataaat                              <  2:4426395/150‑1 (MQ=255)
                                                                                                                         aCGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATc                             >  1:720804/1‑150 (MQ=255)
                                                                                                                             aaTCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACtt                         >  2:8364947/1‑150 (MQ=255)
                                                                                                                             aaTCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACtt                         >  2:5629695/1‑150 (MQ=255)
                                                                                                                                 acacCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCAt                     >  1:8697943/1‑150 (MQ=255)
                                                                                                                                  cacCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATg                    <  2:1526032/150‑1 (MQ=255)
                                                                                                                                    ccTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg                  <  2:8719601/150‑1 (MQ=255)
                                                                                                                                    ccTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATggg                  <  2:1753358/150‑1 (MQ=255)
                                                                                                                                     cTCCAGGGGCCCACTCGGGCCCTTGTTCGCGTCGGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAAGGGAAATAAATCACTTGCCTgggg                 <  2:8177818/150‑1 (MQ=255)
                                                                                                                                     cTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGAATgggg                 >  1:5418415/1‑150 (MQ=255)
                                                                                                                                      tCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGgt                >  1:651235/1‑150 (MQ=255)
                                                                                                                                       aaaagggCCCGATTGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCAAAAATTACGAAAACACGTCATTTCAACAAATGGAAATAAATCACTTGCATGGGgtg               >  1:1928214/5‑150 (MQ=255)
                                                                                                                                         aGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTc             <  2:2505033/150‑1 (MQ=255)
                                                                                                                                          ggggCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCg            >  1:1345911/1‑150 (MQ=255)
                                                                                                                                           gggCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGt           >  2:7353656/1‑150 (MQ=255)
                                                                                                                                               ccAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGACGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAcc       >  2:3100375/1‑150 (MQ=255)
                                                                                                                                               ccAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAcc       <  1:48077/150‑1 (MQ=255)
                                                                                                                                               ccAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGCGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTAcc       >  2:5120954/1‑150 (MQ=255)
                                                                                                                                                cAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGATTTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTACca      <  2:341288/150‑1 (MQ=255)
                                                                                                                                                cAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTACca      <  2:1361422/150‑1 (MQ=255)
                                                                                                                                                  aTCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTACcaca    >  1:7620787/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGTTCGGGGATGCGTGTGGAGGTGTTCCAACAGAATGAGTACGCGGTCACCTTCCTTTCTCCGAGTGGTCGAGCCTTCTACACCTACACCAAGTTCGTGAGGAAGCGCTGACGCAAGTCATACGCAATCACACCTCAAGGGGCCCAATCGGGCCCTTTTTCGCGTCTGTACGTCGATTCTGGAGAATTCGCAGGGGCTAAGCCCGACAGAATCCGCGTAGAAGCGCGCCAAAATTACGAAAACACGTCATTTCCACAAATGGAAATAAATCACTTGCATGGGGTGTCGTTACCACATA  >  NC_002947/2628046‑2628343

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: