Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,895,859 T→C 8.8% W156R (TGG→CGG)  PP_3438 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,895,8590TC 8.8% 264.1 / 9.9 95W156R (TGG→CGG) PP_3438membrane protein
Reads supporting (aligned to +/- strand):  ref base T (44/42);  new base C (4/5);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 107 total aligned reads displayed.

CCTGGCACGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAGCGCCCA  >  NC_002947/3895710‑3896004
                                                                                                                                                      |                                                                                                                                                 
ccTGGCACGGCGCAGAGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCctgct                                                                                                                                                   >  2:1403168/1‑150 (MQ=255)
   ggCACGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctggg                                                                                                                                                <  1:182002/150‑1 (MQ=255)
      aCGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCctgct                                                                                                                                                   >  2:661771/1‑144 (MQ=255)
      aCGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCctgct                                                                                                                                                   <  1:661771/144‑1 (MQ=255)
      aCGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGctgggcct                                                                                                                                             >  1:1463361/1‑150 (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:1830328/1‑150 (MQ=255)
                                cccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTAc                                                                                                                   >  1:1728499/1‑150 (MQ=255)
                                cccAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTAc                                                                                                                   <  2:587153/150‑1 (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)
                                                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                                                                                >  2:1552316/1‑142 (MQ=255)
                                                                           ccGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCg                                                                                <  1:1552316/142‑1 (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:1000720‑M2/1‑72 (MQ=255)
                                                                                          aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                              >  1:1446588‑M2/1‑72 (MQ=255)
                                                                                          aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                              <  2:1000720‑M2/97‑26 (MQ=255)
                                                                                          aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                              >  1:1446687‑M2/1‑72 (MQ=255)
                                                                                          aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                              <  2:1446588‑M2/97‑26 (MQ=255)
                                                                                          aCTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCCGGCCCAGGCCCAgcaggctaccgagatgaaacagatc                                                                                                              <  2:1446687‑M2/97‑26 (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                                                                                                                >  2:1583173/1‑83 (MQ=255)
                                                                                                      tGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTg                                                                                                                <  1:1583173/83‑1 (MQ=255)
                                                                                                         ttACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTAcc                                          <  1:1403168/150‑1 (MQ=255)
                                                                                                             cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCCCGCCCACTGCGCCTGATCCATCGGTGACTGGCCGTTTTACTGTTCGCTATCGCGGTTTACCTTAt                                      <  2:1728499/150‑1 (MQ=255)
                                                                                                             cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCCCCTCCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGAt                                      >  2:2157639/1‑150 (MQ=255)
                                                                                                             cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCg                                                         <  1:149028/131‑1 (MQ=255)
                                                                                                             cTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCg                                                         >  2:149028/1‑131 (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:836059/1‑150 (MQ=255)
                                                                                                                                   tttCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTc                 <  2:1868400/150‑1 (MQ=255)
                                                                                                                                     tCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCt               <  2:1810936/150‑1 (MQ=255)
                                                                                                                                     tCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCt               >  2:346702/1‑150 (MQ=255)
                                                                                                                                        tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTg                         >  2:72224/1‑137 (MQ=255)
                                                                                                                                        tctcGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTg                         <  1:72224/137‑1 (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                                     >  1:3385‑M2/1‑18 (MQ=255)
                                                                                                                                                 ccTGCTGGGCCTGGGCCGgcaggctaccgagatgaaacagatcggcccgtagatccaggtaaacagccacagcaagtgctgttcaccgccggtgtaggcacccaccccggcaacgg                                     <  2:3385‑M2/116‑99 (MQ=255)
                                                                                                                                                 ccTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTa                                <  1:368166/121‑1 (MQ=255)
                                                                                                                                                 ccTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTa                                >  2:368166/1‑121 (MQ=255)
                                                                                                                                                  ctgctgGGCCTGGGCCGgcaggctaccgagatgaaacagatcggcccgtagatccaggtaaacagccacagcaagtgctgttcaccgccggtgtaggcacccaccccggcgtcggctgccaaccaggctttggggttcagccactgcatc  <  1:1706998‑M2/150‑134 (MQ=255)
                                                                                                                                                   tgctgGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCgt                                                                    >  2:1978776/1‑83 (MQ=255)
                                                                                                                                                   tgctgGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCgt                                                                    <  1:1978776/83‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CCTGGCACGGCGCAGCGATGCAGTGGCTGAACCCCAAAGCCTGGTTGGCAGCCGTTGCCGGGGTGGGTGCCTACACCGGCGGTGAACAGCACTTGCTGTGGCTGTTTACCTGGATCTACGGGCCGATCTGTTTCATCTCGGTAGCCTGCTGGGCCTGGGCCGGCAGCGTGATTCGCCAGTACCTGGACAACCCACGCCACCTGCGCATGCTCAATCGGAGCCTGGCCGTGTTACTGGTCGCCAGCGCGGTTTACCTGATCCTTTAAGGCTTGCACGGCTCCTGCAGGAGCGCCCA  >  NC_002947/3895710‑3896004

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