Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,732,549 G→T 9.0% D565Y (GAT→TAT)  malS → alpha‑amylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,732,5490GT 9.0% 711.0 / 27.9 199D565Y (GAT→TAT) malSalpha‑amylase
Reads supporting (aligned to +/- strand):  ref base G (78/103);  new base T (10/8);  total (88/111)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 205 total aligned reads displayed.

CGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGA  >  CP009273/3732402‑3732694
                                                                                                                                                   |                                                                                                                                                 
cGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATg                                                                                                                                                   <  2:293466/148‑1 (MQ=255)
cGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATg                                                                                                                                                   <  2:292869/148‑1 (MQ=255)
cGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATg                                                                                                                                                   <  1:1718288/148‑1 (MQ=255)
 gCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGa                                                                                                                                                  <  1:2510283/148‑1 (MQ=255)
 gCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGa                                                                                                                                                  <  1:496977/148‑1 (MQ=255)
   cACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATa                                                                                                                                                <  2:1362572/148‑1 (MQ=255)
    aCGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAc                                                                                                                                               >  2:2218175/1‑148 (MQ=255)
     cGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATAcc                                                                                                                                              >  1:2556739/1‑148 (MQ=255)
       gCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCg                                                                                                                                            >  1:940465/1‑148 (MQ=255)
        cTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGc                                                                                                                                           >  2:598622/1‑148 (MQ=255)
        cTTCGATGCGATGATCAATGTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGc                                                                                                                                           <  2:244585/148‑1 (MQ=255)
           cgatgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTg                                                                                                                                        <  2:477545/148‑1 (MQ=255)
             atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                      <  2:1932890/148‑1 (MQ=255)
             atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                      >  2:183771/1‑148 (MQ=255)
             atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                      <  2:1800605/148‑1 (MQ=255)
             atgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGtt                                                                                                                                      >  1:1541623/1‑148 (MQ=255)
              tgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTc                                                                                                                                     >  1:2650488/1‑148 (MQ=255)
              tgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTc                                                                                                                                     >  2:3309648/1‑148 (MQ=255)
              tgcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTc                                                                                                                                     <  2:2669502/148‑1 (MQ=255)
               gcgatgATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTcc                                                                                                                                    <  2:1351778/148‑1 (MQ=255)
                    gatCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGaa                                                                                                                               >  1:304426/1‑148 (MQ=255)
                     atCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAg                                                                                                                              <  1:1709740/148‑1 (MQ=255)
                      tCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgg                                                                                                                             <  1:223348/148‑1 (MQ=255)
                      tCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgg                                                                                                                             <  2:1787615/148‑1 (MQ=255)
                      tCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAgg                                                                                                                             <  2:1036208/148‑1 (MQ=255)
                       cAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAggg                                                                                                                            >  1:2746274/1‑148 (MQ=255)
                           ttCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCg                                                                                                                        >  1:841054/1‑148 (MQ=255)
                              gATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACa                                                                                                                     <  2:76417/148‑1 (MQ=255)
                                ttATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACaaa                                                                                                                   >  2:206102/1‑148 (MQ=255)
                                 tATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAag                                                                                                                  >  1:1599381/1‑148 (MQ=255)
                                  aTCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagc                                                                                                                 <  1:565102/148‑1 (MQ=255)
                                   tCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAagca                                                                                                                <  1:1809217/148‑1 (MQ=255)
                                      ggAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAgcagca                                                                                                             <  1:2218175/148‑1 (MQ=255)
                                       gAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGcagcag                                                                                                            <  2:7234/148‑1 (MQ=255)
                                        aGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCaga                                                                                                           >  1:1440842/1‑148 (MQ=255)
                                         gCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCagag                                                                                                          <  1:316998/148‑1 (MQ=255)
                                           aGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtt                                                                                                        <  1:97187/148‑1 (MQ=255)
                                           aGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGtt                                                                                                        <  1:1660648/148‑1 (MQ=255)
                                                 cgAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACt                                                                                                  >  1:2894435/1‑148 (MQ=255)
                                                 cgAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACt                                                                                                  <  1:1863135/148‑1 (MQ=255)
                                                   aaaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAt                                                                                                >  1:2996322/1‑148 (MQ=255)
                                                    aaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                               >  2:2406958/1‑148 (MQ=255)
                                                    aaGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTAtt                                                                                               >  2:358194/1‑148 (MQ=255)
                                                     aGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTa                                                                                              <  1:1257016/148‑1 (MQ=255)
                                                      gCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGGGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAg                                                                                             >  2:2562351/1‑148 (MQ=255)
                                                         gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgcgc                                                                                          >  1:875846/1‑148 (MQ=255)
                                                         gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgcgc                                                                                          <  1:3540471/148‑1 (MQ=255)
                                                         gTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAgcgc                                                                                          <  1:297785/148‑1 (MQ=255)
                                                               tGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGgcg                                                                                    <  1:1488027/148‑1 (MQ=255)
                                                                  cTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCgg                                                                                 <  2:1211699/148‑1 (MQ=255)
                                                                        cAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACaaa                                                                           <  1:2622831/148‑1 (MQ=255)
                                                                         aGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAAt                                                                          <  1:1735596/148‑1 (MQ=255)
                                                                              gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                     >  1:728090/1‑148 (MQ=255)
                                                                              gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                     >  1:3342925/1‑148 (MQ=255)
                                                                              gATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCtttt                                                                     <  2:3230029/148‑1 (MQ=255)
                                                                               aTACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTa                                                                    >  1:1018383/1‑148 (MQ=255)
                                                                               aTACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTa                                                                    >  1:2007243/1‑148 (MQ=255)
                                                                                tACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTAt                                                                   >  2:413450/1‑148 (MQ=255)
                                                                                 aCGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATg                                                                  <  1:2065264/148‑1 (MQ=255)
                                                                                         gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                          >  2:3262258/1‑148 (MQ=255)
                                                                                         gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGaa                                                          >  1:1226254/1‑148 (MQ=255)
                                                                                         gcagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTGTTATGGTGATGaa                                                          >  2:1842976/1‑148 (MQ=255)
                                                                                          cagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAAt                                                         <  1:2967998/148‑1 (MQ=255)
                                                                                          cagcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAAt                                                         <  2:610965/148‑1 (MQ=255)
                                                                                           agcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATc                                                        <  1:2786030/148‑1 (MQ=255)
                                                                                           agcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATc                                                        <  2:256797/148‑1 (MQ=255)
                                                                                           agcaAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATc                                                        <  2:2556739/148‑1 (MQ=255)
                                                                                             caAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCt                                                      >  2:1385700/1‑148 (MQ=255)
                                                                                                  ggCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGt                                                 <  2:1165030/148‑1 (MQ=255)
                                                                                                     ggAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCg                                              <  1:605172/148‑1 (MQ=255)
                                                                                                     ggAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCg                                              <  2:3477548/148‑1 (MQ=255)
                                                                                                      gagaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGt                                             >  2:460327/1‑148 (MQ=255)
                                                                                                       agaAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGtt                                            <  1:1691139/148‑1 (MQ=255)
                                                                                                          aaTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCgg                                         >  1:2677457/1‑148 (MQ=255)
                                                                                                           aTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGt                                        >  1:1941113/1‑148 (MQ=255)
                                                                                                           aTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGt                                        <  2:2130555/148‑1 (MQ=255)
                                                                                                           aTTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGt                                        <  2:2129905/148‑1 (MQ=255)
                                                                                                            ttGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTc                                       <  2:1458658/148‑1 (MQ=255)
                                                                                                               cAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTa                                    <  2:826717/148‑1 (MQ=255)
                                                                                                                     ttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                              >  1:3154763/1‑148 (MQ=255)
                                                                                                                     ttCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGtt                              <  2:1599381/148‑1 (MQ=255)
                                                                                                                       cAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGGTCGGTCCTACAGGTTCt                            >  2:2017041/1‑148 (MQ=255)
                                                                                                                        aaCGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTg                           >  1:1524798/1‑148 (MQ=255)
                                                                                                                          cGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGAt                         <  2:1103130/148‑1 (MQ=255)
                                                                                                                           gtgtTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATc                        <  2:2405222/148‑1 (MQ=255)
                                                                                                                              ttGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGc                     >  1:109831/1‑148 (MQ=255)
                                                                                                                               tGAGCTACCTCTCGTCGCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCt                    <  2:875846/148‑1 (MQ=255)
                                                                                                                                gAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTg                   <  1:2403550/148‑1 (MQ=255)
                                                                                                                                  gCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCa                 >  2:2998925/1‑148 (MQ=255)
                                                                                                                                  gCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCa                 >  2:745506/1‑148 (MQ=255)
                                                                                                                                  gCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCa                 >  1:1998406/1‑148 (MQ=255)
                                                                                                                                   cTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTCCACATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCaa                >  2:1253128/1‑148 (MQ=255)
                                                                                                                                     aCCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAgg              <  2:2773138/148‑1 (MQ=255)
                                                                                                                                      ccTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGt             <  1:769069/148‑1 (MQ=255)
                                                                                                                                        tctcGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTac           >  2:3363024/1‑148 (MQ=255)
                                                                                                                                        tctcGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTac           >  1:3374554/1‑148 (MQ=255)
                                                                                                                                           cgtcgCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACg        <  2:3342974/148‑1 (MQ=255)
                                                                                                                                             tcgCATTATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGtt      >  2:2473009/1‑148 (MQ=255)
                                                                                                                                              cgCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTc     <  2:2748964/148‑1 (MQ=255)
                                                                                                                                                 aTGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGa  <  2:1541623/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCCACGGCTTCGATGCGATGATCAATTTCGATTATCAGGAGCAGGCGGCGAAAGCAGTCGACTGTCTGGCGCAGATGGATACGACCTGGCAGCAAATGGCGGAGAAATTGCAGGGTTTCAACGTGTTGAGCTACCTCTCGTCGCATGATACCCGCCTGTTCCGTGAAGGGGGCGACAAAGCAGCAGAGTTATTACTATTAGCGCCAGGCGCGGTACAAATCTTTTATGGTGATGAATCCTCGCGTCCGTTCGGTCCTACAGGTTCTGATCCGCTGCAAGGTACACGTTCGGA  >  CP009273/3732402‑3732694

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: ‑