Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 296,486 C→T 5.5% G507S (GGT→AGT) ‡ gshA ← glutamate‑‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947296,4860CT 5.5% 1208.7 / 35.5 379G507S (GGT→AGT) ‡gshAglutamate‑‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (187/171);  new base T (12/9);  total (199/180)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 390 total aligned reads displayed.

CGGCCAACAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCT  >  NC_002947/296337‑296635
                                                                                                                                                     |                                                                                                                                                     
cGGCCAACAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGATCAGGTCGAAGTCAcc                                                                                                                                                       <  2:1351308/150‑1 (MQ=255)
 ggCCAACAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCg                                                                                                                                                      <  1:4668334/150‑1 (MQ=255)
   ccAACAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGAc                                                                                                                                                    <  2:219914/150‑1 (MQ=255)
    cAACAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACt                                                                                                                                                   >  1:4523406/1‑150 (MQ=255)
      aCAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCTGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTc                                                                                                                                                 <  1:3863036/150‑1 (MQ=255)
         gCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGTGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCGGTGACTTCgac                                                                                                                                              <  1:755388/150‑3 (MQ=255)
         gCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGtt                                                                                                                                              <  1:1945510/150‑1 (MQ=255)
         gCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTGTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCGGTGACTTCgac                                                                                                                                              <  1:699010/150‑3 (MQ=255)
         gCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTGGGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCGGTGACTTCgac                                                                                                                                              <  1:3454410/150‑3 (MQ=255)
          cTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                                                                                                                                             <  1:957940/150‑1 (MQ=255)
          cTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGttt                                                                                                                                             <  1:2526439/150‑1 (MQ=255)
            cACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTg                                                                                                                                           <  1:3963666/150‑1 (MQ=255)
             aCGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGc                                                                                                                                          >  1:7305938/1‑150 (MQ=255)
             aCGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCCGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGc                                                                                                                                          <  2:8592370/150‑1 (MQ=255)
               gcgcAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTc                                                                                                                                        >  2:7018416/1‑150 (MQ=255)
               gcgcAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTc                                                                                                                                        <  1:2758919/150‑1 (MQ=255)
                 gcAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCa                                                                                                                                      <  1:5150469/150‑1 (MQ=255)
                      cgccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCg                                                                                                                                 <  2:4796878/150‑1 (MQ=255)
                      cgccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCg                                                                                                                                 <  2:3917327/150‑1 (MQ=255)
                       gccAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCg                                                                                                                                 >  2:4100267/1‑149 (MQ=255)
                       gccAGGATCGAGGCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCgg                                                                                                                                >  1:3326321/1‑150 (MQ=255)
                            gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc                                                                                                                           >  2:8678285/1‑150 (MQ=255)
                            gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTGTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc                                                                                                                           <  1:684784/150‑1 (MQ=255)
                            gATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTGGGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGc                                                                                                                           <  1:7642280/150‑1 (MQ=255)
                                 ccTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAGCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                      >  2:674122/1‑150 (MQ=255)
                                 ccTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                      <  1:1770711/150‑1 (MQ=255)
                                 ccTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                      <  2:3041789/150‑1 (MQ=255)
                                 aggcGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                      >  2:1643644/4‑150 (MQ=255)
                                 aggcGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgc                                                                                                                      >  2:6262852/4‑150 (MQ=255)
                                  cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcg                                                                                                                     >  2:2936223/1‑150 (MQ=255)
                                  cTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcg                                                                                                                     >  2:3387258/1‑150 (MQ=255)
                                   tCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcga                                                                                                                    <  1:4149501/150‑1 (MQ=255)
                                    cGATCCTGGCCTGGCTGCTGGCGCGTGGTGTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcgag                                                                                                                   <  1:7143720/150‑1 (MQ=255)
                                       tCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGgcgagcga                                                                                                                >  1:7981613/1‑150 (MQ=255)
                                               tggcTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCggg                                                                                                        <  2:5928545/150‑1 (MQ=255)
                                                ggcTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGc                                                                                                       >  1:5335621/1‑150 (MQ=255)
                                                      tgGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTg                                                                                                                                           >  2:5501033/1‑108 (MQ=255)
                                                      tgGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTg                                                                                                                                           <  1:5501033/108‑1 (MQ=255)
                                                        gcgcGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGtct                                                                                               >  1:2830204/1‑150 (MQ=255)
                                                           cgTGGTGTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGa                                                                                            <  1:5039451/150‑1 (MQ=255)
                                                              ggTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAAtg                                                                                         <  1:2067797/150‑1 (MQ=255)
                                                                 ttGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAAtgctt                                                                                      <  2:2476090/150‑1 (MQ=255)
                                                                   gcgcGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAAtgcttgc                                                                                    <  1:5256403/150‑1 (MQ=255)
                                                                       ggggTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGGTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCCGCGGCTCGAATGCTTgctgct                                                                                >  2:4201439/1‑150 (MQ=255)
                                                                           tggtggtgCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTc                                                                            <  1:4862995/150‑1 (MQ=255)
                                                                           agcaggtgCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTc                                                                            <  1:3199864/146‑1 (MQ=255)
                                                                                  gCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTc                                                                     >  2:4154695/1‑150 (MQ=255)
                                                                                  gCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTc                                                                     >  1:8573332/1‑150 (MQ=255)
                                                                                     gTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGc                                                                  >  2:6412937/1‑150 (MQ=255)
                                                                                         atcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGtt                                                              >  1:1517276/1‑150 (MQ=255)
                                                                                         atcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGtt                                                              >  2:3490003/1‑150 (MQ=255)
                                                                                          tcaACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTg                                                             >  1:2990337/1‑150 (MQ=255)
                                                                                           caACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGc                                                            <  1:1071349/150‑1 (MQ=255)
                                                                                           caACTGCTGATTGCCAGGATGCGGGCCTGGTAGGCGCCGACGACCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGc                                                            <  1:8195363/150‑1 (MQ=255)
                                                                                                gctgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTcgcg                                                       <  1:5963506/150‑1 (MQ=255)
                                                                                                  tgATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGa                                                     >  2:4428575/1‑150 (MQ=255)
                                                                                                   gATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTtgctgc                                                                                 >  2:8153988/1‑121 (MQ=255)
                                                                                                    aTTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGaaa                                                   <  2:4468991/150‑1 (MQ=255)
                                                                                                        ccAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGtct                                               >  2:1533829/1‑150 (MQ=255)
                                                                                                         cAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGtctc                                              >  2:5424641/1‑150 (MQ=255)
                                                                                                         cAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGtctc                                              >  1:1418586/1‑150 (MQ=255)
                                                                                                         cAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGtctc                                              <  1:6228886/150‑1 (MQ=255)
                                                                                                         cAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGtctc                                              <  2:1011650/150‑1 (MQ=255)
                                                                                                         cAGGATGCTGGCCACGTAGCCGCCGACGAACACGTCCAAGTCACCGACTTCGTTTAGCTCCAGCTCGGACTGCCAGGCGAGCCACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGATGCTCAGCGCATGCACGCGCAAAGtctc                                              >  2:1881092/1‑150 (MQ=255)
                                                                                                            gATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGc                                           >  2:2735880/1‑150 (MQ=255)
                                                                                                              tGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGt                                         <  2:8769857/150‑1 (MQ=255)
                                                                                                              tGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGt                                         <  2:5147082/150‑1 (MQ=255)
                                                                                                               gCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTg                                        >  2:6487138/1‑150 (MQ=255)
                                                                                                                 tGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGca                                      >  2:5630947/1‑150 (MQ=255)
                                                                                                                   gCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGcaca                                    >  1:1564864/1‑150 (MQ=255)
                                                                                                                   gCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGcaca                                    <  1:2352962/150‑1 (MQ=255)
                                                                                                                     cTGATAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACg                                  >  2:1272848/1‑150 (MQ=255)
                                                                                                                      tGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACgg                                 <  1:4961236/150‑1 (MQ=255)
                                                                                                                         tAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGctc                              >  2:7093378/1‑150 (MQ=255)
                                                                                                                             cgcCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGgc                          <  1:2804839/150‑1 (MQ=255)
                                                                                                                             cgcCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGgc                          <  2:7556672/150‑1 (MQ=255)
                                                                                                                                cgacgaACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCa                       <  2:6906419/150‑1 (MQ=255)
                                                                                                                                cgacgaACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCa                       <  1:8267646/150‑1 (MQ=255)
                                                                                                                                  acgaACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGctct                             >  2:2004381/1‑142 (MQ=255)
                                                                                                                                  acgaACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGc                     >  2:2355939/1‑150 (MQ=255)
                                                                                                                                  acgaACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGc                     <  1:8737111/150‑1 (MQ=255)
                                                                                                                                     aaCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGTGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgag                  <  2:367600/150‑1 (MQ=255)
                                                                                                                                      aCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgaga                 >  2:5715363/1‑150 (MQ=255)
                                                                                                                                      aCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgaga                 >  1:4254188/1‑150 (MQ=255)
                                                                                                                                      aCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgaga                 >  2:3950016/1‑150 (MQ=255)
                                                                                                                                      aCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgaga                 >  2:7864679/1‑150 (MQ=255)
                                                                                                                                      aCAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCgaga                 >  2:7650305/1‑150 (MQ=255)
                                                                                                                                       cAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGaa                >  2:719589/1‑150 (MQ=255)
                                                                                                                                         ggTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACt              <  2:7321120/150‑1 (MQ=255)
                                                                                                                                         aaacGAAGTCGGTGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACt              >  1:7133785/4‑150 (MQ=255)
                                                                                                                                          gTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTg             <  2:5822728/150‑1 (MQ=255)
                                                                                                                                          aacGAAGTCGGTGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTg             >  1:6605721/3‑150 (MQ=255)
                                                                                                                                            cGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGa           <  1:561567/150‑1 (MQ=255)
                                                                                                                                            cGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGa           <  1:546817/150‑1 (MQ=255)
                                                                                                                                            cGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGa           <  1:4803821/150‑1 (MQ=255)
                                                                                                                                                gTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGaa       <  1:3702126/150‑1 (MQ=255)
                                                                                                                                                 tCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAg      >  2:5020178/1‑150 (MQ=255)
                                                                                                                                                 tCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAg      >  2:1779026/1‑150 (MQ=255)
                                                                                                                                                   aCCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGct    >  1:2919329/1‑150 (MQ=255)
                                                                                                                                                     cGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGctct  >  1:749226/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCCAACAGCTCACGCGCAGCCGCCAGGATCGCCTCGATCCTGGCCTGGCTGCTGGCGCGTGGTTTGCGCGGGGTGGTGGTGCGGTTCATCAACTGCTGATTGCCAGGATGCTGGCCTGGTAGGCGCCGACGAACAGGTCGAAGTCACCGACTTCGTTTTGCTCCAGCTCGGACTGCCTGGCGAGCGACTCGCGGGCCAGCGTCTCGAATGCTTGCTGCTTTTCGTTGCTCAGCGGTTGCTCGCGGAAAGTCTCGGCGTGCACACGGCTCTGGCGCAGCGAGAACTGGACGAAGCTCT  >  NC_002947/296337‑296635

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