Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,987,235 G→A A26V (GCG→GTG)  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,987,2350GA96.5% 152.4 / NA 57A26V (GCG→GTG) unknown16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/31);  minor base C (0/2);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGATGTTGTCTGTTGGCATAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGACGTTGTAAATCGGCGAGAAAATCATCGTTCGCCGCCCCGCTGAGCAGCAGGGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAAAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGCAATTGG  >  1/1986993‑1987454
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
ccaGTAAAACTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCAcccc                                                                                                                                                                                                                         >  1:748261/1‑247 (MQ=255)
 cagtaaaaCTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGa                                                                                                                                                                                                                    >  1:70827/5‑251 (MQ=255)
   gtaaaaCTTCCTGTCCGGCGGTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGa                                                                                                                                                                                                                    >  1:69385/3‑249 (MQ=255)
                    gcggTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCg                                                                                                                                                                                                   >  1:374102/1‑249 (MQ=255)
                    gcggTACATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCg                                                                                                                                                                                                   >  1:621028/1‑249 (MQ=255)
                            ttGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGa                                                                                                                                                                                                                    <  1:53771/224‑1 (MQ=255)
                                gagagaAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCaa                                                                                                                                                                                     >  1:161844/1‑251 (MQ=255)
                                    gaAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCgg                                                                                                                                                                                                                       <  1:500441/213‑2 (MQ=255)
                                         tGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTc                                                                                                                                                                            <  1:111770/251‑1 (MQ=255)
                                                 cTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAacac                                                                                                                                                                      <  1:696085/249‑1 (MQ=255)
                                                 cTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACtt                                                                                                                                                                    <  1:833122/251‑1 (MQ=255)
                                                   ggAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCg                                                                                                                                                                  >  1:227864/1‑251 (MQ=255)
                                                                cTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACaaa                                                                                                                                                     <  1:798096/251‑1 (MQ=255)
                                                                 tGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAAt                                                                                                                                                    >  1:22942/1‑251 (MQ=255)
                                                                 tGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAAt                                                                                                                                                    >  1:161682/1‑251 (MQ=255)
                                                                  gTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGcagcagag                                                                                                                                                                                                           <  1:17806/195‑3 (MQ=255)
                                                                              tCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:657975/251‑1 (MQ=255)
                                                                              tCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:288961/251‑1 (MQ=255)
                                                                              tCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:380476/251‑1 (MQ=255)
                                                                                cAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCt                                                                                                                                     <  1:284892/251‑1 (MQ=255)
                                                                                cAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCt                                                                                                                                     <  1:117471/251‑1 (MQ=255)
                                                                                    aaaaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGc                                                                                                                                                                                 <  1:790000/203‑1 (MQ=255)
                                                                                    aaaaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATgg                                                                                                                                 <  1:842067/251‑1 (MQ=255)
                                                                                     aaaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGc                                                                                                                                <  1:268589/251‑1 (MQ=255)
                                                                                      aaCAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCt                                                                                                                               <  1:282305/251‑1 (MQ=255)
                                                                                           gATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGAc                                                                                                                                                                                                         <  1:753793/172‑1 (MQ=255)
                                                                                            aTCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCgg                                                                                                                                                                                                                       >  1:102828/1‑156 (MQ=255)
                                                                                                          gCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAACAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGc                                                                                                                                >  1:442945/1‑230 (MQ=255)
                                                                                                             aTGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCAt                                                                                                        <  1:205206/251‑1 (MQ=255)
                                                                                                                aaaTTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATtat                                                                                                     <  1:325650/251‑1 (MQ=255)
                                                                                                                  aTTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   <  1:560422/251‑1 (MQ=255)
                                                                                                                  aTTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   <  1:594207/251‑1 (MQ=255)
                                                                                                                                 ccACCGCACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACaaa                                                                                                                                                     >  1:481550/1‑186 (MQ=255)
                                                                                                                                       cACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAac                                                                                                                                                                        >  1:462963/1‑161 (MQ=255)
                                                                                                                                       cACTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCt                                                                                                                                     <  1:687501/196‑1 (MQ=255)
                                                                                                                                         cTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGc                                                                                                                                                                                 <  1:742007/150‑1 (MQ=255)
                                                                                                                                         cTTTGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGc                                                                                                                                >  1:382188/1‑199 (MQ=255)
                                                                                                                                           ttGTGGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       <  1:393692/190‑1 (MQ=255)
                                                                                                                                             gtgGCGTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTg                                                                                                 >  1:549478/1‑226 (MQ=255)
                                                                                                                                                  gTTTGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   >  1:656655/1‑219 (MQ=255)
                                                                                                                                                   tttGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAAc                                                                   <  1:673120/250‑1 (MQ=255)
                                                                                                                                                   tttGGTGTTGTCGGTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAAc                                                                   <  1:503310/250‑1 (MQ=255)
                                                                                                                                                              cggTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTgg                                                                                                                         >  1:739912/3‑185 (MQ=255)
                                                                                                                                                              cggTTGGCGTAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGacacac                                                       <  1:564288/249‑1 (MQ=255)
                                                                                                                                                                   ggcgtAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTcc                                                   >  1:10092/5‑250 (MQ=255)
                                                                                                                                                                         gTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGAtt                                                                                    >  1:598067/1‑211 (MQ=255)
                                                                                                                                                                           cAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATcc                                                                                                                                      <  1:753925/159‑1 (MQ=255)
                                                                                                                                                                                     gCTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       >  1:666166/1‑148 (MQ=255)
                                                                                                                                                                                      cTCCATGCGCGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATc                                                                                                                                       >  1:717476/1‑147 (MQ=255)
                                                                                                                                                                                            gcgcGGGCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATAc                                                                                                   <  1:583273/177‑1 (MQ=255)
                                                                                                                                                                                                 ggCGTTGAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCa                    >  1:450767/1‑251 (MQ=255)
                                                                                                                                                                                                       gAATCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGc               >  1:536005/1‑250 (MQ=255)
                                                                                                                                                                                                         aTCCCTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCaa            >  1:185177/1‑251 (MQ=255)
                                                                                                                                                                                                           cccTGGCGTTGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGata                                                                             <  1:133275/184‑1 (MQ=255)
                                                                                                                                                                                                             ctggcGTTGTAAATCGGCGAGAAAATCATCATTCACCCCCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGc        <  1:262236/247‑1 (MQ=255)
                                                                                                                                                                                                             ctggcGTTGTAAATCGGCGAGAAAATCATCATTCACCCCCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGc        <  1:81522/247‑1 (MQ=255)
                                                                                                                                                                                                                   ctGTAAATCGGCGAGAAAATCATCATTCACCACCCCGGTGAGCAGCAGAGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAGAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGCAATTgg  <  1:894461/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
CCAGCAAAACCTCCTGCCCGGCGGTGCATTGTGAGAGAAGCTGCATGAGCTGGAAATGGCATTCCTGTTTATTCTTCGTCCAATAAAACAGGATCACATCGGCGCTGCCATGAAATTCAAGGGAAAAATCCACCGCACTTTGTGGCGTTTGATGTTGTCTGTTGGCATAGTCAAAATACCAGCTCCATGCGCGGGCGTTGAATCCCTGACGTTGTAAATCGGCGAGAAAATCATCGTTCGCCGCCCCGCTGAGCAGCAGGGACTTCCCGCTTAAAAAATCCAAATGCCGTTCCAACACTTCGCTTTCTGCACAAATCACGCTATTAATCCTATGGCTTATTGGTTAAAATGGGCGCCCATTATACTGGATTGAATTGATTTTAGATATGCAAAAAACGCGCCGACACACTTCCCGCCGAGATTTATTATTGCAGCAAATGGGCATCAGCCAATGGCAATTGG  >  1/1986993‑1987454

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