Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,895,862 G→C 7.1% A157P (GCC→CCC)  PP_3438 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,895,8620GC 7.1% 302.6 / 12.9 98A157P (GCC→CCC) PP_3438membrane protein
Reads supporting (aligned to +/- strand):  ref base G (45/46);  new base C (4/3);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 107 total aligned reads displayed.

GGCACGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAGCGCCCAGGCTAGC  >  NC_002947/3895713‑3896011
                                                                                                                                                      |                                                                                                                                                     
ggCACGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctggg                                                                                                                                                       <  1:182002/150‑1 (MQ=255)
   aCGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctgggcct                                                                                                                                                    >  1:1463361/1‑150 (MQ=255)
     ggCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctgggcctgg                                                                                                                                                  <  2:403100/150‑1 (MQ=255)
        gcAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCtgggcctgggcc                                                                                                                                               <  1:1034503/150‑1 (MQ=255)
                gCAGTGGCTGAACCCCAAAGCCTGGTTGGCCGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTAAATGGATCAAGTCTCTCTGTGATGTCAACTCTGTAGCCTGCTGGGCCTGGGCCGGCAGCGt                                                                                                                                       <  1:46658/150‑1 (MQ=255)
                gCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGt                                                                                                                                       <  1:1130969/150‑1 (MQ=255)
                  aGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGa                                                                                                                                     >  1:2236570/1‑150 (MQ=255)
                   gTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTATGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGAt                                                                                                                                    >  1:274263/1‑150 (MQ=255)
                   gTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATATCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGAt                                                                                                                                    >  1:617270/1‑150 (MQ=255)
                        tGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGcc                                                                                                                               >  2:1066025/1‑150 (MQ=255)
                         gAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCa                                                                                                                              <  2:215780/150‑1 (MQ=255)
                          aaCCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctgggc                                                                                                                                                      >  2:587120/1‑125 (MQ=255)
                          aaCCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctgggc                                                                                                                                                      <  1:587120/125‑1 (MQ=255)
                             cccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTAc                                                                                                                          <  2:587153/150‑1 (MQ=255)
                             cccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTAc                                                                                                                          >  1:1728499/1‑150 (MQ=255)
                             cccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTAc                                                                                                                          >  2:1830328/1‑150 (MQ=255)
                             cccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGAGCCTGGGCCGGCAGCGTGATTCGCCAGTAc                                                                                                                          >  2:621887/1‑150 (MQ=255)
                              ccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACGGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTAcc                                                                                                                         <  1:2191225/150‑1 (MQ=255)
                               cAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCt                                                                                                                        >  1:155534/1‑150 (MQ=255)
                               cAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCt                                                                                                                        <  1:1788313/150‑1 (MQ=255)
                                  aGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGa                                                                                                                     <  1:1806129/150‑1 (MQ=255)
                                    ccTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCGTCTGGGGATCCGTCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACa                                                                                                                   >  2:1881806/1‑150 (MQ=255)
                                     cTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACaa                                                                                                                  >  1:2211344/1‑150 (MQ=255)
                                        gTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGAGCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAAccc                                                                                                               >  2:574602/1‑150 (MQ=255)
                                         ttGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCa                                                                                                              <  2:484210/150‑1 (MQ=255)
                                          tGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCAc                                                                                                             <  1:453651/150‑1 (MQ=255)
                                             cAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGcc                                                                                                          >  1:2015087/1‑150 (MQ=255)
                                              aGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGg                                                                                                                                             <  2:989417‑M2/114‑2 (MQ=255)
                                              aGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGg                                                                                                                                             >  1:989417‑M2/1‑113 (MQ=255)
                                              aGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCa                                                                                                         >  1:1007492/1‑150 (MQ=255)
                                               gCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCAc                                                                                                        <  2:155534/150‑1 (MQ=255)
                                                 cGTTGCCGAGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGc                                                                                                           >  1:1329196/1‑145 (MQ=255)
                                                     gCCGTGGTGGGTGCCTACACCGGCGGTGAACAGAACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTgcgc                                                                                                  >  1:2410734/1‑150 (MQ=255)
                                                       cGGGGTGGTTGGTTACACCGGGGGTGAACAGCCCTGGCGGTGGCTGTTTATTTGTTTGTACGGGCCGTTCTGTCTCAGCTAGGTATCCTGCTGGGACTGGGTCGGCAGAGTGATGTGCCAGGACCTGGACAACCCACGCCACCTGCGCAt                                                                                                <  2:1152923/150‑1 (MQ=255)
                                                              ggtgCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTATGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAAt                                                                                         <  2:274263/150‑1 (MQ=255)
                                                              ggtgCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAAt                                                                                         >  1:1958624/1‑150 (MQ=255)
                                                                tgCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCg                                                                                       >  2:404486/1‑150 (MQ=255)
                                                                    tACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTTCTCTCGGTAGCCTGCCGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGc                                                                                   <  2:835168/150‑1 (MQ=255)
                                                                        ccGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCg                                                                                       <  1:1552316/142‑1 (MQ=255)
                                                                        ccGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCg                                                                                       >  2:1552316/1‑142 (MQ=255)
                                                                         cggcggTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCAACGGGCCGCTCTGTCTCATCTCGGTAGCCTCCTGGGCCTAGGACACAACCGAACTCCAGTAACGCTAGTGGAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGc                                                                              >  1:2380483/1‑150 (MQ=255)
                                                                              gTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCgtgt                                                                         >  1:1370778/1‑150 (MQ=255)
                                                                                gAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCAc                                                                                                             <  1:27322/112‑1 (MQ=255)
                                                                                gAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCAc                                                                                                             >  2:27322/1‑112 (MQ=255)
                                                                                  aCAGCACTTGCTGTGGCTGTTTACCTGGNTCTACGGGCCGATCTNTNTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCCTGCTCAATCGGAGCCTGGCCGTGTNACt                                                                     >  1:1558209/1‑150 (MQ=255)
                                                                                  aCAGCACTTGCTGTGGCTGTTTACCTGGATCTAAGAGCCGAGCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACt                                                                     <  2:2410734/150‑1 (MQ=255)
                                                                                    aGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTgg                                                                   >  1:1913921/1‑150 (MQ=255)
                                                                                      cACTTGGTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTc                                                                 <  2:1016236/150‑1 (MQ=255)
                                                                                      cACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGGTTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTc                                                                 <  2:1651123/150‑1 (MQ=255)
                                                                                      cACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCCTCTACCATCGGAGCCTGGCCGTGTTACTGGTc                                                                 <  1:1881806/150‑1 (MQ=255)
                                                                                      cACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTc                                                                 <  1:404486/150‑1 (MQ=255)
                                                                                       aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                                     >  1:1446687‑M2/1‑72 (MQ=255)
                                                                                       aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                                     <  2:1446687‑M2/97‑26 (MQ=255)
                                                                                       aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                                     <  2:1446588‑M2/97‑26 (MQ=255)
                                                                                       aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                                     <  2:1000720‑M2/97‑26 (MQ=255)
                                                                                       aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                                     >  1:1000720‑M2/1‑72 (MQ=255)
                                                                                       aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                                     >  1:1446588‑M2/1‑72 (MQ=255)
                                                                                           gCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAg                                                            >  1:102473/1‑150 (MQ=255)
                                                                                            cTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAg                                                                                                                             >  1:569287/1‑84 (MQ=255)
                                                                                            cTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAg                                                                                                                             <  2:569287/84‑1 (MQ=255)
                                                                                                ggCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCgg                                                       <  2:2015087/150‑1 (MQ=255)
                                                                                                   tGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTg                                                                                                                       <  1:1583173/83‑1 (MQ=255)
                                                                                                   tGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTg                                                                                                                       >  2:1583173/1‑83 (MQ=255)
                                                                                                      ttACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTAcc                                                 <  1:1403168/150‑1 (MQ=255)
                                                                                                          cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCCCGCCCACTGCGCCTGATCCATCGGTGACTGGCCGTTTTACTGTTCGCTATCGCGGTTTACCTTAt                                             <  2:1728499/150‑1 (MQ=255)
                                                                                                          cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCCCCTCCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGAt                                             >  2:2157639/1‑150 (MQ=255)
                                                                                                          cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCg                                                                >  2:149028/1‑131 (MQ=255)
                                                                                                          cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCg                                                                <  1:149028/131‑1 (MQ=255)
                                                                                                          cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGAt                                             >  1:2205095/1‑150 (MQ=255)
                                                                                                              aTCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCtt                                         <  2:2236570/150‑1 (MQ=255)
                                                                                                                caaCGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCCCGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTa                                       <  2:2380483/148‑1 (MQ=255)
                                                                                                                  aCGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAg                                     <  2:1007492/150‑1 (MQ=255)
                                                                                                                      gCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCtt                                 <  2:1544049/150‑1 (MQ=255)
                                                                                                                         gATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAg                                     >  1:994003/1‑143 (MQ=255)
                                                                                                                           tCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTc                                                                 >  2:112366/1‑113 (MQ=255)
                                                                                                                           tCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTc                                                                 <  1:112366/113‑1 (MQ=255)
                                                                                                                           tCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACg                            >  1:1810936/1‑150 (MQ=255)
                                                                                                                             tGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGc                          >  1:1501928/1‑150 (MQ=255)
                                                                                                                                tttCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTaa                                      <  1:960026/136‑1 (MQ=255)
                                                                                                                                tttCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTaa                                      >  2:960026/1‑136 (MQ=255)
                                                                                                                                tttCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTc                        <  2:1868400/150‑1 (MQ=255)
                                                                                                                                tttCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTc                        >  2:836059/1‑150 (MQ=255)
                                                                                                                                  tCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCt                      <  2:1810936/150‑1 (MQ=255)
                                                                                                                                  tCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCt                      >  2:346702/1‑150 (MQ=255)
                                                                                                                                     tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTg                                <  1:72224/137‑1 (MQ=255)
                                                                                                                                     tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTg                                >  2:72224/1‑137 (MQ=255)
                                                                                                                                     tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCa                   >  1:1648849/1‑150 (MQ=255)
                                                                                                                                     tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCa                   >  2:697005/1‑150 (MQ=255)
                                                                                                                                     tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCGTTAAGGCTTGCACGGCTCCTGCa                   <  1:1066025/150‑1 (MQ=255)
                                                                                                                                         ggTAGCCTGCCGGGCCTGGGCCGGCAGGGTGATTCGGCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAg               <  2:2211344/150‑1 (MQ=255)
                                                                                                                                          gTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGAt                                             >  2:2021300/1‑119 (MQ=255)
                                                                                                                                          gTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGAt                                             <  1:2021300/119‑1 (MQ=255)
                                                                                                                                              ccTGCTGGGCCTGGGCCGgcaggctaccgagatgaaacagatcggcccgtagatccaggtaaacagccacagcaagtgctgttcaccgccggtgtaggcacccaccccggcaacgg                                            <  2:3385‑M2/116‑99 (MQ=255)
                                                                                                                                              ccTGCTGGGCCTGGGCCGgcaggctaccgagatgaaacagatcggcccgtagatccaggtaaacagccacagcaagtgctgttcaccgccggtgtaggcacccaccccggcaacgg                                            >  1:3385‑M2/1‑18 (MQ=255)
                                                                                                                                              ccTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTa                                       <  1:368166/121‑1 (MQ=255)
                                                                                                                                               ctgctgGGCCTGGGCCGgcaggctaccgagatgaaacagatcggcccgtagatccaggtaaacagccacagcaagtgctgttcaccgccggtgtaggcacccaccccggcgtcggctgccaaccaggctttggggttcagccactgcatc         <  1:1706998‑M2/150‑134 (MQ=255)
                                                                                                                                                tgctgGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCgt                                                                           >  2:1978776/1‑83 (MQ=255)
                                                                                                                                                tgctgGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCgt                                                                           <  1:1978776/83‑1 (MQ=255)
                                                                                                                                                      gcctgggccGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAGCGCCCAGGCTAGc  <  2:2447003/150‑1 (MQ=255)
                                                                                                                                                      gcctgggccGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAGCGCCCAGGCTAGc  <  2:1648849/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGCACGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAGCGCCCAGGCTAGC  >  NC_002947/3895713‑3896011

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