Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,084,281 T→C N68D (AAT→GAT)  unknown ← N‑acetylmuramoyl‑L‑alanine amidase (EC 3.5.1.28)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,084,2810TC98.5% 179.1 / NA 68N68D (AAT→GAT) unknownN‑acetylmuramoyl‑L‑alanine amidase (EC 3.5.1.28)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/36);  minor base A (0/1);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATTTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCAATAATCTGT  >  1/1084042‑1084521
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
gCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCtt                                                                                                                                                                                                                                       >  1:192374/1‑251 (MQ=255)
   cAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTa                                                                                                                                                                                                                                    <  1:835339/251‑1 (MQ=255)
    aaCAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTaa                                                                                                                                                                                                                                   <  1:500321/251‑1 (MQ=255)
                       cATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAgg                                                                                                                                                                                                                <  1:245241/251‑1 (MQ=255)
                            tCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTa                                                                                                                                                                                                            >  1:679812/1‑250 (MQ=255)
                             cTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACa                                                                                                                                                                                                          >  1:174828/1‑251 (MQ=255)
                             cTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACa                                                                                                                                                                                                          >  1:517433/1‑251 (MQ=255)
                              ttCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACAt                                                                                                                                                                                                         <  1:178048/251‑1 (MQ=255)
                              ttCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACAt                                                                                                                                                                                                         <  1:463528/251‑1 (MQ=255)
                                                    tCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGt                                                                                                                                                                                    <  1:559203/250‑1 (MQ=255)
                                                                   tttGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAgcgc                                                                                                                                                                     <  1:302204/250‑1 (MQ=255)
                                                                      gACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCggg                                                                                                                                                                 >  1:231651/1‑251 (MQ=255)
                                                                             cccATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGtt                                                                                                                                                                                                             <  1:550068/200‑1 (MQ=255)
                                                                             cccATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCttt                                                                                                                                                           >  1:607068/1‑250 (MQ=255)
                                                                                aTACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGaaa                                                                                                                                                                                                                  <  1:73/192‑1 (MQ=255)
                                                                                aTACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGa                                                                                                                                                                                                                    <  1:82982/190‑1 (MQ=255)
                                                                                aTACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCg                                                                                                                                                       >  1:897432/1‑251 (MQ=255)
                                                                                              cgTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGAt                                                                                                                                          <  1:160546/250‑1 (MQ=255)
                                                                                              cgTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATc                                                                                                                                         <  1:848267/251‑1 (MQ=255)
                                                                                                 aaaTTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATGGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAAc                                                                                                                                      <  1:515537/251‑1 (MQ=255)
                                                                                                    ttCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATCCGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGc                                                                                                                                   <  1:782953/251‑1 (MQ=255)
                                                                                                      cGGGTTTTCTGAGAAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTCCGCGCAATTTCGGACCGTTGCGGCACGAAAATAAAGTAATCGCCGTTGCGGGTAAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGTTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCaa                                                                                                                                 <  1:23404/251‑1 (MQ=255)
                                                                                                        gggTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAggg                                                                                                                                                                                                               <  1:847361/171‑1 (MQ=255)
                                                                                                          gTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGt                                                                                                                             >  1:632981/1‑251 (MQ=255)
                                                                                                             tCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTcc                                                                                                                           >  1:211650/1‑250 (MQ=255)
                                                                                                                     aTCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTAc                                                                                                                  <  1:147716/251‑1 (MQ=255)
                                                                                                                              aaTGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCgggg                                                                                                                                                                >  1:466052/1‑196 (MQ=255)
                                                                                                                                      ccAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGTAATAAAGTAATCGCCGTTGCGGGTTAATACAGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCCCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCa                                                                                                   >  1:561739/1‑249 (MQ=255)
                                                                                                                                                      gtgCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAgg                                                                                                                                                                                                                >  1:521719/1‑124 (MQ=255)
                                                                                                                                                             cgcgAAATTTCGGAACGTTGCGGCACGAAAATAAAGTAATCGCCGTTGCGGGTTAATCCGGCATGAAAACGCGGATCCTTATCAAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGAc                                                                          <  1:573731/251‑1 (MQ=255)
                                                                                                                                                                cggggagatgtgaaaaagagacaggAAAAAAAGTAATCGCCGTGGGGGGTTAATACGGCATGAAACCGCGGATCTTTATATAAAAGCGTTTTTAATTCCGGGGGGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCGGATGGCGCCGGGGTTTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAg                                                                       <  1:894159/228‑1 (MQ=255)
                                                                                                                                                                   aTTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGaaa                                                                    >  1:900472/1‑251 (MQ=255)
                                                                                                                                                                    tttCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAaat                                                                   >  1:91357/1‑251 (MQ=255)
                                                                                                                                                                    tttCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAACTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAaat                                                                   >  1:313331/1‑251 (MQ=255)
                                                                                                                                                                         gAACGTTGCGGCACGAAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGTTCCTTATCTAAAAGCGCTTTTAATTCCCGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAAtt                                                               <  1:707271/250‑1 (MQ=255)
                                                                                                                                                                                    aaagaaaaTAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTGGCGCCCGATGGCGCCGGGGTTTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGc                                                   <  1:534376/246‑1 (MQ=255)
                                                                                                                                                                                      cGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGa                                                         <  1:676312/243‑1 (MQ=255)
                                                                                                                                                                                        gAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTc                                                                                                                                                              >  1:236998/1‑140 (MQ=255)
                                                                                                                                                                                        gAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCa                                                                                                                                                    >  1:626060/1‑150 (MQ=255)
                                                                                                                                                                                        gAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAg                                               >  1:817228/1‑251 (MQ=255)
                                                                                                                                                                                            aaaaGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAGGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCGGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTc                                           <  1:829495/250‑1 (MQ=255)
                                                                                                                                                                                                gTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAAttt                                                              <  1:162973/228‑1 (MQ=255)
                                                                                                                                                                                                gTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                                       >  1:832392/1‑251 (MQ=255)
                                                                                                                                                                                                gTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATACGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                                       >  1:441856/1‑251 (MQ=255)
                                                                                                                                                                                                 tAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCCGGGGCGATAGAAAGGTTTACTTTTTTTTCATAATTTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                                       <  1:245403/250‑1 (MQ=255)
                                                                                                                                                                                                 tAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                                       <  1:198748/250‑1 (MQ=255)
                                                                                                                                                                                                 tAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCGGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                                       <  1:662164/250‑1 (MQ=255)
                                                                                                                                                                                                   aTCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCTTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAGGGGTCACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTgg                                    <  1:424880/251‑1 (MQ=255)
                                                                                                                                                                                                       ccGTTGCGGGTTAATCCGGCATGAAACCGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAAAAAGGTTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTAc                                 <  1:642493/250‑1 (MQ=255)
                                                                                                                                                                                                       ccGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAg                                                                       >  1:724244/1‑212 (MQ=255)
                                                                                                                                                                                                        cGTTGCGGGTTAATACGGCATAAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTAcat                               <  1:796432/251‑1 (MQ=255)
                                                                                                                                                                                                           tGCGGTTTAATACGGCATGAAACCGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTAcatcat                            <  1:743525/251‑1 (MQ=255)
                                                                                                                                                                                                              gggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGc                                                                                                                                                                                                                          >  1:206137/1‑58 (MQ=255)
                                                                                                                                                                                                              gggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTAc                                                                                                                  >  1:401383/1‑162 (MQ=255)
                                                                                                                                                                                                               ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAg                                                                                                                                                                                      <  1:567668/93‑1 (MQ=255)
                                                                                                                                                                                                               ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAgcgc                                                                                                                                                                     <  1:613256/110‑1 (MQ=255)
                                                                                                                                                                                                               ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCg                                                                                                                      >  1:362065/1‑157 (MQ=255)
                                                                                                                                                                                                               ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAg                                                          >  1:417177/1‑217 (MQ=255)
                                                                                                                                                                                                               ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGtt                        >  1:369307/1‑251 (MQ=255)
                                                                                                                                                                                                               ggTTAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGtt                        >  1:137743/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttAATACGGCATGAAAACGCGGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTAt                      <  1:197249/251‑1 (MQ=255)
                                                                                                                                                                                                                      aCGGCATGAAAACGCGGTTCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAATTTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAAt                 <  1:727814/251‑1 (MQ=255)
                                                                                                                                                                                                                         gCATGAAAACGCGGATCCTTATCTAAAAGCGTTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCAAAAATTTTTGGGTGGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATtg                                       <  1:374669/226‑1 (MQ=255)
                                                                                                                                                                                                                                  cgcgGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCaataat      >  1:742953/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cgcgGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCaataat      >  1:725436/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cgcgGATCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCAATAAtc     <  1:192068/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCAATAAtctgt  >  1:250436/1‑248 (MQ=255)
                                                                                                                                                                                                                                        tCCTTATCTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCAATAAtctgt  >  1:742078/1‑248 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
GCCCAACAGCTCCGAACGTTTTTCATGATCTTCCAGCCATTGCCCCATTTCGTCATTGGCGCGACGGTTTGACAACACCCATACGGAAGCGCCGCGTAAATTCGGGGTTTCTGAGGAATCTGCATGAATGGAAACCAAATAATTGGCTTTGTGCTTACGCGCAATTTCGGAACGTTGCGGCACGGAAATAAAGTAATCGCCGTTGCGGGTTAATACGGCATGAAAACGCGGATCCTTATTTAAAAGCGCTTTTAATTCACGGGCGATAGAAAGGGTTACATTTTTTTCATAAATTTTTAGGTTGCGCCCGATAGCGCCGGGGTCTTTCCCGCCATGACCGGGATCAACGGCAATTGTCCATGCGTTACCGGCGAGAACGGGCATTGCAAAGAAAAAAACGCTGTTGACAAGAAATAATTTAAAGATATTGCCGAGTTTCATTGTGGTACATCATCGTTATTTAATTGTTCAATAATCTGT  >  1/1084042‑1084521

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