Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 593,330 T→C intergenic (+201/+216) unknown → / ← unknown hypothetical protein/UPF0758 protein YicR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1593,3300TC90.9% 162.9 / 4.7 66intergenic (+201/+216)unknown/unknownhypothetical protein/UPF0758 protein YicR
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base C (21/39);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GGGGCGTTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGACAGTTAAAAATAAGTTGGCTGCCGGCAATTAAAAGGAGTAAGATTTAAAACGCCAGAGCAATCTAGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAA  >  1/593085‑593488
                                                                                                                                                                                                                                                     |                                                                                                                                                              
ttagatttatttgagcaagctgtagcggaaaaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTAc                                                                                                                                                           <  1:168102‑M2/140‑1 (MQ=255)
  agatttatttgagcaagctgtagcggaaaaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACt                                                                                                                                                          >  1:348389‑M2/110‑250 (MQ=255)
    atttatttgagcaagctgtagcggaaaaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTActc                                                                                                                                                         <  1:858124‑M2/142‑2 (MQ=255)
    atttatttgagcaagctgtagcggaaaaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTActc                                                                                                                                                         <  1:590268‑M2/142‑2 (MQ=255)
        atttgagcaagctgtagcggaaaaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCAc                                                                                                                                                   <  1:23606‑M2/148‑1 (MQ=255)
                   ctgtagcggaaaaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAg                                                                                                                                         <  1:802821‑M2/158‑1 (MQ=255)
                              aaagattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAAc                                                                                                                             <  1:21087‑M2/170‑1 (MQ=255)
                                 gattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAAcgg                                                                                                                           >  1:489507‑M2/79‑250 (MQ=255)
                                 gattaatgatagttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCCGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAacggg                                                                                                                          >  1:788641‑M2/79‑250 (MQ=255)
                                             ttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAcc                                                                                                               <  1:505072‑M2/184‑1 (MQ=255)
                                             ttaaaaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAc                                                                                                                <  1:460813‑M2/183‑1 (MQ=255)
                                                 aaataagttgactgccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGATCAGGCAGACAAAAACGGCGGGGTAATACCcgtcg                                                                                                          <  1:388202‑M2/189‑1 (MQ=255)
                                                               ccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTtgtttatgtt                                                                                            <  1:662066‑M2/203‑1 (MQ=255)
                                                               ccggcaattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTtgtttatgtt                                                                                            <  1:823710‑M2/203‑1 (MQ=255)
                                                                   caattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTtata                                                                                        <  1:639172‑M2/207‑1 (MQ=255)
                                                                    aattaaaaggattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTatat                                                                                       <  1:295208‑M2/208‑1 (MQ=255)
                                                                      ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTAc                                                                                                                                                           >  1:455571/10‑175 (MQ=255)
                                                                      ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTAc                                                                                                                                                           <  1:409961/166‑1 (MQ=255)
                                                                      ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTCATTTACTGCCGAATAGGCAgcttagaaa                                                                                                                             >  1:104402/10‑197 (MQ=255)
                                                                             gattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATtc                                                                               <  1:179525‑M2/216‑1 (MQ=255)
                                                                             gattaagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCGTCGTTTTTTTTTATGTTTATATGATGAtt                                                                                >  1:716652‑M2/35‑249 (MQ=255)
                                                                                aaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTCATTTACTGCCGAATAGGCAgcttagaaa                                                                                                                             >  1:462206/1‑187 (MQ=255)
                                                                                 aagatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGt                                                                          <  1:442944‑M2/221‑1 (MQ=255)
                                                                                   gatttaaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTtgtttatg                                                                                              <  1:224828‑M2/201‑1 (MQ=255)
                                                                                       taaaatgccagagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAAtt                                                                    <  1:436222‑M2/227‑1 (MQ=255)
                                                                                                 gagcaatctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCGTCGTTTTTGTTTTTGTTTATATGATGATTTTGGGTTAAATttgagggcgat                                                          >  1:421794‑M2/15‑247 (MQ=255)
                                                                                                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTtaaatttgagggcgattagtc                                                     <  1:681563‑M2/242‑10 (MQ=255)
                                                                                                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTtaaatttgagggcgattagtc                                                     <  1:824457‑M2/242‑10 (MQ=255)
                                                                                                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAAt                                                                     <  1:290052‑M2/226‑1 (MQ=255)
                                                                                                      atctggcggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCCTCGGTTTTGTTTTTGTTTTTATTATGGTTCTGCGGTtaattttagggcgattagt                                                      >  1:279353‑M2/10‑242 (MQ=255)
                                                                                                            cggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTctgcgttaaatttgagggcgattagtcgttatt                                               <  1:859078‑M2/248‑16 (MQ=255)
                                                                                                            cggTGGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTctgcgttaaatttgagggcgattagtcgttatt                                               <  1:321418‑M2/248‑16 (MQ=255)
                                                                                                            aTCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACtc                                                                                                                                                         >  1:9770/1‑138 (MQ=255)
                                                                                                              gtggCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttattc                                              <  1:603904/247‑17 (MQ=255)
                                                                                                               tggCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttattca                                             <  1:767941/248‑18 (MQ=255)
                                                                                                                    gTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCAc                                                                                                                                                   >  1:465789/1‑143 (MQ=255)
                                                                                                                    gTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttatt                                               <  1:462008/243‑16 (MQ=255)
                                                                                                                    gTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGAGGGAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAgggcgattagtcgttatt                                               >  1:432932/1‑228 (MQ=255)
                                                                                                                              gCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAggg                                                              <  1:691838/218‑1 (MQ=255)
                                                                                                                              gCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCGTCGTTTTTGTTTTTGTTTttat                                                                                       >  1:37977/1‑193 (MQ=255)
                                                                                                                              gCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTTAgggcgattagtcgttatt                                               >  1:647138/1‑218 (MQ=255)
                                                                                                                                cTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGt                                                                          <  1:13219/204‑1 (MQ=255)
                                                                                                                                cTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAg                           >  1:793873/1‑251 (MQ=255)
                                                                                                                                      ggCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTtgtttatg                                                                                              >  1:56762/1‑178 (MQ=255)
                                                                                                                                        cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                   <  1:622206/251‑5 (MQ=255)
                                                                                                                                        cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                   <  1:126024/251‑5 (MQ=255)
                                                                                                                                             tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta              <  1:821071/251‑10 (MQ=255)
                                                                                                                                             tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta              <  1:783339/251‑10 (MQ=255)
                                                                                                                                                   aGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTTTTTATGTTTATATGATGATTCTGCGTTAAATTTTAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAgtgt                        >  1:322279/1‑235 (MQ=255)
                                                                                                                                                   aGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTTAgggcgattagtcgttatt                                               >  1:348518/1‑197 (MQ=255)
                                                                                                                                                     caaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattc      <  1:543917/251‑18 (MQ=255)
                                                                                                                                                      aaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAACCGTCGCTTTTGTTTTTGTTTATATGATGATTCTGCGTTAAATTTGAgggcgatttgtcgttatt                                               >  1:649193/1‑194 (MQ=255)
                                                                                                                                                      aaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCATCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTTAGGGGGGTTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                   >  1:158898/1‑233 (MQ=255)
                                                                                                                                                         caaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattcatac  <  1:859192/251‑22 (MQ=255)
                                                                                                                                                            gacgtAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTCATTTACTGCCGaat                                                                                                                                           <  1:100385/107‑2 (MQ=255)
                                                                                                                                                                   ctgctgAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAc                                                                                                                >  1:642986/1‑131 (MQ=255)
                                                                                                                                                                     gctgAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaatta         <  1:482745/232‑15 (MQ=255)
                                                                                                                                                                        gAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTTTGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                   >  1:826309/1‑215 (MQ=255)
                                                                                                                                                                                 gcTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATa                                                                                                                 <  1:460553/116‑1 (MQ=255)
                                                                                                                                                                                 gcTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATatgat                                                                                    >  1:155728/1‑145 (MQ=255)
                                                                                                                                                                                    tAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATatgat                                                                                    <  1:51227/142‑1 (MQ=255)
                                                                                                                                                                                           tGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgtt                    >  1:775368/1‑196 (MQ=255)
                                                                                                                                                                                                                    gCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta              >  1:440256/1‑171 (MQ=255)
                                                                                                                                                                                                                                       aaaTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                   <  1:258670/156‑5 (MQ=255)
                                                                                                                                                                                                                                         aTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat           <  1:284319/162‑13 (MQ=255)
                                                                                                                                                                                                                                               gtTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATtc                                                                               <  1:191995/88‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                              
GGGGCGTTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGACAGTTAAAAATAAGTTGGCTGCCGGCAATTAAAAGGAGTAAGATTTAAAACGCCAGAGCAATCTAGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAA  >  1/593085‑593488

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