Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,367,944 G→A A40T (GCG→ACG)  yhcG → DUF1016 domain‑containing protein YhcG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,367,9440GA94.6% 106.7 / ‑1.4 37A40T (GCG→ACG) yhcGDUF1016 domain‑containing protein YhcG
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (17/18);  total (18/19)
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.56e-01

TTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTC  >  NC_000913/3367699‑3368185
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
ttCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCacgacg                                                                                                                                                                                                                                              >  2:445943/1‑251 (MQ=255)
           tctTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAAt                                                                                                                                                                                                                                   >  1:23283/1‑251 (MQ=255)
                  atatCTTGGGGAAAGGAAATGTAATACGCTGTTCTTAATTGCTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCTACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGCAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCCTCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTTACGCGTTAATGACCacgacggaaccggaacatggccg                                                                                                                                                                                                                             >  2:212628/1‑233 (MQ=255)
                         gggtaaaGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAACCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTCCCAGGAAATTGGCCGACGAAtt                                                                                                                                                                                                                      <  1:45372/246‑1 (MQ=255)
                             aaaGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGAATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGAAGACAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCg                                                                                                                                                                                                                  >  2:461548/1‑250 (MQ=255)
                              aaGTAAATGTAATCCGTGGTTCTTTATTGTTCAGGGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACTCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGaa                                                                                                                                                                                                                <  2:340106/251‑1 (MQ=255)
                                      gTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCAGATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAAc                                                                                                                                                                                                         <  1:434227/250‑1 (MQ=255)
                                         aTACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGCCGAATTTGAACAAgg                                                                                                                                                                                                     <  2:111334/251‑1 (MQ=255)
                                                gTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCg                                                                                                                                                                                                >  1:215895/1‑249 (MQ=255)
                                                 ttCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTAAACAAGGTGGCGAgg                                                                                                                                                                                             >  1:30908/1‑251 (MQ=255)
                                                              agttcTCACCCCCCAGATTATAGCTCATTCCATTATTCAACGTAGCCCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGTAACCACAGCAGGCGACCAGCAAATCCACTACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGTCCAGGGCTGCGTa                                                                                                                                                                                <  2:460006/247‑1 (MQ=255)
                                                                      ccccccAGATTATAGCTCCTGCCATTATTCAACGTAGCTCCATATAAAGGGACAGATAATCGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGTCCAGTGCTGCGTATGGTgc                                                                                                                                                                          <  2:105366/249‑1 (MQ=255)
                                                                                tataGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCa                                                                                                                                                               >  1:52635/1‑250 (MQ=255)
                                                                                 ataGCTCATCCCATTATTCAACTTAGCGCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAAGCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTTCGTATGGTGCGCAGCTAATCAAg                                                                                                                                                             <  2:25627/251‑1 (MQ=255)
                                                                                      tCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGAAGCTACTCAAGCGACt                                                                                                                                                        >  1:446593/1‑251 (MQ=255)
                                                                                          tCCATTATTCAACTAAGCACCATAATAAGGGACAGATAATGGAATCGCTCTCAGAATGAACCACAGCAGGCTACCAGCAAATCCAAGACGGTATTATTCATCTGGTTGATAGCGCCCGGACGGAAACGGGACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATc                                                                                                                                                     <  2:240014/250‑1 (MQ=255)
                                                                                                attCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAt                                                                                                                                              >  2:257900/1‑251 (MQ=255)
                                                                                                 ttCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGAtt                                                                                                                                             >  1:365750/1‑251 (MQ=255)
                                                                                                            tCCATATAAAGGGACAGATAATGGAATAGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGAGCCCGGACGGACACGGTACGTAGCGATACCGCGTTAATGACCACGACGCACCAGGAAATTGGCCGACGAATTGTCGAATTTGAAAAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCATCTAATAAAGCGACTATCAAGGTATTTATGTCTAAgg                                                                                                                                  >  2:186201/1‑251 (MQ=255)
                                                                                                                             aTAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTtctc                                                                                                                 <  2:41596/251‑1 (MQ=255)
                                                                                                                               aaTGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTg                                                                                                               <  1:257900/251‑1 (MQ=255)
                                                                                                                                            tcAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGAATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCCGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGGATAAGCGTGGGTTCTCTGCAAAAAACTTACg                                                                                                  >  1:154281/1‑251 (MQ=255)
                                                                                                                                                aaGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTTGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCa                                                                                               >  1:74894/1‑250 (MQ=255)
                                                                                                                                                          agcagGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCttt                                                                                    >  2:81446/1‑251 (MQ=255)
                                                                                                                                                           gcagGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCtttt                                                                                   <  1:249271/251‑1 (MQ=255)
                                                                                                                                                           gcagGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAAAGAGGCtttt                                                                                   >  2:249271/1‑251 (MQ=255)
                                                                                                                                                             agGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCtttttt                                                                                 <  1:79361/251‑1 (MQ=255)
                                                                                                                                                                   ccAGCAAATCCACGATGGTATTATTCATCTGGTCTATAGCTCCCGGACTGAAACGGTACGTAGCGTTAACGCGTTAATGCCCACGACTTACCAGGAAATTGTCCGACGAATTGGCGAAGTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGTTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCATATGATGCTTTTTTACctct                                                                           <  2:349699/251‑1 (MQ=255)
                                                                                                                                                                     aGCAAATCCACTACAGTATTATTCATCTGGTCGATAGCGCTCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGTCCACGACGTACCAGGAAATTTGCCGACGAATTGTCGAAGTTGAACAAGGGGGCGAGGCCAGGGCTGTTTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAATGTATAAGCGTGGTTTCTCTGCAAAAAACTTACCCCAAATGAGGCTTTTTTACCTCttt                                                                         <  2:85803/251‑1 (MQ=255)
                                                                                                                                                                             cACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACTGAAACGGTACGGTGCGTTAACGCGTTAATTACCACGACTTGCCACGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGTTAATCATGCGACTATCAAAGTATTTATGTCTAAGGTATAAGTGTGGTTTATCTGCAAAAAACTTACGCCAAATGAGGCTTTTATGCCGCTTTTTTctata                                                                 <  2:159799/251‑5 (MQ=255)
                                                                                                                                                                             cACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAAc                                                                  >  1:410514/1‑250 (MQ=255)
                                                                                                                                                                                gacgGAATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGAACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGTCTTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTCTGTCTAGGGAAAAAGCGTGGGTTCTCTGCAAACAACTTCCGACACATGAGGCTTTTTTACCTCTTTTTTCAAAATGt                                                              >  2:60438/1‑251 (MQ=255)
                                                                                                                                                                                           tCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCAc                                                    <  1:81446/250‑1 (MQ=255)
                                                                                                                                                                                                                gacggaAACGGTACGTATCGTTCACGCGTTAATGACCACGACGTACCAGGGAATTGGCCGACGAATTGGCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTTATCAAGCGACTATAAAAGGATTTATGTCTAAGGTATAAGAGTGGTTTCTCTGCAAAAAACTTACGCCTAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTaa                              <  2:392879/251‑1 (MQ=255)
                                                                                                                                                                                                                            aCGTAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGa                   >  1:256079/1‑250 (MQ=255)
                                                                                                                                                                                                                               tAGCGTTAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTc               >  1:10752/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tAACGCGTTAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAAc                                                                                                       <  1:231508/157‑1 (MQ=255)
                                                                                                                                                                                                                                     tAACGCGTTAATGACCAAGACGTACCAGGCAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGGGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTATCTGCAAAAAAc                                                                                                       >  2:231508/1‑157 (MQ=255)
                                                                                                                                                                                                                                             tAATGACCACGACGTACCAGGAAATTGGCCGACGAATTGTCGAAATTGAACAAGGTGTCGAGGCCAGGGTTGCGTATGGTGCGCAGCTAATCAAGCAACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTc  <  2:173831/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
TTCAGGAAAGCTCTTCGCATATCTTGGGGAAAGTAAATGTAATACGCTGTTCTTAATTGTTCAGTGCTCACCCCCCAGAGTATAGCTCATTCCATTATTCAACTTAGCTCCATATAAAGGGACAGATAATGGAATCGCTCTCAGAAGGAACCACAGCAGGCTACCAGCAAATCCACGACGGTATTATTCATCTGGTCGATAGCGCCCGGACGGAAACGGTACGTAGCGTTAACGCGTTAATGACCGCGACGTACCAGGAAATTGGCCGACGAATTGTCGAATTTGAACAAGGTGGCGAGGCCAGGGCTGCGTATGGTGCGCAGCTAATCAAGCGACTATCAAAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTC  >  NC_000913/3367699‑3368185

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