Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 581,414 G→A S125N (AGC→AAC)  unknown → SSU ribosomal protein S1p

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1581,4140GA96.5% 145.0 / NA 59S125N (AGC→AAC) unknownSSU ribosomal protein S1p
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/27);  minor base C (0/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCCGAATCAGCTATTCCTGCCGAAGAATTTTTTAATGCACAAGGTGAATTGGAAGTTCAAGTGGGTGATGTTGTTGATGTCGTGCTTAAAGCGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAGCGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAGGTACCGTTAAAAACTTAAC  >  1/581167‑581648
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
aaaTCCGAATCAGCTATTCCTGCCGAAGAATTTTTTAATGCACAAGGTGAATTGGAAGTTCAAGTGGGTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACgg                                                                                                                                                                                                                                          <  1:296163/251‑1 (MQ=255)
                                    aaTGCACAAGGTGAATTGGAAGTTCAAGTGGGTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGATCAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATAc                                                                                                                                                                                                      <  1:347937/251‑1 (MQ=255)
                                          caAGGTGAATTGGAAGTTCAAGTGGGTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTcc                                                                                                                                                                                                <  1:814076/251‑1 (MQ=255)
                                             ggTGAATTGGAAGTTCAAGTGGGTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACgg                                                                                                                                                                                                                                          <  1:300373/206‑1 (MQ=255)
                                                                 gggTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATc                                                                                                                                                                               <  1:14357/245‑1 (MQ=255)
                                                                 gggTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCactta                                                                                                                                                                          >  1:249546/1‑250 (MQ=255)
                                                                  ggTGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCacttac                                                                                                                                                                         <  1:670721/250‑1 (MQ=255)
                                                                    tGATGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCacttactt                                                                                                                                                                       <  1:864596/250‑1 (MQ=255)
                                                                      aTGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTgg                                                                                                                                                                     >  1:210566/1‑250 (MQ=255)
                                                                      aTGTTGTTGATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTgg                                                                                                                                                                     >  1:679917/1‑250 (MQ=255)
                                                                         ttgttgATGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCaaa                                                                                                                                                                 >  1:545091/1‑251 (MQ=255)
                                                                               aTGTCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGtt                                                                                                                                                                                                                                      <  1:519744/176‑1 (MQ=255)
                                                                                  tCGTGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAagaattagaa                                                                                                                                                        <  1:149041/251‑1 (MQ=255)
                                                                                     tGCTTAAAGTGGTGGATGACTGTTTCGGTGAATCTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCAGACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTc                                                                                                                                                     <  1:702876/251‑1 (MQ=255)
                                                                                     tGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAagaattagaatt                                                                                                                                                      >  1:338624/1‑250 (MQ=255)
                                                                                     tGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAagaattagaatt                                                                                                                                                      >  1:248570/1‑250 (MQ=255)
                                                                                     tGCTTAAAGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAagaattagaatt                                                                                                                                                      >  1:135292/1‑250 (MQ=255)
                                                                                          aaaGTGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGt                                                                                                                                                <  1:315335/251‑1 (MQ=255)
                                                                                           aagtggTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTa                                                                                                                                               >  1:161266/5‑251 (MQ=255)
                                                                                                   gATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAAc                                                                                                                                        <  1:494485/250‑1 (MQ=255)
                                                                                                    aTGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACt                                                                                                                                       <  1:847408/250‑1 (MQ=255)
                                                                                                    aTGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACt                                                                                                                                       <  1:112740/250‑1 (MQ=255)
                                                                                                                                  gTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACgttgttg                                                                                                             <  1:124960/246‑1 (MQ=255)
                                                                                                                                               aCGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACgttg                                                                                                                >  1:605693/1‑230 (MQ=255)
                                                                                                                                                        aGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACt                                                                                                                                       <  1:430843/198‑1 (MQ=255)
                                                                                                                                                        aGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCg                                                                                        >  1:321002/1‑245 (MQ=255)
                                                                                                                                                         gCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCgccgtgccg                                                                                              <  1:647459/238‑1 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCacttac                                                                                                                                                                         >  1:60504/1‑154 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACgt                                                                                                                  >  1:695730/1‑209 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCgccgtgccg                                                                                              <  1:333164/229‑1 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCgccgtgcc                                                                                               >  1:544527/1‑228 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCgccgtgc                                                                                                >  1:558040/1‑227 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGaaa                                                                               <  1:143139/244‑1 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCa                                                                        >  1:876330/1‑251 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCa                                                                        >  1:301302/1‑251 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCa                                                                        >  1:302303/1‑251 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCa                                                                        >  1:351169/1‑251 (MQ=255)
                                                                                                                                                                  gTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGACGCAGATc                                                                                                                                                                               <  1:297059/148‑1 (MQ=255)
                                                                                                                                                                       ggAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGaaa                                                                               >  1:331408/1‑239 (MQ=255)
                                                                                                                                                                                     acgaACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGTTGGCTTCACAGTTGAGTTAATCGGTGTTCGCGCATTCTTACCGGTTTCATTGGTTGATCCCCGTCCTTTTCGTGATGCAAATCACTTACTTGGCAAAGAATTAGAATTCAAAGAAATCAAACTTGATCAAAAACGTAACAACTTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGAt                                                                    <  1:18536/236‑1 (MQ=255)
                                                                                                                                                                                           aaGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAACCGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAAc                                              <  1:203498/251‑1 (MQ=255)
                                                                                                                                                                                                     tCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTgg                                                   >  1:693713/1‑236 (MQ=255)
                                                                                                                                                                                                               gttAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAg                           <  1:317969/250‑1 (MQ=255)
                                                                                                                                                                                                               gttAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAg                           <  1:208609/250‑1 (MQ=255)
                                                                                                                                                                                                                      gTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAACCGTAACAACGTTTTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAgg                    <  1:542603/250‑1 (MQ=255)
                                                                                                                                                                                                                       tAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAGGt                   <  1:469674/250‑1 (MQ=255)
                                                                                                                                                                                                                          aGTTAAAGGTGGCTTCACCGTTGAGTTAAACGGTATTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGAAAACTTACTTGGGAACAAATTAGAATTCAACGTAATCAAACTTGATCACAAACGTAACAGCGTGATTGTTTCCCGAAGTGCCGTTATAGCATACGAAAGCAGCCACGATCGCGAAGCAAG‑ATTGGCGAAATTGTCTGAAGGTGCAGAAGTTAAAGGTAc                 >  1:867524/1‑249 (MQ=255)
                                                                                                                                                                                                                          aGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACt                                                                                                                                       >  1:458350/1‑132 (MQ=255)
                                                                                                                                                                                                                          aGTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAg                                                                      >  1:590279/1‑197 (MQ=255)
                                                                                                                                                                                                                           gTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACgttgttg                                                                                                             >  1:789301/1‑157 (MQ=255)
                                                                                                                                                                                                                           gTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAgg                    <  1:25057/245‑1 (MQ=255)
                                                                                                                                                                                                                           gTTAAAGGTGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAGGTACCGt              >  1:753690/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tGGCTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAGGTACCGTTAAAAAc       <  1:32594/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cTTCACAGTTGAGTTAACCGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATAAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAgg                    <  1:809134/234‑1 (MQ=255)
                                                                                                                                                                                                                                      cTTCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACgt                                                                                                                  >  1:358124/1‑141 (MQ=255)
                                                                                                                                                                                                                                        tCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAgg                                   >  1:239296/1‑217 (MQ=255)
                                                                                                                                                                                                                                        tCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAg                                    >  1:632941/1‑216 (MQ=255)
                                                                                                                                                                                                                                        tCACAGTTGAGTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAGGTACCGTTAAAAACTTAAc  >  1:237852/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gTTAAACGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACg                                                                                                                   <  1:680848/128‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           
AAATCCGAATCAGCTATTCCTGCCGAAGAATTTTTTAATGCACAAGGTGAATTGGAAGTTCAAGTGGGTGATGTTGTTGATGTCGTGCTTAAAGCGGTGGATGACGGTTTCGGTGAAACTGTAGCTTCCCGTGCTGATGCAAAACGCAACGAAGCTTGGATTGTATTGGAAAAAGCATACGACGAACAAGCGACTGTTCTTGGTTTGGTTAACGGTAAAGTTAAAGGTGGCTTCACAGTTGAGTTAAGCGGTGTTCGCGCATTCTTACCGGGTTCATTGGTTGATACCCGTCCTGTTCGTGATGCAGATCACTTACTTGGCAAAGAATTAGAATTCAAAGTAATCAAACTTGATCAAAAACGTAACAACGTTGTTGTTTCCCGCCGTGCCGTTATCGAATCCGAAAGCAGCCAAGATCGCGAAG‑AAGTATTGGCGAACTTGTCTGAAGGTGCAGAAGTTAAAGGTACCGTTAAAAACTTAAC  >  1/581167‑581648

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