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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1593,1830CT93.0% 214.2 / 3.1 86intergenic (+54/+363)unknown/unknownhypothetical protein/UPF0758 protein YicR
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base T (38/42);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

ATGACATTGCGGACATTATCGCTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTTTACAAGATGCTGGAGCTAAATTTAATAGAGATGGGGCGTTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGACAGTTAAAAATAAGTTGGCTGCCGGCAATTAAAAGGAGTAAGATTTAAAACGCCAGAGCAATCTAGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGG  >  1/592948‑593428
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
atgaCATTGCGGACATTATCGCTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtgg                                                                                                                                                                                                                                        >  1:260116‑M1/1‑248 (MQ=255)
             cATTATCGCTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggtttt                                                                                                                                                                                                                           >  1:355614‑M1/1‑235 (MQ=255)
                 aTCGCTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgc                                                                                                                                                                                                                       <  1:488391‑M1/251‑21 (MQ=255)
                  tCGCTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGc                                                                                                                                                                                                                                               <  1:645136/226‑1 (MQ=255)
                     cTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTCCAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTCATGATAGTTAAAAATACGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagc                                                                                                                                                                                                                    >  1:222643‑M1/1‑227 (MQ=255)
                       gCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgc                                                                                                                                                                                                                       >  1:240334‑M1/1‑225 (MQ=255)
                          gATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaa                                                                                                                                                                                                              >  1:278964‑M1/1‑222 (MQ=255)
                               cGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgc                                                                                                                                                                                                                       <  1:309368‑M1/237‑21 (MQ=255)
                                 cTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgag                                                                                                                                                                                                       >  1:366439‑M1/1‑215 (MQ=255)
                                 cTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctg                                                                                                                                                                                                         >  1:144640‑M1/1‑215 (MQ=255)
                                          cATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaag                                                                                                                                                                                               <  1:238959‑M1/250‑45 (MQ=255)
                                          cATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaag                                                                                                                                                                                               <  1:547566‑M1/250‑45 (MQ=255)
                                                aagaaAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgc                                                                                                                                                                                                                       >  1:537335‑M1/1‑200 (MQ=255)
                                                    aaGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaag                                                                                                                                                                                      >  1:811874‑M1/1‑196 (MQ=255)
                                                         ggTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaa                                                                                                                                                                                 >  1:111492‑M1/1‑191 (MQ=255)
                                                           tATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCt                                                                                                                                                                                                                                          <  1:77592‑M1/190‑2 (MQ=255)
                                                                   cTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataag                                                                                                                                                                     >  1:464298‑M1/1‑181 (MQ=255)
                                                                    tATTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATGTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagc                                                                                                                                                                    <  1:450445‑M1/251‑72 (MQ=255)
                                                                     aTTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagca                                                                                                                                                                   >  1:557638‑M1/1‑179 (MQ=255)
                                                                     aTTCGAGAGAATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagca                                                                                                                                                                   >  1:869732‑M1/1‑179 (MQ=255)
                                                                            agaATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttag                                                                                                                                                             <  1:15288‑M1/250‑79 (MQ=255)
                                                                             gaATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgc                                                                                                                                                                             >  1:746604‑M1/1‑171 (MQ=255)
                                                                             gaATGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaa                                                                                                                                                                                        <  1:136009‑M1/222‑52 (MQ=255)
                                                                              aaTGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaa                                                                                                                                                          >  1:709535‑M1/1‑170 (MQ=255)
                                                                              aaTGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaa                                                                                                                                                          <  1:830664‑M1/251‑82 (MQ=255)
                                                                              aaTGTGCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaa                                                                                                                                                          <  1:793022‑M1/251‑82 (MQ=255)
                                                                                 gtgCAGACGCACACTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataag                                                                                                                                                                     <  1:498614‑M1/237‑71 (MQ=255)
                                                                                   aagagacatcacttggattttCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggt                                                                                                                                                     <  1:412469‑M1/242‑87 (MQ=255)
                                                                                            cacTTGGATTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagc                                                                                                                                                                                     >  1:571215‑M1/1‑156 (MQ=255)
                                                                                                   aTTTTCTACAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATGAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttacggaaag                                                                                                                                      <  1:448613‑M1/250‑102 (MQ=255)
                                                                                                        ctaCAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttca                                                                                                                                  >  1:652735‑M1/2‑144 (MQ=255)
                                                                                                        ctaCAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaa                                                                                                                                                          >  1:368206‑M1/2‑144 (MQ=255)
                                                                                                        ctaCAAGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctaaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagc                                                                                                                                                                    >  1:453487‑M1/2‑144 (MQ=255)
                                                                                                             aGATGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcct                                                                                                                                                                                                                  <  1:435185‑M1/164‑26 (MQ=255)
                                                                                                               aTGCGGGAGCTAAATTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaac                                                                                                                         <  1:729806‑M1/251‑115 (MQ=255)
                                                                                                                          aaaTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCagag                                                                                                                                                                                                                                                >  1:484286/1‑121 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaa                                                                                                            >  1:74929‑M1/1‑124 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaa                                                                                                            >  1:367057‑M1/1‑124 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaa                                                                                                            >  1:28088‑M1/1‑124 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaa                                                                                                            >  1:513870‑M1/1‑124 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctc                                                                                                                                                >  1:84616‑M1/1‑124 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctg                                                                                                                                                                           <  1:476758‑M1/188‑65 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaag                                                                                                                                                                                               <  1:816764‑M1/168‑45 (MQ=255)
                                                                                                                            aTTTAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaagg                                                                                                                                                                                                            >  1:676878‑M1/1‑124 (MQ=255)
                                                                                                                             tttAATAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaa                                                                                                           <  1:730728‑M1/251‑129 (MQ=255)
                                                                                                                                 aTAGAGATGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatag                                                                                                       <  1:720344‑M1/251‑133 (MQ=255)
                                                                                                                                     agaTGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataa                                                                                                                                                                      >  1:182497‑M1/1‑115 (MQ=255)
                                                                                                                                      gaTGAGGCATTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaac                                                                                                                         >  1:419665‑M1/1‑114 (MQ=255)
                                                                                                                                            gcattagaTTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagct                                                                                             <  1:379077‑M1/247‑143 (MQ=255)
                                                                                                                                            gcattagaTTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaat                                                                                                         >  1:251016‑M1/4‑108 (MQ=255)
                                                                                                                                             cattagATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctc                                                                                            <  1:156422‑M1/248‑144 (MQ=255)
                                                                                                                                             cattagATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaa                                                                                                                                        <  1:115508‑M1/204‑100 (MQ=255)
                                                                                                                                               ttAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctac                                                                                         <  1:168102‑M1/251‑147 (MQ=255)
                                                                                                                                                 aGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctact                                                                                        >  1:348389‑M1/1‑103 (MQ=255)
                                                                                                                                                   atttatttGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactc                                                                                       <  1:590268‑M1/249‑149 (MQ=255)
                                                                                                                                                   atttatttGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactc                                                                                       <  1:858124‑M1/249‑149 (MQ=255)
                                                                                                                                                       atttGAGCAAGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcac                                                                                 <  1:23606‑M1/251‑155 (MQ=255)
                                                                                                                                                                aGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagc                                                                                              >  1:24663‑M1/1‑88 (MQ=255)
                                                                                                                                                                aGCTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgc                                                                                                                              >  1:719171‑M1/1‑88 (MQ=255)
                                                                                                                                                                  cTGTAGCGGAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacag                                                                       <  1:802821‑M1/250‑165 (MQ=255)
                                                                                                                                                                         ggAAAAAGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaac                                                                                                                         >  1:137816‑M1/1‑79 (MQ=255)
                                                                                                                                                                             aaaGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaac                                                           <  1:21087‑M1/251‑177 (MQ=255)
                                                                                                                                                                               aGATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagc                                                                                                                  >  1:234151‑M1/1‑73 (MQ=255)
                                                                                                                                                                                gATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacgg                                                         >  1:489507‑M1/1‑72 (MQ=255)
                                                                                                                                                                                gATTAATGATAGTTAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggccgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggg                                                        >  1:788641‑M1/1‑72 (MQ=255)
                                                                                                                                                                                            ttAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacc                                             <  1:505072‑M1/250‑191 (MQ=255)
                                                                                                                                                                                            ttAAAAATAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatac                                              <  1:460813‑M1/249‑190 (MQ=255)
                                                                                                                                                                                                aaaTAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgatcaggcagacaaaaacggcggggtaatacccgtcg                                        <  1:388202‑M1/251‑196 (MQ=255)
                                                                                                                                                                                                 aaTAAGTTGACTGCCGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaa                                                                                                                                                          >  1:292646‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                          gaccgttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAAc                                                                                                                               >  1:125852/15‑148 (MQ=255)
                                                                                                                                                                                                              ccGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtt                          <  1:823710‑M1/251‑210 (MQ=255)
                                                                                                                                                                                                              ccGGCAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtt                          <  1:662066‑M1/251‑210 (MQ=255)
                                                                                                                                                                                                               ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGc                                                                                                                           >  1:426702/10‑147 (MQ=255)
                                                                                                                                                                                                               ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTAc                                                                                               <  1:409961/166‑1 (MQ=255)
                                                                                                                                                                                                               ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTAc                                                                                               >  1:455571/10‑175 (MQ=255)
                                                                                                                                                                                                               ttactaatcaaaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTCATTTACTGCCGAATAGGCAgcttagaaa                                                                 >  1:104402/10‑197 (MQ=255)
                                                                                                                                                                                                                  cAATTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtttata                      <  1:639172‑M1/251‑214 (MQ=255)
                                                                                                                                                                                                                   aaTTAAAAGGATTAAGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtttatat                     <  1:295208‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                         aaGGAGTAAGATTCAAACCGTCCGAGCAATCTGGCGGTTTTATGCAGCC‑GAAAGTCTGAGTGTTG‑‑‑‑AAACAGA‑CGTAACTGCTGAACAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTCATTTACTGCCGAATAGGCAgcttagaaa                                                                 >  1:462206/1‑187 (MQ=255)
                                                                                                                                                                                                                           ggattaagaTTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaa                                                                                                           <  1:383722‑M1/153‑129 (MQ=255)
                                                                                                                                                                                                                           ggattaagaTGTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtag                                                                                                    <  1:518242‑M1/160‑136 (MQ=255)
                                                                                                                                                                                                                            gattaagATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtttatatgatgattc             <  1:179525‑M1/247‑223 (MQ=255)
                                                                                                                                                                                                                            gattaagATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcgggggaatacccgtcgtttttttttatgtttatatgatgatt              >  1:716652‑M1/4‑28 (MQ=255)
                                                                                                                                                                                                                                aaGATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtttatatgatgattctgcgt        <  1:442944‑M1/251‑228 (MQ=255)
                                                                                                                                                                                                                                  gATTTAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatg                            <  1:224828‑M1/229‑208 (MQ=255)
                                                                                                                                                                                                                                      tAAAATGCCAGAGCAATCtggcggtggcggttttatgcagcctgaaaggctgagtgttgaagcaaacaaagcttaactgctgaataagcagcttagaaaatggtagctcttactgaaagttcactgccgaacaggcagcttagaaaatagtagttaagctctactcattcactgccgaacaggcagacaaaaacggcggggtaatacccgtcgtttttgtttatgtttatatgatgattctgcgttaaatt  <  1:436222‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
ATGACATTGCGGACATTATCGCTGCAGATTACGCTTTAATCGCATGCGAAGAAAGCGGGTATGATGACTATTCGAGAGAATGTGCAGACGCACACTTGGATTTTTTACAAGATGCTGGAGCTAAATTTAATAGAGATGGGGCGTTAGATTTATTTGAGCAAGCTGTAGCGGAAAAAGATTAATGACAGTTAAAAATAAGTTGGCTGCCGGCAATTAAAAGGAGTAAGATTTAAAACGCCAGAGCAATCTAGCGGTTTTATGCAGCCTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATACCCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGG  >  1/592948‑593428

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