Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,856,968 T→A I268F (ATT→TTT)  ydjI ← L‑glycero‑L‑galacto‑octuluronate‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,856,9680TA95.1% 129.5 / ‑0.1 41I268F (ATT→TTT) ydjIL‑glycero‑L‑galacto‑octuluronate‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (20/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GCGGCATCTTTACCAATACGACTCATTAATGCTGTGCGATGGCCCAGACGAGAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAATTTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTGCAAAACCAGCGGAACAGAAACGGCATCGCGTACGCGTTTGACGACCTCAAAGTTC  >  NC_000913/1856723‑1857209
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
gCGGCATCTTTACCAATACGACTCATTAATGCTGTGCGATGGCCCAGACGAGAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAattt                                                                                                                                                                                                                                                <  1:461721/249‑1 (MQ=255)
         ttACCAATACGACTCATTAATGCTGTGCGATGGCCCAGACGAGAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTCTTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGccc                                                                                                                                                                                                                                     >  2:107454/1‑251 (MQ=255)
           aCCAATACGACTCATTAATGCTGTGCGATGGCCCAGACGAGAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGt                                                                                                                                                                                                                                   >  1:537589/1‑251 (MQ=255)
                                  tgCGATGGCCCAGACGAGAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTc                                                                                                                                                                                                            >  2:476379/1‑251 (MQ=255)
                                    caagggCCCAGACGATAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGTAATGGAATATCAATAATAGCGGCACCTAGACAAATAACGTCGAGATTAGCCATATCATTAGTCCGCTTTGCCATCAGAAACGAACAGTTTAAATTTTTCCAATGCCCGTTCGTTTACAGCTTTACGTACTTCAc                                                                                                                                                                                                          <  2:166824/247‑1 (MQ=255)
                                                      ataatTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAAGCCACAGCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCGTTACGTACTTCACGTTCTAAATGCAGGAAAg                                                                                                                                                                                        <  2:486421/251‑1 (MQ=255)
                                                         atTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATATCGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCt                                                                                                                                                                                     >  2:139659/1‑251 (MQ=255)
                                                               gcgtcagttATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATcc                                                                                                                                                                                <  2:21403/243‑1 (MQ=255)
                                                                       tATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTtttc                                                                                                                                                                        <  2:419275/250‑1 (MQ=255)
                                                                           gCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGACAGGCTGATCCTGGTtttctttcc                                                                                                                                                                   >  2:299006/1‑250 (MQ=255)
                                                                                ccccccGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGTCTCCTATCCAAATAACGTCGAGATTATCCATATTATGATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAAGGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGGTGTCTTTCACTGCa                                                                                                                                                              <  1:168013/250‑1 (MQ=255)
                                                                                 cccccGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTGCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTg                                                                                                                                                                <  1:442546/248‑1 (MQ=255)
                                                                                  ccccGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAAc                                                                                                                                                            <  1:575962/251‑1 (MQ=255)
                                                                                      ggtggtCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGCCCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGAATATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATTCATTAAAGGCTGCTCCTGGTTTTCTTTCACTGCAACCAtt                                                                                                                                                        <  1:314332/251‑1 (MQ=255)
                                                                                       gtggtCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGACCCTGGTTTTCTTTCACTGCAACCAt                                                                                                                                                         >  2:386355/1‑249 (MQ=255)
                                                                                           tCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGc                                                                                                                                                   >  2:261050/1‑251 (MQ=255)
                                                                                                 cGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACa                                                                                                                                             <  1:215838/251‑1 (MQ=255)
                                                                                                 cGATTCTTTCAAGAGGATAAGAAGCCACATCAAAGATATTTTTACTGACCTGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACa                                                                                                                                             <  2:328031/251‑1 (MQ=255)
                                                                                                     tCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCTGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCAATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAATTTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGc                                                                                                                                         >  2:212959/1‑251 (MQ=255)
                                                                                                       tttCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTc                                                                                                                                       >  2:444058/1‑251 (MQ=255)
                                                                                                                gATAAGACTCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGTTCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAAt                                                                                                                              <  1:443728/251‑1 (MQ=255)
                                                                                                                    aGAATCCACATCAAAGATATTTTTACTGTCCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCACTATTCCGCTTTGCCATCAGAACCGAACATTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACTTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTTACTGCAACCAGTGCTGCCTGACATAGCTCCGTATGAATGTTg                                                                                                                          <  1:217790/251‑1 (MQ=255)
                                                                                                                                aaaGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTg                                                                                                                                                       >  2:304491/1‑210 (MQ=255)
                                                                                                                                aaaGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTg                                                                                                                                                       <  1:304491/210‑1 (MQ=255)
                                                                                                                                   gATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGAGCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATTAATGTTGATTTTCGCGATAccc                                                                                                           <  2:44819/251‑1 (MQ=255)
                                                                                                                                                aCCGGTGGCAATGGAATATCAACAATAGCGGCTCCGATACAAATAACGTCGATATTATCCATATCATTATTCCGCTTTGCCATCGGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCAGTGCGGCCTGACTGAGCTCCGTATGAATGTTGATTGTCGCGATACCCATCGAAATTgtag                                                                                              <  2:497181/251‑4 (MQ=255)
                                                                                                                                                 ccGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGt                                                                                             <  2:486010/251‑1 (MQ=255)
                                                                                                                                                  cGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGtt                                                                                            <  1:107454/251‑1 (MQ=255)
                                                                                                                                                    gTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGtttt                                                                                          <  1:332088/251‑1 (MQ=255)
                                                                                                                                                           tGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGCCAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGt                                                                                     >  2:3832/1‑249 (MQ=255)
                                                                                                                                                            ggAATATCACCAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTCCGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCgg                                                                                  >  2:548516/1‑251 (MQ=255)
                                                                                                                                                                        aTAGCGGCACCTATACAAATAACGTCGAGATTATCAATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAATTTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATcc                                                                      <  1:212959/251‑1 (MQ=255)
                                                                                                                                                                                     taCAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTg                                                         >  1:394194/1‑251 (MQ=255)
                                                                                                                                                                                               gTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGCACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATTCACCGTGCAAAACCAGc                                               >  2:131625/1‑251 (MQ=255)
                                                                                                                                                                                                       tATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGCTTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTGCAAAACCAGCGGAAc                                          >  2:564064/1‑248 (MQ=255)
                                                                                                                                                                                                                  ttattCAGCTTTGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCCGGAAAGGCTGATCCTGGTTTTCTTTCACTGCACCCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAAAGAAATTGCAGTTTTACTGTCGGCATCACTGATCCACGATGCACAGTGCCAAACCAGCGGAACAGAAACGGCATcgc                            >  2:436370/1‑251 (MQ=255)
                                                                                                                                                                                                                           tttGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTGCAAAACCAGCGGAACAGAAACGGCATCGCGTACGCGtt                   >  1:459786/1‑251 (MQ=255)
                                                                                                                                                                                                                           tttGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATAGCGGCCTGACAGAGCTCCGTATGACTGTTGATTTTCGCGATACCCAACCAAATTGCAGTTTTAATGTCGGCATCACTGATCCCAGATGCACCGAGTAAAACCAGCGGAAAAGAAACGACATCGCGTACGCGtt                   >  1:75247/1‑251 (MQ=255)
                                                                                                                                                                                                                             tGCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTGCAAAACCAGCGGAACAGAAACGGCATCGCGTACGCGttt                  <  1:397364/250‑1 (MQ=255)
                                                                                                                                                                                                                              gCCATCAGAACCGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTTCTCCCCGATGCACCGTGCACAACCAGCGGAACAGAAACGGCATCGCGTACGCTTTTg                 >  2:45687/1‑250 (MQ=255)
                                                                                                                                                                                                                                        ccGAACAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTGCAAAACCAGCGGAACAGAAACGGCATCGCGTACGCGTTTGACGACCTCaa       >  1:151058/1‑250 (MQ=255)
                                                                                                                                                                                                                                             cAGTTTAAATTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCCCTGTTCCCCGATGCACCGTGCAAAACCAGCGGAACAGAAACGGCATCGCGTACGCGTTTGACGACCTCAAAGTTc  >  2:104371/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
GCGGCATCTTTACCAATACGACTCATTAATGCTGTGCGATGGCCCAGACGAGAAATAATTGTTGCTTCATTTATTGCATCACCCCCGGTGGTCATTGCGATTCTTTCAAGAGGATAAGAATCCACATCAAAGATATTTTTACTGACCGGTTGCAATGGAATATCAACAATAGCGGCACCTATACAAATAACGTCGAGATTATCCATATCATTATTCCGCTTTGCCATCAGAACCGAACAGTTTAATTTTTTCCAATGCCCGTTCTTTTACAGCTTTACGTACTTCACGTTCTAAATGCAGGAAAGGCTGATCCTGGTTTTCTTTCACTGCAACCATTGCGGCCTGACAGAGCTCCGTATGAATGTTGATTTTCGCGATACCCAACGAAATTGCAGTTTTAATGTCGGCATCACTGATCCCCGATGCACCGTGCAAAACCAGCGGAACAGAAACGGCATCGCGTACGCGTTTGACGACCTCAAAGTTC  >  NC_000913/1856723‑1857209

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