Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,490,910 G→C 5.0% P175R (CCG→CGG) ‡ clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,490,9100GC 5.0% 318.3 / 6.3 100P175R (CCG→CGG) ‡clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base G (47/48);  new base C (2/3);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 103 total aligned reads displayed.

TGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGG  >  NC_002947/3490763‑3491056
                                                                                                                                                    |                                                                                                                                                  
tGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTCCCGCGGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGCCAGCCCAGCggc                                                                                                                                                  <  1:559407/150‑2 (MQ=255)
 gCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGc                                                                                                                                                 >  2:236968/1‑150 (MQ=255)
   cGTACTTCGAACTCACGGCCCAGCACCTGTTTTATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               <  1:1020270/150‑1 (MQ=255)
   cGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               <  2:1328989/150‑1 (MQ=255)
   cGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGAGGCGTCCGCGAGTTTACTCATGGACTGGTCCAGCGCCGGGGTCTTGCCGCTGGACTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               <  1:285358/150‑1 (MQ=255)
     tACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGATCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgc                                                                                                                                             <  1:844737/150‑1 (MQ=255)
     tACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgc                                                                                                                                             <  1:1576971/150‑1 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCTATGCGGCCTTCGCGGGCGCTTTGGGTAGGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGGGCCGGTGGGCTGGCCAGCGGACTGgcctg                                                                                                                                          >  2:2364486/1‑150 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                                          >  1:1110425/1‑150 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                                          >  2:1586983/1‑150 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccag                                                                                                                                          >  1:167267/1‑148 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccag                                                                                                                                          <  2:167267/150‑3 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccag                                                                                                                                          <  2:1380948/150‑3 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTACTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                                          >  1:1512611/1‑150 (MQ=255)
        ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCGCGACGTGCTCACTGGTGCCGCTGGGCTGGCCCGCGGGCTGgcctc                                                                                                                                          >  2:1708558/1‑149 (MQ=255)
           gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCCGCTGGCCAGCCCAGCggcgc                                                                                                                                                >  2:1914811/1‑139 (MQ=255)
           gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCCGCTGGCCAGCCCAGCggcgc                                                                                                                                                <  1:1914811/142‑4 (MQ=255)
            aaCTCACCGCCCAGCACCGGGTCGATGCGTCCTTCGCGGGCGCTTTGGGTAAGGTTGACCTTGTACTGGTCCAGCTCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGcctgcctg                                                                                                                                      <  1:1968945/150‑1 (MQ=255)
               tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctgcccag                                                                                                                                     >  1:1458363/1‑145 (MQ=255)
               tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccagccc                                                                                                                                       <  1:250361/146‑2 (MQ=255)
               tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccagccc                                                                                                                                       >  2:250361/1‑145 (MQ=255)
                 aCGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTc                                                                                                                                 <  1:1428585/150‑1 (MQ=255)
                   ggCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCgg                                                                                                                               <  1:662898/150‑1 (MQ=255)
                      ccAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                            <  1:1586983/150‑1 (MQ=255)
                      ccAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                            >  2:1324444/1‑150 (MQ=255)
                       cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcg                                                                                                                             <  2:1281092/148‑1 (MQ=255)
                       cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcg                                                                                                                             >  1:1281092/1‑148 (MQ=255)
                       cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTCTGGGTAAGGCTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGCCGTGCTCCCCGGCGCCGCTGCGCCGGCCCGCGGGCTGGCCTGCCCGGTCGCGgcgcc                                                                                                                           >  1:1018058/1‑149 (MQ=255)
                        aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                                          <  2:1152510/150‑1 (MQ=255)
                        aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                                          <  1:836655/150‑1 (MQ=255)
                        aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACATGCTCACTCGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                                          >  2:2244913/1‑150 (MQ=255)
                                       cGGCCTTCGCCCGCTCTTTGGCTAAGCTTGACCTTGTACTTGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTCCTCCACGTGCTCTTTCTCGCCGCTGGGCTGGCCAGCGGCCTGGCCTGCCTTTACTCGGCGCTACCCTCCACCAGcgcg                                                                                                           <  2:815483/150‑1 (MQ=255)
                                           cTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCaa                                                                                                       >  1:1460015/1‑150 (MQ=255)
                                             tCGCGGGCGCTTTGGGTCTCGGTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAgg                                                                                                     <  1:1328555/150‑1 (MQ=255)
                                              cgcgGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCTCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGt                                                                                                    <  2:346625/150‑1 (MQ=255)
                                                 gggCGCTTTGGGTACGGTTGACAGGGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCa                                                                                                 <  2:1758810/150‑1 (MQ=255)
                                                 gggCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCa                                                                                                 >  1:2313370/1‑150 (MQ=255)
                                                 gggCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCa                                                                                                 <  2:232607/150‑1 (MQ=255)
                                                       tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCCGGCGCa                                                                                           >  2:538640/1‑150 (MQ=255)
                                                       tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCa                                                                                           >  1:1478079/1‑150 (MQ=255)
                                                       tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCa                                                                                           <  1:1426240/150‑1 (MQ=255)
                                                         tGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGACGATGGAATTGCGCGAGGTGCGCACTGGCACACGTGTGCTGGCCAGCGGGCTTGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAgg                                                                                         <  1:1741155/150‑1 (MQ=255)
                                                             tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                            >  1:1238335/1‑111 (MQ=255)
                                                             tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                            <  2:1238335/111‑1 (MQ=255)
                                                             tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCAt                                                                                     <  2:2025524/150‑1 (MQ=255)
                                                              aaGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATc                                                                                    <  1:951967/150‑1 (MQ=255)
                                                                 gTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                                >  1:1367953/1‑135 (MQ=255)
                                                                 gTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                                <  2:1367953/135‑1 (MQ=255)
                                                                    gACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGtt                                                                              >  2:1859590/1‑150 (MQ=255)
                                                                      ccGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGa                                                                            <  1:333477/150‑1 (MQ=255)
                                                                       cGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAt                                                                           >  1:339286/1‑150 (MQ=255)
                                                                        gtgtACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtt                                                                          >  2:1402228/1‑150 (MQ=255)
                                                                             cTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                                     <  1:236968/150‑1 (MQ=255)
                                                                             cTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                                     <  1:1423559/150‑1 (MQ=255)
                                                                                gTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAg                                                                  >  1:1821614/1‑150 (MQ=255)
                                                                                     gcgcCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCg                                                             <  2:1302446/150‑1 (MQ=255)
                                                                                      cgcCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCgg                                                            >  1:1935511/1‑150 (MQ=255)
                                                                                       gcCGGGGTCTTGCCGCTGGCCTTGCTCGACTTGCTCACTGGCGCCGCTGGGCTGGCCAGAGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCTCGGCAAGGTTCAGGCGCAGTTCATCGGCGTTGATTTTTTCCAGCTCCGGt                                                           <  2:811415/150‑1 (MQ=255)
                                                                                       gcCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGt                                                           >  1:72728/1‑150 (MQ=255)
                                                                                         cGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGc                                                         >  2:66925/1‑150 (MQ=255)
                                                                                         cGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGc                                                         >  1:1695976/1‑150 (MQ=255)
                                                                                           gggTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACGAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCg                                                       >  2:1277774/1‑150 (MQ=255)
                                                                                             gTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAg                                                     >  2:585499/1‑150 (MQ=255)
                                                                                             gTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAg                                                     <  2:1821614/150‑1 (MQ=255)
                                                                                                ttGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCg                                                  >  1:2379933/1‑150 (MQ=255)
                                                                                                  gCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGAt                                                >  2:2270789/1‑150 (MQ=255)
                                                                                                   ccGCTGCCCGTGGTCGGCGGGCTCCCTGGCGCCGCGGGGCTGGCCAGCGGGCTGCCCTGCCTGGACCCGGCGCTACCCCCCCCCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTGGATTTTTTCCAGCGCCGGTGCCGAGGCGATc                                               <  2:1018058/150‑1 (MQ=255)
                                                                                                    cGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               <  2:1016607/53‑1 (MQ=255)
                                                                                                    cGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               >  1:1016607/1‑53 (MQ=255)
                                                                                                                tCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCa                                  <  1:585499/150‑1 (MQ=255)
                                                                                                                  gACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGt                                >  2:186323/1‑150 (MQ=255)
                                                                                                                    cGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAATGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCAcgc                                        >  1:362542/1‑140 (MQ=255)
                                                                                                                    cGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAATGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCAcgc                                        <  2:362542/140‑1 (MQ=255)
                                                                                                                     gTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCa                             <  2:2313370/150‑1 (MQ=255)
                                                                                                                        cTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCa                          >  2:1925996/1‑150 (MQ=255)
                                                                                                                           aCTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                                       >  2:1707088/1‑102 (MQ=255)
                                                                                                                           aCTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                                       <  1:1707088/102‑1 (MQ=255)
                                                                                                                           aCTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCg                       >  1:1578232/1‑150 (MQ=255)
                                                                                                                           aCTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGGTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCg                       >  2:607337/1‑150 (MQ=255)
                                                                                                                            cTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACAATACACCAGCGAGGCAAGGTTACGGCCAAGGTCATCCGCGTTGATTTTTTCTAGCTCCGGTTCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGt                      <  2:1479427/150‑1 (MQ=255)
                                                                                                                               gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGGGCGTTACGCGACGCAGTTCAGCATCGTCaa                   <  1:2244913/150‑1 (MQ=255)
                                                                                                                               gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCaa                   >  2:1300961/1‑150 (MQ=255)
                                                                                                                                cgccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg                  <  2:1578232/150‑1 (MQ=255)
                                                                                                                                cgccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg                  <  1:1324444/150‑1 (MQ=255)
                                                                                                                                     gctgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCaga              <  1:1088096/150‑1 (MQ=255)
                                                                                                                                       tgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGc                                                                                 >  2:348229/1‑81 (MQ=255)
                                                                                                                                       tgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGc                                                                                 <  1:348229/81‑1 (MQ=255)
                                                                                                                                         ggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc                     <  1:1945199/139‑1 (MQ=255)
                                                                                                                                         ggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc                     >  2:1945199/1‑139 (MQ=255)
                                                                                                                                             cgCCATCGGGCGGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgcag      <  2:1110425/149‑1 (MQ=255)
                                                                                                                                               ccAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtttt                                                                        >  2:1898413/1‑82 (MQ=255)
                                                                                                                                               ccAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtttt                                                                        <  1:1898413/82‑1 (MQ=255)
                                                                                                                                                cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGc                                                   <  2:2283643/102‑1 (MQ=255)
                                                                                                                                                cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGc                                                   >  1:2283643/1‑102 (MQ=255)
                                                                                                                                                cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGcagcag   >  1:242740/1‑150 (MQ=255)
                                                                                                                                                cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGcagcag   >  1:926036/1‑150 (MQ=255)
                                                                                                                                                 aGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                                     >  1:1209434/1‑83 (MQ=255)
                                                                                                                                                 aGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                                     <  2:1209434/83‑1 (MQ=255)
                                                                                                                                                 aGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAgg  <  1:875363/150‑1 (MQ=255)
                                                                                                                                                 aGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAgg  >  2:2247733/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGG  >  NC_002947/3490763‑3491056

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