Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,986,831 G→A P161S (CCA→TCG)  unknown ← 16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,986,8310GA98.6% 209.1 / ‑4.8 72P161S (CCA→TCG) unknown16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (47/24);  total (47/25)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCATGAATATCTGTCATGATTAATCGGGTGTTCGGGTTGTGCCGTTGAATATAACTGCCGATCACTCCCGCGCCGCAGCCCACGTCCAACACCCCGCCTGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGTTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAATTTGTGGCGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAACTGAAAATGATACAGACTGCAACGGCGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACGCCGCTACGATTTTCGCCCACCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTC  >  1/1986588‑1987066
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
gATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGa                                                                                                                                                                                                                                      >  1:392892/1‑251 (MQ=255)
 aTTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGAt                                                                                                                                                                                                                                     <  1:646635/251‑1 (MQ=255)
     tCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGa                                                                                                                                                                                                                                 <  1:555053/251‑1 (MQ=255)
      ccTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGAt                                                                                                                                                                                                                                >  1:544524/1‑251 (MQ=255)
           aaaGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCa                                                                                                                                                                                                                           >  1:252635/1‑251 (MQ=255)
                     ttGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGtttt                                                                                                                                                                                                                  >  1:857387/1‑250 (MQ=255)
                                           gtGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGAtttt                                                                                                                                                                                           >  1:174457/2‑251 (MQ=255)
                                              aaTATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTa                                                                                                                                                                                        >  1:152147/1‑251 (MQ=255)
                                                       cATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAAt                                                                                                                                                                               <  1:618059/251‑1 (MQ=255)
                                                              aaTCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAg                                                                                                                                                                        <  1:730085/251‑1 (MQ=255)
                                                                      gtTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCaa                                                                                                                                                                 >  1:732606/1‑250 (MQ=255)
                                                                      gtTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCaa                                                                                                                                                                 >  1:850159/1‑250 (MQ=255)
                                                                       ttCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAAc                                                                                                                                                                <  1:257430/250‑1 (MQ=255)
                                                                        tCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAAcga                                                                                                                                                              <  1:792504/251‑2 (MQ=255)
                                                                           ggTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGTTGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATAc                                                                                                                                                                          >  1:304567/1‑236 (MQ=255)
                                                                                    gTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAAcgacg                                                                                                                                                            >  1:129192/1‑238 (MQ=255)
                                                                                    gTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCg                                                                                                                                                  >  1:151337/1‑251 (MQ=255)
                                                                                    gTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGAGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCg                                                                                                                                                  >  1:510831/1‑251 (MQ=255)
                                                                                       gAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGAtt                                                                                                                                               >  1:180815/1‑251 (MQ=255)
                                                                                       gAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGAtt                                                                                                                                               >  1:658031/1‑251 (MQ=255)
                                                                                         atatAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCGGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGAtttt                                                                                                                                             >  1:615956/1‑251 (MQ=255)
                                                                                                 gccaaTCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAatat                                                                                                                                     >  1:791005/5‑251 (MQ=255)
                                                                                                  ccGATCACTCCCGCGCCGCAGCCCACGTCCAGCACCTCGCCTGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGTTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAATTTGTGGCGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAACTGAAAATGATACAGACTGCAACGGCGGGCGGAGTCGATTTTGTGAatat                                                                                                                                     <  1:11417/250‑1 (MQ=255)
                                                                                                    aaTCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCg                                                                                                                                   >  1:393609/2‑250 (MQ=255)
                                                                                                       cacccctgcgcCGCAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGt                                                                                                                               >  1:839009/8‑251 (MQ=255)
                                                                                                                   gcAGCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCa                                                                                                                   >  1:374079/1‑251 (MQ=255)
                                                                                                                      gCCCACGTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAatat                                                                                                                                     >  1:381684/1‑230 (MQ=255)
                                                                                                                            gTCCAACGCCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCAtttttt                                                                                                             <  1:751717/248‑1 (MQ=255)
                                                                                                                            gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGa                                                                                                                                         <  1:726661/220‑1 (MQ=255)
                                                                                                                            gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAata                                                                                                                                      <  1:716407/223‑1 (MQ=255)
                                                                                                                            gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTcgg                                                                                                          >  1:406534/1‑251 (MQ=255)
                                                                                                                            gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTcgg                                                                                                          >  1:874023/1‑251 (MQ=255)
                                                                                                                            gTCCAACACCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTcgg                                                                                                          >  1:364997/1‑251 (MQ=255)
                                                                                                                                  cacCTCGCCCGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGTCGGaa                                                                                                     >  1:276473/1‑250 (MQ=255)
                                                                                                                                      ccgcccgcGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcac                                                                                                <  1:782685/251‑1 (MQ=255)
                                                                                                                                         cccgCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcacacc                                                                                             <  1:685270/248‑4 (MQ=255)
                                                                                                                                         cccgCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcacacc                                                                                             <  1:235148/248‑4 (MQ=255)
                                                                                                                                         cccgCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcacac                                                                                              <  1:831075/247‑3 (MQ=255)
                                                                                                                                                 gTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAAcga                                                                                                                                                              >  1:68569/1‑177 (MQ=255)
                                                                                                                                                 gTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTaa                                                                                                                             >  1:109320/1‑211 (MQ=255)
                                                                                                                                                 gTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGa                                                                                      >  1:645805/1‑250 (MQ=255)
                                                                                                                                                  ttttttAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGAtt                                                                                    <  1:7695/251‑1 (MQ=255)
                                                                                                                                                   tttttAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCAtttttt                                                                                                             <  1:391161/225‑1 (MQ=255)
                                                                                                                                                              gAAAGCAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAAcga                                                                                                                                                              >  1:744681/1‑164 (MQ=255)
                                                                                                                                                                   cAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGTTAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTaaaa                                                                   >  1:327340/1‑251 (MQ=255)
                                                                                                                                                                   cAATAATGCCGTGCCGCCATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTaaaa                                                                   >  1:267767/1‑251 (MQ=255)
                                                                                                                                                                                ccgccATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCcgg                                                       >  1:726689/1‑250 (MQ=255)
                                                                                                                                                                                  gccATCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGAc                                                                                                                                                                      <  1:827607/137‑1 (MQ=255)
                                                                                                                                                                                     aTCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGAc                                                                                                                                                                      <  1:429329/134‑1 (MQ=255)
                                                                                                                                                                                     aTCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGa                                                                                      >  1:366719/1‑214 (MQ=255)
                                                                                                                                                                                     aTCCAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCCCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGta                                                  >  1:768600/1‑249 (MQ=255)
                                                                                                                                                                                       ccAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                               <  1:385221/251‑1 (MQ=255)
                                                                                                                                                                                       ccAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                               <  1:899841/251‑1 (MQ=255)
                                                                                                                                                                                       ccAGCTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                               <  1:452063/251‑1 (MQ=255)
                                                                                                                                                                                                     cTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAg                                    >  1:150425/1‑248 (MQ=255)
                                                                                                                                                                                                       gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAAcga                                                                                                                                                              >  1:632610/1‑123 (MQ=255)
                                                                                                                                                                                                       gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACAt                                               <  1:599727/235‑1 (MQ=255)
                                                                                                                                                                                                       gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCa                               >  1:110471/1‑251 (MQ=255)
                                                                                                                                                                                                       gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCa                               >  1:699636/1‑251 (MQ=255)
                                                                                                                                                                                                       gAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCa                               >  1:194344/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCAACGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTgag                                         <  1:275620/230‑1 (MQ=255)
                                                                                                                                                                                                                     aCGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                 >  1:99561/1‑251 (MQ=255)
                                                                                                                                                                                                                      cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGAc                                                                                                                                                                      <  1:765871/101‑1 (MQ=255)
                                                                                                                                                                                                                      cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGa                                                                                      >  1:352346/1‑181 (MQ=255)
                                                                                                                                                                                                                      cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                 >  1:749738/1‑250 (MQ=255)
                                                                                                                                                                                                                      cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                 >  1:845888/1‑250 (MQ=255)
                                                                                                                                                                                                                      cGCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATgg                 >  1:212149/1‑250 (MQ=255)
                                                                                                                                                                                                                       gCATAGATAGTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCa               >  1:169108/1‑251 (MQ=255)
                                                                                                                                                                                                                             atagTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAacgcgcaccactacg                                                                                       >  1:217985/1‑162 (MQ=255)
                                                                                                                                                                                                                              tagTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTc            >  1:783633/1‑247 (MQ=255)
                                                                                                                                                                                                                              tagTTAAATCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTc            >  1:529807/1‑247 (MQ=255)
                                                                                                                                                                                                                                    aaTCGGCAAGTTGCGACGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAATTGAAAATGATACAGACTGCAACGACGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACACCACTACGATTTTCGCCCACCAGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTAttc  <  1:244789/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
GATTTTCCTGTAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCATGAATATCTGTCATGATTAATCGGGTGTTCGGGTTGTGCCGTTGAATATAACTGCCGATCACTCCCGCGCCGCAGCCCACGTCCAACACCCCGCCTGCGATGTTTTTTAGCGTAGAAAGCAATAATGCCGTGCCGCCATCCAGTTCGGCGGCGCTGAATACGCCGGGCAACGCATAGATAGTTAAATCGGCAATTTGTGGCGATTGATAGGATTTCCAGAAAGTTTTCAGTTCAAAGTGCGGTGGATTTTTTAACTGAAAATGATACAGACTGCAACGGCGGGCGGAGTCGATTTTGTGAATATCGCCGTAATCGGCAAGCATTTTTTCGGCGGAACGCACGCCGCTACGATTTTCGCCCACCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTC  >  1/1986588‑1987066

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