Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 168,933 A→G intergenic (‑33482/+44764) unknown ← / ← unknown SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑Arg‑ACG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1168,9330AG98.6% 186.7 / NA 69intergenic (‑33482/+44764)unknown/unknownSSU rRNA ## 16S rRNA, small subunit ribosomal RNA/tRNA‑Arg‑ACG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/38);  minor base T (1/0);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCTTTCACCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTACCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTGCTTCCAGCTTG  >  1/168700‑169167
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
gCTTTCACCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAAc                                                                                                                                                                                                                            <  1:49625/250‑1 (MQ=255)
gCTTTCACCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAAc                                                                                                                                                                                                                            >  1:531495/1‑250 (MQ=255)
   ttCACCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCata                                                                                                                                                                                                                        >  1:325901/1‑251 (MQ=255)
   ttCACCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCata                                                                                                                                                                                                                        >  1:343822/1‑251 (MQ=255)
      aCCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTg                                                                                                                                                                                                                     <  1:749367/251‑1 (MQ=255)
                  tcttTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTg                                                                                                                                                                                                          >  1:493390/1‑250 (MQ=255)
                    tttCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTgcgc                                                                                                                                                                                                       <  1:648294/251‑1 (MQ=255)
                           gCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCat                                                                                                                                                                                                                         >  1:658799/1‑226 (MQ=255)
                               gggATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGa                                                                                                                                                                                             >  1:465071/1‑250 (MQ=255)
                                     gattATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGTTGCCGGCGGCAATGACGATGCCTTTTCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGGTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGAGCGCCTAAGCCCCACGtg                                                                                                                                                                                         >  1:819426/1‑248 (MQ=255)
                                     gattATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTg                                                                                                                                                                                                                     <  1:464077/220‑1 (MQ=255)
                                     gattATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGtg                                                                                                                                                                                         >  1:307929/1‑248 (MQ=255)
                                     gattATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGt                                                                                                                                                                                          <  1:387345/247‑1 (MQ=255)
                                      attATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGtgctg                                                                                                                                                                                      <  1:415205/250‑1 (MQ=255)
                                             ccATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACaa                                                                                                                                                                              <  1:187829/251‑1 (MQ=255)
                                                tACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACg                                                                                                                                                                           <  1:346597/251‑1 (MQ=255)
                                                       ttAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGTAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTcc                                                                                                                                                                    <  1:429552/251‑1 (MQ=255)
                                                            tCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTc                                                                                                                                                               <  1:86703/251‑1 (MQ=255)
                                                            tCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTc                                                                                                                                                               <  1:464513/251‑1 (MQ=255)
                                                                 gTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCg                                                                                                                                                          >  1:531111/1‑251 (MQ=255)
                                                                   cGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGc                                                                                                                                                        <  1:302423/251‑1 (MQ=255)
                                                                       aTGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACg                                                                                                                                                                                           >  1:536960/1‑212 (MQ=255)
                                                                       aTGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCCCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGt                                                                                                                                                    >  1:265355/1‑251 (MQ=255)
                                                                         gTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTc                                                                                                                                                                                                                 <  1:735412/188‑1 (MQ=255)
                                                                          tttCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACg                                                                                                                                                 <  1:163803/251‑1 (MQ=255)
                                                                                 aTCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATg                                                                                                                                          >  1:117658/1‑251 (MQ=255)
                                                                                   cACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATggg                                                                                                                                        <  1:618957/251‑1 (MQ=255)
                                                                                   cACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATggg                                                                                                                                        <  1:167430/251‑1 (MQ=255)
                                                                                   cACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATggg                                                                                                                                        <  1:125595/251‑1 (MQ=255)
                                                                                   cACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAAAGACGATGCCTTTTCCGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGCTGACAACTTTTCCCGCCACCGGGAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGCTTTTCCTGCTCCGTTGGCTTGGACGGTTTATggg                                                                                                                                        >  1:114417/1‑251 (MQ=255)
                                                                                     ccATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTTCGGGTTATggggt                                                                                                                                      >  1:50409/1‑249 (MQ=255)
                                                                                         cAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACg                                                                                                                                                                           >  1:372785/1‑210 (MQ=255)
                                                                                          aaGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGt                                                                                                                                 <  1:178491/251‑1 (MQ=255)
                                                                                            gCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAAc                                                                                                                                                                                <  1:899333/202‑1 (MQ=255)
                                                                                            gCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTcc                                                                                                                                                                    <  1:522967/214‑1 (MQ=255)
                                                                                            gCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTTCTTCTCCGTCGGCTTGTACGGGTTATGGGGTTTGGt                                                                                                                                 >  1:35287/1‑249 (MQ=255)
                                                                                             ccATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGcc                                                                                                                                                                                                                                      <  1:183539/147‑1 (MQ=255)
                                                                                             ccATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACaa                                                                                                                                                                              >  1:263390/1‑203 (MQ=255)
                                                                                                  gCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGtt                                                                                                                                                                         <  1:75151/203‑1 (MQ=255)
                                                                                                           ccAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCgttgtt                                                                                                                <  1:56492/251‑1 (MQ=255)
                                                                                                           ccAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCgttgtt                                                                                                                <  1:54915/251‑1 (MQ=255)
                                                                                                                 ggCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTg                                                                                                                                                                                                                     <  1:563337/144‑1 (MQ=255)
                                                                                                                                cccgTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTgcgc                                                                                                         <  1:751048/237‑1 (MQ=255)
                                                                                                                                   gTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCa                                                                                         >  1:14459/1‑250 (MQ=255)
                                                                                                                                     cGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGt                                                                                      <  1:796151/251‑1 (MQ=255)
                                                                                                                                          aTGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGATCCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAATCCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAGCGCTTCGCGCTCACGTTGctc                                                                                 <  1:620564/251‑1 (MQ=255)
                                                                                                                                            ggCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCtt                                                                               >  1:138938/1‑251 (MQ=255)
                                                                                                                                                tttgccGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCttgttg                                                                           <  1:332454/247‑1 (MQ=255)
                                                                                                                                                                tGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTcc                                                           <  1:151016/251‑1 (MQ=255)
                                                                                                                                                                       ggcAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGtt                                                    >  1:480683/1‑251 (MQ=255)
                                                                                                                                                                        gcAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTg                                                   >  1:163924/1‑251 (MQ=255)
                                                                                                                                                                                 aTGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCt                                           <  1:791421/250‑1 (MQ=255)
                                                                                                                                                                                  tGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTctcctc                                                                                                                           <  1:319253/169‑1 (MQ=255)
                                                                                                                                                                                       ttCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGtc                                                                                                                                >  1:546551/1‑159 (MQ=255)
                                                                                                                                                                                       ttCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCaa                                    >  1:640351/1‑251 (MQ=255)
                                                                                                                                                                                       ttCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCaa                                    >  1:245058/1‑251 (MQ=255)
                                                                                                                                                                                               gAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGt                                                                                      <  1:787345/193‑1 (MQ=255)
                                                                                                                                                                                                    cGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGctc                                                                                 >  1:709587/1‑193 (MQ=255)
                                                                                                                                                                                                     gCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGtttt                                                                                                              >  1:894309/1‑163 (MQ=255)
                                                                                                                                                                                                      cccATTTGGGTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCt                                                                                                                                                       >  1:470523/1‑121 (MQ=255)
                                                                                                                                                                                                        cATTTGGGTCGAACCAAAGCGGTTGACACCTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCtt                   <  1:208807/251‑1 (MQ=255)
                                                                                                                                                                                                            tGGGTCGACCCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGTCGTGCTGGTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTTCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAg               <  1:508572/251‑1 (MQ=255)
                                                                                                                                                                                                              ggTCGAACCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTg             >  1:151487/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaCCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCGCGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAg               >  1:451845/1‑244 (MQ=255)
                                                                                                                                                                                                                   aaCCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTGCTTc         >  1:82331/1‑250 (MQ=255)
                                                                                                                                                                                                                   aaCCAAAGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTGCTTc         >  1:825963/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaaGCGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTGCTTCCAgct    <  1:103093/251‑1 (MQ=255)
                                                                                                                                                                                                                           cGGTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTGCTTCCAgcttg  <  1:59066/249‑1 (MQ=255)
                                                                                                                                                                                                                             gTTGACAACTTTGCCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCt                                                                                                                                                       <  1:282737/98‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
GCTTTCACCAACTGCCCTTCTTTCACCGCCAGGGATTGATTATACCCATACAAACTTAAATCGCTGTCGCCATGTTTCACAATCACCATCAAGCCATAGCCCTGCAACCAGTTGGCTAAAATCACCCGCCCGTCGGCAATGGCTTTTACCGGCGTGCCGCTGCCGGCGGCAATGACGATGCCTTTCCAGCGGAGTTCGCCCATTTGGGTCGAACCAAAGCGGTTGACAACTTTACCCGCCACCGGAAAACCATATTGTTTCTGCGGTGCGCCTAAGCCCGACGTGCTGCTTAACAAACGTTTTTCCTGCTCCGTCGGCTTGTACGGTTTATGGGTTTTGGTCTCCTCCGCCCGTTGTTTTTGCGCCAACGCTTCGCGCTCACGTTGCTCTTGTTGTTTCACCGCCTGTTCCGCCCGTTGAATTTCCTGACGCAAAGCGTTTTCATTGGCTTTCAGTGCTTCCAGCTTG  >  1/168700‑169167

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