Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,018,917 T→C intergenic (‑12461/‑16302) unknown ← / → unknown tRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,018,9170TC98.6% 174.7 / ‑6.0 73intergenic (‑12461/‑16302)unknown/unknowntRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/41);  minor base G (1/0);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGTCTGGTTTAATAACTTCATAATTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGTGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTT  >  1/1018670‑1019146
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
tAGGTCTGGTTTAATAACTTCATAATTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTggcggt                                                                                                                                                                                                                                          <  1:652412/251‑1 (MQ=255)
       ggTTTAATAACTTCATAATTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGTATTGCCg                                                                                                                                                                                                                                   >  1:561229/1‑251 (MQ=255)
                      tAATTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCg                                                                                                                                                                                                                                   <  1:377091/236‑1 (MQ=255)
                        aTTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGGCGGATCTGCAATCgggggg                                                                                                                                                                                                                  >  1:719506/1‑247 (MQ=255)
                        aTTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTTCAAACggtggg                                                                                                                                                                                                                  >  1:183352/1‑250 (MQ=255)
                        aTTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCggtggt                                                                                                                                                                                                                  <  1:397606/251‑1 (MQ=255)
                          tttAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCggtggtgg                                                                                                                                                                                                                <  1:190934/251‑1 (MQ=255)
                          tttAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCggtggtgg                                                                                                                                                                                                                <  1:117182/251‑1 (MQ=255)
                                   cATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTag                                                                                                                                                                                                       <  1:213436/251‑1 (MQ=255)
                                     ttttATGAAAAAAATCGTATTAGCTGTGTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTagag                                                                                                                                                                                                     <  1:490150/251‑1 (MQ=255)
                                           gAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGGATTGCCGGGTCTGCCATCGGTGGGGGGACAGGGagag                                                                                                                                                                                                     >  1:242313/1‑245 (MQ=255)
                                                            cTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGGATTGCCGGATATGCAATCGGGGGGGGGACCGGGAGAGCGCTTGCAACGGCT‑ggttgggga                                                                                                                                                                             >  1:388156/1‑243 (MQ=255)
                                                              gTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGCGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGGGGGACAGGGAGAGGGCTTGCAACgggcgttgggggaaa                                                                                                                                                                            >  1:19934/1‑238 (MQ=255)
                                                                 ttAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑c                                                                                                                                                                       <  1:560823/251‑1 (MQ=255)
                                                                 ttAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGGGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTgggggaattt                                                                                                                                                                         >  1:54619/1‑242 (MQ=255)
                                                                       ggCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGGGGGACAGGGAGAGGGGTTGCaa                                                                                                                                                                                            >  1:668986/1‑226 (MQ=255)
                                                                           tAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAgcg                                                                                                                                                                                                   <  1:738897/215‑1 (MQ=255)
                                                                                 ttGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGc                                                                                                                                                       <  1:863107/251‑1 (MQ=255)
                                                                                          gtgCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTag                                                                                                                                                                                                       <  1:526272/196‑1 (MQ=255)
                                                                                          gtgCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCCGATCTGCAATTGGTGGGGGTACAGGGAGAGCGCTTGCCACGGCT‑GGTGGGGCAAATG‑CCCGGGCCGTTAACCGCAGTaaaggt                                                                                                                                              >  1:35780/1‑247 (MQ=255)
                                                                                                    aGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGt                                                                                                                                              <  1:228415/240‑1 (MQ=255)
                                                                                                       tGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGGGGGACAGGTAGAGCGCTTGCAACGGCT‑GGTGGAggcaatgccggggc                                                                                                                                                                  >  1:156327/1‑206 (MQ=255)
                                                                                                       tGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGGTCTGCAATCgggggtgggac                                                                                                                                                                                                             >  1:246644/1‑174 (MQ=255)
                                                                                                        gTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGTGGGGGTACAGGGAGAgc                                                                                                                                                                                                    >  1:642017/1‑185 (MQ=255)
                                                                                                           tataGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCggtggtgggac                                                                                                                                                                                                             >  1:545598/1‑170 (MQ=255)
                                                                                                             taGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTTCCGGATCTGCAAACGGTGGGGGGACAggg                                                                                                                                                                                                         >  1:770225/1‑174 (MQ=255)
                                                                                                               gTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGt                                                                                                                                                                                                         <  1:514474/173‑1 (MQ=255)
                                                                                                                     gtgtgtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGGGGTGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTGGGGCaa                                                                                                                                                                            >  1:215070/1‑195 (MQ=255)
                                                                                                                     gtgtgtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGTATTGCCGGGTCTGCAATCGGTGGGGGGACAGGGAGGGCGCTTGCAACGGCT‑GTTGGGGGAATTG‑CCGGGGCAGTTACCGGCCGGAA‑GGTTGAAGAAAAAGTCAGCCCAGTTAgcg‑c                                                                                                                  >  1:752443/1‑250 (MQ=255)
                                                                                                                      tgtgtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAgcg‑c                                                                                                                  <  1:778471/249‑1 (MQ=255)
                                                                                                                      tgtgtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAgcg‑c                                                                                                                  <  1:536949/249‑1 (MQ=255)
                                                                                                                       gtgtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGc                                                                                                                                                                              <  1:695742/191‑1 (MQ=255)
                                                                                                                       gtgtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAg                                                                                                                                                               <  1:492672/205‑1 (MQ=255)
                                                                                                                         gtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGGAGAGCGCTTGCAACGGCT‑GTTGGGGCAATTG‑TCGGTGCTGGTACCGGGAGGAA‑AGGTGAAGAAAAAGACAGCCAAggtttgtctcttt                                                                                                              >  1:303261/1‑239 (MQ=255)
                                                                                                                         gtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTTGGGCAATTGCCCGGGGC‑GGTACCGGCAGGAA‑GGGTGAAGAAAAAGGCAGGCAAGGTAGCG‑CCCTgg                                                                                                             >  1:9602/1‑251 (MQ=255)
                                                                                                                         gtATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGGGGGACAGGGAGAGCGCTTGCAACGGGT‑GGTGGGGCAATTG‑CCGGGGCAGTTACCGGCAGGAA‑GGTTGGAGaaatgtctcttttaaaaatctcccag                                                                                                               >  1:469203/1‑227 (MQ=255)
                                                                                                                            ttcttcGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGaaa                                                                                                          <  1:374867/251‑1 (MQ=255)
                                                                                                                                   ggACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGc                                                                                                                                                                                        <  1:898644/170‑1 (MQ=255)
                                                                                                                                        aGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGt                                                                                                                                                                                                         <  1:455020/148‑1 (MQ=255)
                                                                                                                                        aGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATgg                                                                                                       <  1:876433/242‑1 (MQ=255)
                                                                                                                                           aaaaGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAg                                                                                          <  1:261418/251‑1 (MQ=255)
                                                                                                                                                   cgcgTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGa                                                                                 <  1:233844/251‑1 (MQ=255)
                                                                                                                                                       ttCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGGGGTACAGGt                                                                                                                                                                                                         >  1:41793/1‑133 (MQ=255)
                                                                                                                                                       ttCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGGGGTACAGGGAGAGCGCTTGCAACGGCT‑GTTGGAGCAAATG‑CCGGTGCAGTTACCGGCAGGAAGGGTTG‑AGAAAAAGGCAGCCAAGTTAGCGCCCTGGGAAATGG‑AA‑TTaaaaaaaaagaagggaaaa                                                                              >  1:843098/1‑238 (MQ=255)
                                                                                                                                                       ttCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGTGGTACAGGTAGAGCGCTTTCAACGGCT‑GTTGGAGCAATTT‑CCGGGGCAGTTACCGGCAAGAA‑GGGTGGAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGCAATGGTAA‑TCAAAAAAAAA‑GAAGGAAAAg                                                                             >  1:676049/1‑251 (MQ=255)
                                                                                                                                                         cTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGGAATTGGGGGTATTGCCGGGTCTTTGATCGGTGGTGGTAAAGGTAGTGCGCTTGCAACGGGT‑GTTTGAGGAATTG‑CCGGTGCAGTTTCCGGCAGGAAGGGTTG‑AGGAAAAGTCAGCCCAGGTTGTG‑CCTTGGGAATGGGAA‑TCaaaaaataatgaggaaaaaga                                                                            >  1:61642/1‑236 (MQ=255)
                                                                                                                                                         cTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTa                                                                                                     <  1:550377/227‑1 (MQ=255)
                                                                                                                                                         cTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGGGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGGGCAGGTACCGGCAGGAA‑GGGTGAAGAAAAAGGCAGCCAAGGTTGCG‑GCCTGGAAATGGGAA‑tta                                                                                                >  1:893527/1‑229 (MQ=255)
                                                                                                                                                          tATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGaaa                                                                          <  1:516355/251‑1 (MQ=255)
                                                                                                                                                           aTTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTTCAATTGGTGGTGGGACAGGGAGAGCGCCTGCAACGGCT‑GTTGGGGCAATTT‑CCGGTGCAGGTACCGGCAGTAA‑GGGTGAAGAAAAAAGCAGCCCAGTTAGCG‑CCTTGGGAATGGGAA‑TCAAAAAAAAT‑GACGGAAAAGaaaa                                                                         >  1:37921/1‑250 (MQ=255)
                                                                                                                                                               gTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGa                                                                                         <  1:697737/232‑1 (MQ=255)
                                                                                                                                                               gTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTGGGGCAATTG‑CCGGGGCCGGTACCGGCCGTAA‑GGGTGAAGAAAAAGGCAGCCAAGTTAGCG‑CCCTGGGAATGGGAA‑TCAAAAAAAAA‑GAAGGAAAAGAAAatgtt                                                                     >  1:468825/1‑251 (MQ=255)
                                                                                                                                                               gTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGTGGTGGTACAGGGAGAGCGCTTGCAACGGCT‑GTTGGAGGAATTG‑CCGGGGCAGTTAACGGCAGTAA‑GGGTGAAGAAAAAGGCAGCCAAGGTAGCG‑CCTTGGGAATGGGAA‑TCAAAAAAAAT‑GACGGAAAAGAAAttgtt                                                                     >  1:11096/1‑251 (MQ=255)
                                                                                                                                                                 tATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATtgttgt                                                                   <  1:538485/251‑1 (MQ=255)
                                                                                                                                                                   tGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAg                                                                 <  1:403346/251‑1 (MQ=255)
                                                                                                                                                                   tGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGTGGGGGTACAGGGAGAGCGCTTGCAACGGCT‑GTTGGGGCAATTT‑CCGGGGCAGTTACCCGCAGTAA‑GGGTGGAGAAAAAGGCAGCCAAGTTAACG‑CCCTGGGAATGGGAA‑TCAAAAAAAAA‑GACGGAAAAAAAAATGTTGTAg                                                                 >  1:54046/1‑251 (MQ=255)
                                                                                                                                                                         cATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGCATTGCCGGATCTGCAATCGGGGGGGGGTCAGGGAGAGCGCCTGCTACGGGTGGTTGG‑GTAAATG‑GCGGGGCCGGTACCGGCAGTTA‑GAGTGGAGAAAAAGGCCGCCAAATTAGGG‑CCCTGGGAATGGGAATTCAAAAAA‑ATGGAGGGAAAAGAAAtgttttaggggaaa                                                            >  1:429752/1‑238 (MQ=255)
                                                                                                                                                                             gTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCCGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAg                                                                             <  1:495610/229‑1 (MQ=255)
                                                                                                                                                                             gTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGGGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTGGAGCAAATG‑CCGGGGCAGTTACCGGCAGGAA‑GGGTGAAGAAAAAATCAGCCAAGTTAGCG‑CCCTGGGAATGGGAA‑ACAAAAAAAAA‑GACGGAAAAAAAAatgtctcttataaacat                                                         >  1:206329/1‑236 (MQ=255)
                                                                                                                                                                              tttCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑gg                                                                                                                                               <  1:424696/165‑1 (MQ=255)
                                                                                                                                                                                        cccgGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAAc                                             <  1:148417/250‑1 (MQ=255)
                                                                                                                                                                                            gTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAg                                         <  1:29610/250‑1 (MQ=255)
                                                                                                                                                                                                           gATAACCAGGGCGTGATTGGTACTATTGGCGGTGGGGCAATTGGCGGGATTGCCGGATCTGCAATCGGTGGGGGGACAGGGAGAGCGCTTGCAACGGCT‑GGTGGGGCAATTG‑CCGGTGCAGTTACCGGCAGGAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGGAAAGAAATTGTTGTAGGGCAAAAAgg                                                       >  1:612769/1‑220 (MQ=255)
                                                                                                                                                                                                             tAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAAttgt                                                                      <  1:810768/204‑1 (MQ=255)
                                                                                                                                                                                                              aaCCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACg                      <  1:760412/251‑1 (MQ=255)
                                                                                                                                                                                                               aCCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGTGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGt                     <  1:452186/251‑1 (MQ=255)
                                                                                                                                                                                                               aCCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGt                     <  1:114182/251‑1 (MQ=255)
                                                                                                                                                                                                                ccAGGGCGTGATTGGTACTATT‑GCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTAt                   <  1:130422/251‑1 (MQ=255)
                                                                                                                                                                                                                   gggCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGGATTGCCGGATCTGCAATCGGGGGGGGGACAGGGAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGGGCAGTTACCGGc                                                                                                                                                       >  1:569923/1‑121 (MQ=255)
                                                                                                                                                                                                                    ggCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTAtcg                 <  1:93396/250‑1 (MQ=255)
                                                                                                                                                                                                                    ggCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTAt                   <  1:864158/248‑1 (MQ=255)
                                                                                                                                                                                                                    ggCGTGATTGGTACTATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCAAACCTAACTTAGTGCCGGGTGCGCGGGTACGTAt                   <  1:807343/248‑1 (MQ=255)
                                                                                                                                                                                                                                  tATTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGtt  <  1:395702/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aTTGGCGGTGGTGCAATTGGCGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCggg                                  <  1:16673/218‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
TAGGTCTGGTTTAATAACTTCATAATTTTAAGGAGCATTTTATGAAAAAAATCGTATTAGCTGTTTTAGTGGGCTTAACTGTTGCCGGTTGTGCTAATCAAGATGTTTATAGTGGCAGTGTGTATTCTTCGGGACAAGCAAAAGAAGCGCGTTCTATTAGTTATGGTACCATCGTTTCAACCCGCCCGGTGAAAATTCAGGCAGATAACCAGGGCGTGATTGGTACTATTGGCGGTGGTGCAATTGGTGGTATTGCCGGATCTGCAATCGGTGGTGGTACAGGTAGAGCGCTTGCAACGGCT‑GTTGGAGCAATTG‑CCGGTGCAGTTACCGGCAGTAA‑GGTTGAAGAAAAAGTCAGCCAAGTTAGCG‑CCTTGGAAATGGTAA‑TCAAAAAAGAT‑GACGGAAAAGAAATTGTTGTAGTGCAAAAAGCCGAACCTAACTTAGTGCCGGGTGCGCGGGTACGTATCGTCGGCGGTTCACGTT  >  1/1018670‑1019146

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