Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 484,248 T→C intergenic (+120831/‑87771) unknown → / → unknown tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1484,2480TC98.7% 228.7 / NA 78intergenic (+120831/‑87771)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (42/35);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CAAAGGCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACTCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCAAATCCGCTAACCAAGCAG  >  1/484009‑484491
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
caAAGGCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgc                                                                                                                                                                                                                                          <  1:138856/251‑1 (MQ=255)
    ggCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTg                                                                                                                                                                                                                                              <  1:374253/243‑1 (MQ=255)
    ggCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgcagca                                                                                                                                                                                                                                      >  1:866608/1‑251 (MQ=255)
    ggCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgcagca                                                                                                                                                                                                                                      >  1:257792/1‑251 (MQ=255)
    ggCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGCGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgcagca                                                                                                                                                                                                                                      >  1:33567/1‑251 (MQ=255)
           aCCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTg                                                                                                                                                                                                                                <  1:370188/250‑1 (MQ=255)
              aTAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGAtgc                                                                                                                                                                                                                            <  1:812308/251‑1 (MQ=255)
                 gCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTg                                                                                                                                                                                                                                >  1:879390/1‑244 (MQ=255)
                  cTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGTGCAGCACACCTGATGCTGCa                                                                                                                                                                                                                        >  1:660051/1‑251 (MQ=255)
                      tattTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTc                                                                                                                                                                                                                    <  1:817265/251‑1 (MQ=255)
                                                aCATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACAt                                                                                                                                                                                          >  1:329047/1‑251 (MQ=255)
                                                  atatTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGt                                                                                                                                                                                        >  1:141720/1‑251 (MQ=255)
                                                  atatTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGt                                                                                                                                                                                        >  1:233660/1‑251 (MQ=255)
                                                      tGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGAc                                                                                                                                                                                    <  1:761233/251‑1 (MQ=255)
                                                                     gtgtTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTg                                                                                                                                                                                                                                <  1:834512/192‑1 (MQ=255)
                                                                     gtgtTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTaa                                                                                                                                                                      >  1:169807/1‑250 (MQ=255)
                                                                       gtTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAg                                                                                                                                                                                                      >  1:196449/1‑216 (MQ=255)
                                                                                    atCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTc                                                                                                                                                      <  1:185030/251‑1 (MQ=255)
                                                                                      cGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTccc                                                                                                                                                    <  1:158171/251‑1 (MQ=255)
                                                                                        gTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGAtg                                                                                                                                                                                                                             <  1:691081/176‑1 (MQ=255)
                                                                                        gTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTcc                                                                                                                                                     >  1:570224/1‑248 (MQ=255)
                                                                                            tAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCg                                                                                                                                              <  1:819657/251‑1 (MQ=255)
                                                                                               cGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTAc                                                                                                                                           <  1:526394/251‑1 (MQ=255)
                                                                                               cGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTAc                                                                                                                                           <  1:22334/251‑1 (MQ=255)
                                                                                                 aGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACt                                                                                                                                          <  1:633698/250‑1 (MQ=255)
                                                                                                 aGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTACATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACTCTGATTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACAAGTTGACGGCACGAAATTTAACAATCCGGAAGTACTc                                                                                                                                                      >  1:550149/1‑238 (MQ=255)
                                                                                                          tttCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAaag                                                                                                                                <  1:117390/251‑1 (MQ=255)
                                                                                                              cATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTccc                                                                                                                                                    >  1:687004/1‑227 (MQ=255)
                                                                                                                             aaaTGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTGGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGc                                                                                                              <  1:430630/250‑1 (MQ=255)
                                                                                                                                tGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTc                                                                                                            >  1:390201/1‑249 (MQ=255)
                                                                                                                                  aaGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGt                                                                                                         <  1:537387/250‑1 (MQ=255)
                                                                                                                                       tCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCg                                                                                                                                              >  1:265986/1‑208 (MQ=255)
                                                                                                                                       tCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAaagaag                                                                                                                             >  1:747116/1‑225 (MQ=255)
                                                                                                                                       tCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACa                                                                                                   >  1:577070/1‑251 (MQ=255)
                                                                                                                                       tCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACa                                                                                                   >  1:24937/1‑251 (MQ=255)
                                                                                                                                        ccATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAAc                                                                                                                                                                                            <  1:413058/161‑1 (MQ=255)
                                                                                                                                        ccATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACAt                                                                                                  <  1:187096/251‑1 (MQ=255)
                                                                                                                                        ccATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACAt                                                                                                  >  1:553675/1‑251 (MQ=255)
                                                                                                                                        ccATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACAt                                                                                                  >  1:333768/1‑251 (MQ=255)
                                                                                                                                             gCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgc                                                                                                                                                                                                                                          <  1:28078/110‑1 (MQ=255)
                                                                                                                                             gCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAAc                                                                                                                                                                     >  1:321585/1‑179 (MQ=255)
                                                                                                                                             gCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGaaa                                                                                             >  1:420224/1‑251 (MQ=255)
                                                                                                                                              ccccTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGaa                                                                                              <  1:401612/249‑1 (MQ=255)
                                                                                                                                              ccccTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGaa                                                                                              <  1:460905/249‑1 (MQ=255)
                                                                                                                                                 cTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGa                                                                                         >  1:435206/1‑251 (MQ=255)
                                                                                                                                                   attattTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGAc                                                                                                                                                                                    <  1:573201/158‑1 (MQ=255)
                                                                                                                                                      attTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCaa                                                                                     >  1:738706/1‑250 (MQ=255)
                                                                                                                                                       tttAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAAc                                                                                    >  1:530584/1‑250 (MQ=255)
                                                                                                                                                             gAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCca                                                                                                                                      <  1:191342/194‑1 (MQ=255)
                                                                                                                                                             gAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTcc                                                                              >  1:361442/1‑250 (MQ=255)
                                                                                                                                                              aaTTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCcg                                                                            <  1:531776/251‑1 (MQ=255)
                                                                                                                                                               aTTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgc                                                                                                                                                                                                                                          <  1:376098/92‑1 (MQ=255)
                                                                                                                                                               aTTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCcg                                                                            >  1:132228/1‑250 (MQ=255)
                                                                                                                                                               aTTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCcg                                                                            >  1:569980/1‑250 (MQ=255)
                                                                                                                                                                 tttAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTCCAAAACATGTTGACGGCACGAAATTTAAc                                                                                                                                                                     >  1:631061/1‑159 (MQ=255)
                                                                                                                                                                        aGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTc                                                                                                                                                                                                                    <  1:389529/105‑1 (MQ=255)
                                                                                                                                                                         gTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGgc                                                                                                                                                                                                                                          >  1:837381/1‑82 (MQ=255)
                                                                                                                                                                              ttCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCt                                                            <  1:776586/251‑1 (MQ=255)
                                                                                                                                                                                cATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAg                                                          >  1:585735/1‑251 (MQ=255)
                                                                                                                                                                                cATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAg                                                          >  1:227138/1‑251 (MQ=255)
                                                                                                                                                                                                tGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGaaa                                          <  1:510300/251‑1 (MQ=255)
                                                                                                                                                                                                tGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGaaa                                          <  1:1963/251‑1 (MQ=255)
                                                                                                                                                                                                   cTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGc                                                                                                              >  1:226540/1‑180 (MQ=255)
                                                                                                                                                                                                   cTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGaaaaa                                        >  1:470592/1‑250 (MQ=255)
                                                                                                                                                                                                    tCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCGAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAAcc                                      <  1:429213/251‑1 (MQ=255)
                                                                                                                                                                                                          ttttGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGAcc                                                             <  1:741626/222‑1 (MQ=255)
                                                                                                                                                                                                          ttttGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTg                                >  1:420562/1‑251 (MQ=255)
                                                                                                                                                                                                                  cAGCGGGCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGAtt                         >  1:827604/1‑250 (MQ=255)
                                                                                                                                                                                                                       ggCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCaa                   >  1:246129/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggCTAGGTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCaa                   >  1:321061/1‑251 (MQ=255)
                                                                                                                                                                                                                            ggTGTTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAAc                                                                                    >  1:840823/1‑181 (MQ=255)
                                                                                                                                                                                                                             gtgtTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGAtgct                                                                                                                                                                                                                           >  1:625105/1‑45 (MQ=255)
                                                                                                                                                                                                                               gtTTTTAAACTATCACCCTGTTTGGGGCAGCACGCCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGt                                                                                                         <  1:646733/157‑1 (MQ=255)
                                                                                                                                                                                                                               gtTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCAAATCCGCTa           >  1:808289/1‑251 (MQ=255)
                                                                                                                                                                                                                               gtTTTTAAACTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCAAATCCGCTa           >  1:16380/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACAt                                                                                                  <  1:668807/155‑1 (MQ=255)
                                                                                                                                                                                                                                        cTATCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCAAATCCGCTAACCAAGCAg  >  1:172826/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTCACCCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGAt                                  <  1:384260/217‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
CAAAGGCTACGACCATAGCTATTATTTCATCGCCAGTTTTATTGGGGAACATATTGCATATCACGCGCAGTGTTTAGCTAATTAATCGGTTCTAGCGAGGCTAAATTTTCCATCAAAGTGCAGTGAAATGAAGCGTCCATAGCCCCTATTATTTAAGGAATTTTAATGAGTCTATTCATGAAAATCATCAAATGACTCATGATTTTGGCGCAGCGGGCTAGGTGTTTTTAAACTATCACTCTGTTTGGGGCAGCACACCTGATGCTGCAAGTCTTGCTAAAATGTAGGTTTCAAAACATGTTGACGGCACGAAATTTAACAATCCGGAAGTACTCCCTGTGCGTACTGCCACAAAAGAAGAGCAAACCTCATTGCTCGGTTGGACATGGAAAATGATCAACCTTTCCCGCGGGTAAATTCGACCTAGACAATCACTCGTGAAAAACCCGATTGTTCGATTCACCAAATCCGCTAACCAAGCAG  >  1/484009‑484491

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