Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,044,330 C→G 5.7% A223A (GCC→GCG bglA → 6‑phospho‑beta‑glucosidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,044,3300CG 5.7% 513.3 / 14.6 154A223A (GCC→GCGbglA6‑phospho‑beta‑glucosidase A
Reads supporting (aligned to +/- strand):  ref base C (67/78);  new base G (4/5);  total (71/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 157 total aligned reads displayed.

GGATGACCTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCTTTACCGATG  >  NC_000913/3044185‑3044478
                                                                                                                                                 |                                                                                                                                                    
 gATGACCTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGc                                                                                                                                                 <  2:1432805/150‑1 (MQ=255)
  aTGACCTTCAACGAAATTAACACCCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGGGCGCTGGCg                                                                                                                                                <  1:1162984/150‑1 (MQ=255)
  aTGACCTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCg                                                                                                                                                <  1:824046/150‑1 (MQ=255)
     aCCTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTg                                                                                                                                             <  1:1141776/150‑1 (MQ=255)
     aCCTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTg                                                                                                                                             <  1:1682309/150‑1 (MQ=255)
      ccTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAGGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGa                                                                                                                                            <  1:1726835/150‑1 (MQ=255)
      ccTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGGGCGCTGGCgacaa                                                                                                                                            <  2:1727680/150‑5 (MQ=255)
       cTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGaa                                                                                                                                           >  1:2132638/1‑150 (MQ=255)
        ttCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTTCTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGaaa                                                                                                                                          >  2:43639/1‑150 (MQ=255)
        ttCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGaaa                                                                                                                                          <  2:1969147/150‑1 (MQ=255)
         tCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAg                                                                                                                                         <  1:157135/150‑1 (MQ=255)
         tCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGCTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGaaac                                                                                                                                         >  2:284050/1‑149 (MQ=255)
          cAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGc                                                                                                                                        >  1:2260786/1‑150 (MQ=255)
           aaCGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCt                                                                                                                                       <  1:279573/150‑1 (MQ=255)
           aaCGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCt                                                                                                                                       <  2:1186323/150‑1 (MQ=255)
           aaCGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCCTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCCGCGCCCTGGCGGTGAAAGCt                                                                                                                                       >  1:1062316/1‑150 (MQ=255)
            aCGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTg                                                                                                                                      >  1:2416319/1‑150 (MQ=255)
             cGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTgc                                                                                                                                     <  2:458503/150‑1 (MQ=255)
                  ttAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGGGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTc                                                                                                                                >  2:2240228/1‑150 (MQ=255)
                   tAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGAATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTcc                                                                                                                               >  2:1578382/1‑149 (MQ=255)
                   tAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGGGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCg                                                                                                                               <  1:1674843/150‑1 (MQ=255)
                      caacCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCAt                                                                                                                            >  2:2355515/1‑150 (MQ=255)
                             cGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCg                                                                                                                     >  1:278370/1‑150 (MQ=255)
                              gTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCgg                                                                                                                    <  1:602900/150‑1 (MQ=255)
                              gTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCgg                                                                                                                    <  1:886466/150‑1 (MQ=255)
                                aaCTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGGGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGgag                                                                                                                  <  1:2366427/150‑1 (MQ=255)
                                 aCTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGgaga                                                                                                                 >  1:412031/1‑150 (MQ=255)
                                   tGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATg                                                                                                               <  2:2206010/150‑1 (MQ=255)
                                    ggCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGa                                                                                                              >  1:222667/1‑150 (MQ=255)
                                     gCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGaa                                                                                                             >  2:359505/1‑150 (MQ=255)
                                       gTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAg                                                                                                           <  2:845824/150‑1 (MQ=255)
                                          cACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCg                                                                                                        >  1:2190490/1‑150 (MQ=255)
                                          cACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCg                                                                                                        <  1:1525443/150‑1 (MQ=255)
                                            ccGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGc                                                                                                      >  2:22948/1‑150 (MQ=255)
                                              gCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTg                                                                                                    >  2:927599/1‑150 (MQ=255)
                                                tGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTa                                                                                                  >  2:309478/1‑150 (MQ=255)
                                                     ggTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTg                                                                                             >  1:238235/1‑150 (MQ=255)
                                                      gTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTgg                                                                                            >  2:1662617/1‑150 (MQ=255)
                                                      gTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTgg                                                                                            <  2:378120/150‑1 (MQ=255)
                                                           tgctgcTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATg                                                                                       <  2:491894/150‑1 (MQ=255)
                                                              tgctCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTg                                                                                    <  2:698610/150‑1 (MQ=255)
                                                                       gtggtgTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTAt                                                                           <  2:2190808/150‑1 (MQ=255)
                                                                          gtgtATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCt                                                                        <  2:1193850/150‑1 (MQ=255)
                                                                            gtATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTa                                                                      >  2:218355/1‑150 (MQ=255)
                                                                             tataCCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTAc                                                                     <  1:1334341/150‑1 (MQ=255)
                                                                             tataCCGAGCATGAAAACCCGGAAGAGACGATATATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTAc                                                                     <  2:1229457/150‑1 (MQ=255)
                                                                              ataCCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGCAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACt                                                                    >  2:1311008/1‑150 (MQ=255)
                                                                                 ccGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCt                                                                 >  1:1684538/1‑150 (MQ=255)
                                                                                 ccGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCt                                                                 >  2:315623/1‑150 (MQ=255)
                                                                                  cGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCGGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTg                                                                <  2:717947/150‑1 (MQ=255)
                                                                                    aGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTa                                                              <  1:43639/150‑1 (MQ=255)
                                                                                     gCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGGCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTaa                                                             <  1:840126/150‑1 (MQ=255)
                                                                                      cATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAAc                                                            <  1:296674/150‑1 (MQ=255)
                                                                                       aTGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAAcc                                                           <  2:112587/150‑1 (MQ=255)
                                                                                        tGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAAccc                                                          >  2:185059/1‑150 (MQ=255)
                                                                                        tGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAAccc                                                          >  2:1620356/1‑150 (MQ=255)
                                                                                        tGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACAAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTCTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAAccc                                                          >  1:2007407/1‑150 (MQ=255)
                                                                                         gAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCg                                                         >  1:2149012/1‑150 (MQ=255)
                                                                                             aCCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACg                                                     <  2:2313362/150‑1 (MQ=255)
                                                                                                   aaGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACgatgtga                                               <  1:943586/150‑1 (MQ=255)
                                                                                                    agagACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACgatgtgat                                              >  2:846459/1‑150 (MQ=255)
                                                                                                     tagaCGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACgatgtgatg                                             <  1:2399834/149‑1 (MQ=255)
                                                                                                         cGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCg                                         <  2:1707128/150‑1 (MQ=255)
                                                                                                          gATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGc                                        <  1:1690230/150‑1 (MQ=255)
                                                                                                            tGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTc                                      <  1:2355515/150‑1 (MQ=255)
                                                                                                            acaaactggTGATGCAGCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTc                                      >  2:2097827/8‑150 (MQ=255)
                                                                                                             gTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCa                                     >  1:1608823/1‑150 (MQ=255)
                                                                                                             gTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCa                                     <  2:2282825/150‑1 (MQ=255)
                                                                                                              tATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCTGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCGGGACTATGTGATGTTCGCTCAg                                    <  2:2046693/150‑1 (MQ=255)
                                                                                                              tATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAg                                    <  1:1196039/150‑1 (MQ=255)
                                                                                                              tATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAg                                    <  1:1662617/150‑1 (MQ=255)
                                                                                                               aTCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAgg                                   <  1:676506/150‑1 (MQ=255)
                                                                                                               aTCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAgg                                   <  1:711572/150‑1 (MQ=255)
                                                                                                               aTCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAgg                                   >  1:982826/1‑150 (MQ=255)
                                                                                                               aTCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCGCatc                                   >  1:1048080/1‑148 (MQ=255)
                                                                                                                  aGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGt                                >  2:1428074/1‑150 (MQ=255)
                                                                                                                   ggTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTc                               >  1:2068679/1‑150 (MQ=255)
                                                                                                                       ctgcATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATg                           <  1:22948/150‑1 (MQ=255)
                                                                                                                        tgcATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATgc                          <  1:1413693/150‑1 (MQ=255)
                                                                                                                        tgcATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATgc                          >  1:1710933/1‑150 (MQ=255)
                                                                                                                               ccAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATgg                                                                                      <  1:1343069/83‑1 (MQ=255)
                                                                                                                               ccAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATgg                                                                                      >  2:1343069/1‑83 (MQ=255)
                                                                                                                               ccAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACg                   <  2:2121855/150‑1 (MQ=255)
                                                                                                                               ccAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACg                   >  2:1527467/1‑150 (MQ=255)
                                                                                                                                 ctttcaccGCCAGGGCGCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCt                 >  2:1515775/8‑150 (MQ=255)
                                                                                                                                  tttcaccGCCAGGGCGCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTa                >  2:612219/7‑150 (MQ=255)
                                                                                                                                   tttGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGCTCGCTCAGGAGTCGATGCGCGAACGCTAc               <  2:1062316/150‑1 (MQ=255)
                                                                                                                                    ttGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACg              <  1:1548446/150‑1 (MQ=255)
                                                                                                                                    tcaccGCCAGGGCGCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACg              >  1:1244166/5‑150 (MQ=255)
                                                                                                                                     tGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGt                                >  2:1257178/1‑131 (MQ=255)
                                                                                                                                     tGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGt                                <  1:1257178/131‑1 (MQ=255)
                                                                                                                                      gTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTc            >  1:2004074/1‑150 (MQ=255)
                                                                                                                                      gTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTc            <  1:309478/150‑1 (MQ=255)
                                                                                                                                        cGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCtt          >  2:655811/1‑150 (MQ=255)
                                                                                                                                         gCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCttt         >  2:1338070/1‑150 (MQ=255)
                                                                                                                                         gCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGATGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCAGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCttt         >  2:1677971/1‑150 (MQ=255)
                                                                                                                                           cAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCTTTAc       <  1:1994433/150‑1 (MQ=255)
                                                                                                                                              cgcCCTGGCGGTGAAAGCTGCGCGGCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCGCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCTTTACCGa    <  2:2402090/150‑1 (MQ=255)
                                                                                                                                                cccTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCTTTACCGAtg  <  1:1269494/150‑1 (MQ=255)
                                                                                                                                                cccTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCTTTACCGAtg  <  2:1553205/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGATGACCTTCAACGAAATTAACAACCAGCGTAACTGGCGTGCACCGCTGTTCGGTTACTGCTGCTCCGGCGTGGTGTATACCGAGCATGAAAACCCGGAAGAGACGATGTATCAGGTGCTGCATCACCAGTTTGTCGCCAGCGCCCTGGCGGTGAAAGCTGCGCGTCGCATTAACCCGGAGATGAAAGTCGGCTGTATGCTGGCGATGGTGCCGCTCTATCCTTACTCCTGTAACCCGGACGATGTGATGTTCGCTCAGGAGTCGATGCGCGAACGCTACGTCTTTACCGATG  >  NC_000913/3044185‑3044478

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