Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,597 G→A 8.0% D209N (GAC→AAC)  PP_0414 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,5970GA 8.0% 551.5 / 27.5 186D209N (GAC→AAC) PP_0414polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (98/73);  new base A (9/6);  total (107/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 188 total aligned reads displayed.

CGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACC  >  NC_002947/503448‑503746
                                                                                                                                                     |                                                                                                                                                     
  gCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGAc                                                                                                                                                     >  1:990668/1‑150 (MQ=255)
  gCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCGTTGCCACCTCGTTTGAc                                                                                                                                                     >  1:4770444/1‑150 (MQ=255)
    gACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGa                                                                                                                                                   >  1:5000920/1‑150 (MQ=255)
       gCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggt                                                                                                                                                >  1:2548365/1‑150 (MQ=255)
         tGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtgg                                                                                                                                              <  2:4112618/150‑1 (MQ=255)
         tGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtgg                                                                                                                                              <  2:2911453/150‑1 (MQ=255)
          gCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggt                                                                                                                                             <  1:1918338/150‑1 (MQ=255)
           cAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtg                                                                                                                                            >  1:75289/1‑150 (MQ=255)
           cAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtg                                                                                                                                            <  1:581656/150‑1 (MQ=255)
            aGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtgg                                                                                                                                           <  1:5087933/150‑1 (MQ=255)
             ggggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtggt                                                                                                                                          >  1:3296560/1‑150 (MQ=255)
             ggggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtaaa                                                                                                                                          <  1:2741589/150‑4 (MQ=255)
              gggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGgtggtggtg                                                                                                                                         >  2:3186368/1‑150 (MQ=255)
              ggcttCAACTACAATCTGGTGCGTGCGGCGGCCAGCCAGGGCGCTTCACCTCTGCTGGCTATCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGgtggtggtg                                                                                                                                         >  2:2053267/4‑150 (MQ=255)
                gTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGAc                                                                                                                                       >  2:2122007/1‑150 (MQ=255)
                gTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGAc                                                                                                                                       <  1:3010939/150‑1 (MQ=255)
                gTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggcaaaggc                                                                                                                                       <  1:3116156/150‑8 (MQ=255)
                 ttCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACg                                                                                                                                      <  1:1054095/150‑1 (MQ=255)
                   cAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCt                                                                                                                                    <  1:3844329/150‑1 (MQ=255)
                    aaCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTg                                                                                                                                   <  2:3387671/150‑1 (MQ=255)
                     aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGt                                                                                                                                  >  2:1137090/1‑150 (MQ=255)
                        aCAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTcc                                                                                                                               >  1:655591/1‑150 (MQ=255)
                          aaTCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTc                                                                                                                             >  2:3369202/1‑150 (MQ=255)
                          aaTCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTc                                                                                                                             <  2:4665321/150‑1 (MQ=255)
                              tGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccg                                                                                                                         <  2:913106/150‑1 (MQ=255)
                                 tgcgtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggc                                                                                                                      <  1:4216698/150‑1 (MQ=255)
                                 tgcgtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggc                                                                                                                      <  1:2521937/150‑1 (MQ=255)
                                       cggcggCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGc                                                                                                                <  1:4137734/150‑1 (MQ=255)
                                           ggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGc                                                                                                            <  1:156836/150‑1 (MQ=255)
                                              cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAac                                                                                                         >  1:4627429/1‑150 (MQ=255)
                                              cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAac                                                                                                         >  2:4020431/1‑150 (MQ=255)
                                              cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAac                                                                                                         >  1:686167/1‑150 (MQ=255)
                                               aGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAaca                                                                                                        >  1:2617700/1‑150 (MQ=255)
                                               aGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggt                                                                                                                                                >  2:2506943/1‑110 (MQ=255)
                                                gCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAacac                                                                                                       >  1:5148317/1‑150 (MQ=255)
                                                gCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAacac                                                                                                       >  1:1679483/1‑150 (MQ=255)
                                                gCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAacac                                                                                                       >  1:223120/1‑150 (MQ=255)
                                                   tGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGc                                                                                                    <  2:3540418/150‑1 (MQ=255)
                                                    gggCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGcc                                                                                                   >  2:4481248/1‑150 (MQ=255)
                                                      gcgcTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCcgc                                                                                                 <  1:4553089/150‑1 (MQ=255)
                                                         cTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGc                                                                                              <  1:3980982/150‑1 (MQ=255)
                                                          ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCa                                                                                             >  1:4684500/1‑150 (MQ=255)
                                                          ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCa                                                                                             >  2:81780/1‑150 (MQ=255)
                                                          ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTATTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCa                                                                                             <  2:610565/150‑1 (MQ=255)
                                                             aCCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGAt                                                                                          >  2:4847579/1‑150 (MQ=255)
                                                             aCCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGAt                                                                                          <  1:2499088/150‑1 (MQ=255)
                                                              ccTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATg                                                                                         >  1:4126188/1‑150 (MQ=255)
                                                               ctctGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGt                                                                                        >  2:3918982/1‑150 (MQ=255)
                                                                tctGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGtt                                                                                       >  1:470517/1‑150 (MQ=255)
                                                                 ctgctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGCGCCCTGTTTGCGTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGCGCAGGCAACACTGCCGCGGCCGGTGTTc                                                                                      >  1:4283873/1‑150 (MQ=255)
                                                                 ctgctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTc                                                                                      >  1:2826725/1‑150 (MQ=255)
                                                                     tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGCCGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                                  <  2:666376/150‑1 (MQ=255)
                                                                     tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                                  >  2:2578247/1‑150 (MQ=255)
                                                                     tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                                  >  2:229866/1‑150 (MQ=255)
                                                                     tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                                  >  2:1533844/1‑150 (MQ=255)
                                                                        cTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTa                                                                               >  2:3121743/1‑150 (MQ=255)
                                                                         tttCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTAt                                                                              >  1:795172/1‑150 (MQ=255)
                                                                          ttCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTAtt                                                                             >  1:3526343/1‑150 (MQ=255)
                                                                           tCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTc                                                                            <  1:1329671/150‑1 (MQ=255)
                                                                              gccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTcgcg                                                                         >  1:5197297/1‑150 (MQ=255)
                                                                               ccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                        >  1:1512991/1‑150 (MQ=255)
                                                                               ccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                        <  2:838612/150‑1 (MQ=255)
                                                                                  ggTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGaa                                                                     <  1:1137090/150‑1 (MQ=255)
                                                                                          tGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGcc                                                             >  2:2265030/1‑150 (MQ=255)
                                                                                            ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCg                                                           >  1:96023/1‑150 (MQ=255)
                                                                                             cGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGa                                                          >  1:2672514/1‑150 (MQ=255)
                                                                                             cGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGa                                                          >  1:4080769/1‑150 (MQ=255)
                                                                                             cGCTGAGTGCGCCCGGTGTCATCGCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGa                                                          <  2:4381665/150‑1 (MQ=255)
                                                                                              gCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGAc                                                         <  2:5000920/150‑1 (MQ=255)
                                                                                                            gtgtCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgcggcgg                                           >  1:1563008/1‑150 (MQ=255)
                                                                                                                cATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAc                                       <  1:619989/150‑1 (MQ=255)
                                                                                                                 aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                      >  2:3952104/1‑150 (MQ=255)
                                                                                                                  tctcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCcg                                                                                                  >  2:980490/1‑89 (MQ=255)
                                                                                                                     cGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGc                                  <  2:470517/150‑1 (MQ=255)
                                                                                                                      gggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAc                                       <  2:1548579/144‑1 (MQ=255)
                                                                                                                        gAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGa                               >  2:53021/1‑150 (MQ=255)
                                                                                                                          gCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGAtt                             <  1:3186368/150‑1 (MQ=255)
                                                                                                                             cTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGc                                                                                              >  2:639349/1‑82 (MQ=255)
                                                                                                                             cTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                        <  1:5369087/104‑1 (MQ=255)
                                                                                                                              tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGt                         >  2:3789303/1‑150 (MQ=255)
                                                                                                                                   gcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg                    <  2:231710/150‑1 (MQ=255)
                                                                                                                                   gcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg                    <  1:229866/150‑1 (MQ=255)
                                                                                                                                   gcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg                    >  2:4249103/1‑150 (MQ=255)
                                                                                                                                   gcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg                    <  1:4847579/150‑1 (MQ=255)
                                                                                                                                    cctttgccACCTCGTTTGACGGGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCgg                   <  1:4772634/150‑1 (MQ=255)
                                                                                                                                     ctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgt                  <  1:1533844/150‑1 (MQ=255)
                                                                                                                                      tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtg                 >  1:5022604/1‑150 (MQ=255)
                                                                                                                                        caccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgtt               >  1:1071109/3‑150 (MQ=255)
                                                                                                                                        caccACCTCGTCAAACGAGGTGGCGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgtt               >  2:1136430/3‑150 (MQ=255)
                                                                                                                                           cACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTgttgtt            >  2:5337477/1‑150 (MQ=255)
                                                                                                                                           cACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTgttgtt            >  2:2880888/1‑150 (MQ=255)
                                                                                                                                            aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGttgttg           <  1:2122007/150‑1 (MQ=255)
                                                                                                                                            aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGACAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGttgttg           >  1:2637083/1‑150 (MQ=255)
                                                                                                                                              cTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgct         <  2:1275368/150‑1 (MQ=255)
                                                                                                                                              cTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgct         <  1:1280160/150‑1 (MQ=255)
                                                                                                                                               tCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGTGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctg        >  1:5066104/1‑150 (MQ=255)
                                                                                                                                               tCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctg        >  1:1188739/1‑150 (MQ=255)
                                                                                                                                                cGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctgc       <  2:4529266/150‑1 (MQ=255)
                                                                                                                                                cGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctgc       <  2:1787707/150‑1 (MQ=255)
                                                                                                                                                 gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTgctgct      >  1:2873186/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACC  >  NC_002947/503448‑503746

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: