Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 536,492 A→G intergenic (+6/‑94) hyi → / → glxR hydroxypyruvate isomerase/tartronate semialdehyde reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913536,4920AG97.4% 273.4 / NA 79intergenic (+6/‑94)hyi/glxRhydroxypyruvate isomerase/tartronate semialdehyde reductase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (44/32);  minor base T (0/2);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTAACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCTGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTATGGGTACACCGATGGCCATTAATCTGGCGCGTGCCGGTCATCAATTACATGTCACGACCATTGGACCGGTTGCTGATGAATTACTGTCACTGGGTGCCGTCAGTGTTGAAACTGCT  >  NC_000913/536248‑536732
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
aTCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAACCGTGACGCt                                                                                                                                                                                                                                             >  2:321388/1‑250 (MQ=255)
aTCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTa                                                                                                                                                                                                                                            <  1:81216/251‑1 (MQ=255)
  cATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTAtt                                                                                                                                                                                                                                          >  1:838348/1‑251 (MQ=255)
   atatGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTc                                                                                                                                                                                                                                         >  1:447570/1‑251 (MQ=255)
          gCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAAt                                                                                                                                                                                                                                  >  1:393553/1‑251 (MQ=255)
            ggATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGc                                                                                                                                                                                                                                >  2:832822/1‑251 (MQ=255)
              aTGGACGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTCCGCTGGATGGATCCACACCGTTAAAACGGGACGCTCATCAGACAATGCtt                                                                                                                                                                                                                              >  2:251427/1‑251 (MQ=255)
              aTGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCtt                                                                                                                                                                                                                              >  1:342539/1‑251 (MQ=255)
                  aaGGTGAATTAACCAACACCATGACTCAGTGGGCAGATAAAATTGGTAACCTGCAAATTGCAGCTAATCCGCATCGCGGCGATCCGGGAACCGGAGAAATTAATTATGCTTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTATACGCTGGATGTATCCGTACCGTTAAAACGTCACGCTATTCAGACAATGCtttttt                                                                                                                                                                                                                          >  2:883918/1‑251 (MQ=255)
                     gTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTc                                                                                                                                                                                                                                         >  2:86304/1‑233 (MQ=255)
                     gTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTc                                                                                                                                                                                                                                         <  1:86304/233‑1 (MQ=255)
                              ccatcacCATTACTCAGTGGTCTTATAAAATTGGTCACCTGCTAATTGCCGATAATCCGCATCGAGGCGAAACGGGAACCGGAGAAATTAATTATGATTATCGCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGTGTGGGTGTGAATATAACCCCCAAACCACCACTGAAGCCGTTTTACGCTGGATGGCTCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGATAAGt                                                                                                                                                                                                              <  2:511330/247‑1 (MQ=255)
                                  cacCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTgg                                                                                                                                                                                                           >  1:863612/1‑250 (MQ=255)
                                     cATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAgg                                                                                                                                                                                                       >  1:842682/1‑251 (MQ=255)
                                          cTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATc                                                                                                                                                                                                  >  1:140349/1‑251 (MQ=255)
                                           tCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCg                                                                                                                                                                                                 >  2:662805/1‑251 (MQ=255)
                                                gggCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGtgtt                                                                                                                                                                                             <  1:417723/250‑1 (MQ=255)
                                                        aaaaTTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAAtt                                                                                                                                                                                    <  1:388603/251‑1 (MQ=255)
                                                                 cACCTGCAAATTGCCTATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGGCCGAATTCAGGAAAAg                                                                                                                                                                           <  2:131267/251‑1 (MQ=255)
                                                                   ccTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCg                                                                                                                                                                         <  2:697226/251‑1 (MQ=255)
                                                                     tGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGa                                                                                                                                                                        >  2:145681/1‑250 (MQ=255)
                                                                      gCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGGCAATGCGTTTGTAGGCCGCTAAGTTGGCAGGGGATCTTGTTGTCCGAATTCAGGAAAAGCGaaa                                                                                                                                                                      <  2:641117/251‑1 (MQ=255)
                                                                       cAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAAt                                                                                                                                                                     <  1:625350/251‑1 (MQ=255)
                                                                        aaaTTGCCGATAATCCGCATCGCGGCGACCCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAAtt                                                                                                                                                                    >  2:91652/1‑251 (MQ=255)
                                                                             gCCGATAATCCGCAGCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTaaaa                                                                                                                                                               <  1:882063/251‑1 (MQ=255)
                                                                              ccGATAATCCGCGTCGCGGCGAACAGGTAAAGGGAGAAATTAATTATGATTATCTCATTAAGGCAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACCATGCATTTTTAGGCCGCTAAGTTGGCATGGGATCGTGTTGTACGAATTCAGGAAAAGCGAAAtttcaa                                                                                                                                                                >  2:328200/1‑246 (MQ=255)
                                                                                gATAATCCGCATAGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAag                                                                                                                                                              >  1:45104/1‑249 (MQ=255)
                                                                                   aaTCCGCATCGCGACGCACCGGCACCCGCAGAAACTAATTCAGATTATCCCTTTAAGGGCAACGAAAATTCTGACTCCAACAGTTGAGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGATAACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGGTGGCAGGGGGTCGTGTTGTCCGAATTCAGGCAAAGCGAAATTTAAAAGAGGtt                                                                                                                                                         >  2:704872/1‑251 (MQ=255)
                                                                                      ccGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAAt                                                                                                                                                      <  2:393553/251‑1 (MQ=255)
                                                                                       cGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAAtt                                                                                                                                                     <  1:567904/251‑1 (MQ=255)
                                                                                            cgcgGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAAGCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGa                                                                                                                                                 <  2:20297/250‑1 (MQ=255)
                                                                                               ggcgAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACt                                                                                                                                             <  2:354165/251‑1 (MQ=255)
                                                                                                  gaaccgggaaccggAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTggg                                                                                                                                          >  1:432779/1‑251 (MQ=255)
                                                                                                     ccgggaaccggAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGAtt                                                                                                                                       <  2:342539/251‑1 (MQ=255)
                                                                                                       gggaaccggAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGAttt                                                                                                                                      >  1:398087/1‑250 (MQ=255)
                                                                                                        ggaaccggAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTAt                                                                                                                                    <  1:662805/251‑1 (MQ=255)
                                                                                                        ggaaccggAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTAt                                                                                                                                    >  1:247584/1‑251 (MQ=255)
                                                                                                         gaaccggAGAAATTAATTATGATTTATTCTTTAAGGTAATCGAAAATTCTTACTACCACGGGTGGGTTTGGGGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTAtt                                                                                                                                   <  2:842682/251‑1 (MQ=255)
                                                                                                            ccggAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGc                                                                                                                                >  1:702758/1‑251 (MQ=255)
                                                                                                                  aaaTTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGc                                                                                                                          <  1:529059/251‑1 (MQ=255)
                                                                                                                  aaaTTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGc                                                                                                                          >  1:806290/1‑251 (MQ=255)
                                                                                                                  aaaTTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGc                                                                                                                          >  1:583332/1‑251 (MQ=255)
                                                                                                                  aaaTTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGc                                                                                                                          >  1:568043/1‑251 (MQ=255)
                                                                                                                        aTTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGGTTGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATGGGCTTAGGCATTAgg                                                                                                                    >  2:164502‑M1/1‑249 (MQ=255)
                                                                                                                                   tctcTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgagatt                                                                                                          >  2:566924‑M1/1‑239 (MQ=255)
                                                                                                                                   tctcTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcg                                                                                                               <  1:566924‑M1/246‑8 (MQ=255)
                                                                                                                                     gcttTAAGGTAATCGAAAATGCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGTATCCGTACCGTTAAAACGTTACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcg                                                                                                         <  2:32313‑M1/249‑14 (MQ=255)
                                                                                                                                         tAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcatt                                                                                                       >  1:196592‑M1/1‑233 (MQ=255)
                                                                                                                                            ggTAATCGAAACTTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTACGAATTCAGGACAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacat                                                                                                 >  1:683290‑M1/1‑230 (MQ=255)
                                                                                                                                             gTAATCGAAAATTCTGACTACAACGGTGGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacagt                                                                                                <  2:692671‑M1/251‑23 (MQ=255)
                                                                                                                                              tAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacattt                                                                                               >  1:821347‑M1/1‑228 (MQ=255)
                                                                                                                                               aaTCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGATGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaagcattagtcatttt                                                                                              >  2:131590‑M1/1‑227 (MQ=255)
                                                                                                                                                    aaaaTTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcct                                                                                         >  1:329799‑M1/1‑222 (MQ=255)
                                                                                                                                                       aTTCTGCCTACAAAGGTTGGGCAGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACCATGCTTTTTTAGGCCGCTTAGTTGGCAGGGGATCGTGTTTTACGAATTCAGGAAATGCGAATTTTAAAAGAGGTTAATTATGACACTGGGTTTTATTGGCTTAGGCATTAgtagtcgtaggccttaaccatttgtgcctgat                                                                                      >  2:481148‑M1/1‑219 (MQ=255)
                                                                                                                                                        ttCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatg                                                                                     >  1:207794‑M1/1‑218 (MQ=255)
                                                                                                                                                             aCTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacg                                                                                >  1:411629‑M1/1‑213 (MQ=255)
                                                                                                                                                                      gttgggttggGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcg                                                                       >  1:243618‑M1/1‑204 (MQ=255)
                                                                                                                                                                         gggttgggtGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtct                                                                    <  1:806393‑M1/251‑51 (MQ=255)
                                                                                                                                                                         gggttgggtGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtc                                                                     <  1:825058‑M1/250‑50 (MQ=255)
                                                                                                                                                                             ggggggtgaatagAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGGTCCGTACCGTTAAAACGGGACGCTATTCAGACAATGCTTTTTTAGGCAGCTAAGTTGGCAGGGGATCGTGTTGGCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTTGCTGAGGCATTTgtagtcgtaggcattagacatttgggcctgatgcgacggttgacgcgtcttat                                                                 <  2:255030‑M1/245‑54 (MQ=255)
                                                                                                                                                                                gtgtGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTCTTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcggaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcagg                                                             >  2:854079‑M1/1‑194 (MQ=255)
                                                                                                                                                                                        ataAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATGTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgacgcgacgcttgacgcggcttatcaggcctacaac                                                     <  1:704872‑M1/251‑66 (MQ=255)
                                                                                                                                                                                                        cacGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggc                                     >  1:640448‑M1/1‑170 (MQ=255)
                                                                                                                                                                                                          cGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaa                                   <  1:321388‑M1/251‑84 (MQ=255)
                                                                                                                                                                                                           ggAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaat                                  >  1:732416‑M1/1‑167 (MQ=255)
                                                                                                                                                                                                             aaGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattg                                >  2:58846‑M1/1‑165 (MQ=255)
                                                                                                                                                                                                              aGCCGGTTTACGCTGGATGGATCCGGACCGTTAAAACGTGACGCTAGTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAAGTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACGGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattgg                               <  2:315316‑M1/251‑88 (MQ=255)
                                                                                                                                                                                                                ccGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtg                             <  1:516216‑M1/251‑90 (MQ=255)
                                                                                                                                                                                                                  ggTTTACGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgc                            >  2:660479‑M1/1‑160 (MQ=255)
                                                                                                                                                                                                                       aCGCTGGATGGATCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatg                      <  2:806290‑M1/251‑97 (MQ=255)
                                                                                                                                                                                                                            ggatggatCCGTACCGTTAAATCGTTACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTTGCAGGGGATCGTGTTGTCCGAAAACAGGAAAAGCGATATTTAAAAGAGGTTAATTATGAAACTTGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgttcctgatgcgacgcttgccgcgtcttatcaggcctacaaccggggccgcatccgccaattggtgcacaatgcctga                 <  2:619364‑M1/251‑102 (MQ=255)
                                                                                                                                                                                                                                ggatCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcg             >  2:483033‑M1/1‑146 (MQ=255)
                                                                                                                                                                                                                                  atCCGTACCGTTAAAACGTGACGGTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtaggcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat           <  2:447570‑M1/251‑108 (MQ=255)
                                                                                                                                                                                                                                  atCCGTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat           >  2:549412‑M1/1‑144 (MQ=255)
                                                                                                                                                                                                                                     cGTACCGTTAAAACGTGACGCTATTAAGACCATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGCGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgccggatgcgacgcttgacgcgtcttatcaggcctacaaccggggccgcatccggcaattggtgcacaatgcctgatgcgatgct        >  2:678681‑M1/1‑141 (MQ=255)
                                                                                                                                                                                                                                      gTACCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt       <  1:824599‑M1/251‑112 (MQ=255)
                                                                                                                                                                                                                                        aCCGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttg      <  1:91652‑M1/250‑113 (MQ=255)
                                                                                                                                                                                                                                          cGTTAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacg   <  1:832822‑M1/251‑116 (MQ=255)
                                                                                                                                                                                                                                            ttAAAACGTGACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCCGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTAgtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgc  >  1:111627‑M1/1‑134 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
ATCATATGCAGCGGATGGAAGGTGAATTAACCAACACCATGACTCAGTGGGCTGATAAAATTGGTCACCTGCAAATTGCCGATAATCCGCATCGCGGCGAACCGGGAACCGGAGAAATTAATTATGATTATCTCTTTAAGGTAATCGAAAATTCTGACTACAACGGTTGGGTTGGGTGTGAATATAAACCCCAAACCACCACGGAAGCCGGTTTACGCTGGATGGATCCGTACCGTTAAAACGTAACGCTATTCAGACAATGCTTTTTTAGGCCGCTAAGTTGGCAGGGGATCGTGTTGTCTGAATTCAGGAAAAGCGAAATTTAAAAGAGGTTAATTATGAAACTGGGATTTATTGGCTTAGGCATTATGGGTACACCGATGGCCATTAATCTGGCGCGTGCCGGTCATCAATTACATGTCACGACCATTGGACCGGTTGCTGATGAATTACTGTCACTGGGTGCCGTCAGTGTTGAAACTGCT  >  NC_000913/536248‑536732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: