Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,623,743 C→T V125V (GTG→GTA)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,623,7430CT84.1% 98.7 / 7.2 44V125V (GTG→GTA) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base T (18/19);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACC  >  NC_000913/4623499‑4623989
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
ccTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCac‑tacc                                                                                                                                                                                                                                                      >  1:34242/1‑248 (MQ=255)
       tatagCCTCGCAAACCTGCTCATGGGTGAACCCCGACAACAGCTCGGTGAGGTTGGTCACGCGGGAAAGTTCCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGCTAGTTTGCCAGGCGGCTGAGGTCGGCATTGGGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGAGCTAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCTC‑CACCAGATCGtgac                                                                                                                                                                                                                                            >  2:86431/5‑248 (MQ=255)
        aTGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACg                                                                                                                                                                                                                                           >  1:263535/1‑251 (MQ=255)
           gCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTcc                                                                                                                                                                                                                                        >  1:270000/1‑251 (MQ=255)
            catCGCAAATCTGCTCATGGGTGATCCCCGGCAAAAGCCCGGTGCGGTTCGTCACGCGGGACCGTACCGACGTAATGCCTTTAGCCGCCCGTTTCTTTTTATCCGGCTTGAGCTCGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGCGCCGTGGTGGAAGCAGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGCCTTTGCGGTCGCCTTCGACGGTTTTCACTCAC‑AGATCGTTACGTCCg                                                                                                                                                                                                                                       >  2:267651/3‑251 (MQ=255)
                   aaCCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCtt                                                                                                                                                                                                                                >  1:102238/1‑251 (MQ=255)
                     ccTGCTCATGGGTGATCCCCGGCAACAGCGAGGCGCGGTTGGTCACGCGGGAACGTAACGACGTAATGCCTATCGCCGCCAGTTTCTTTTTATCCGGGATGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTTAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGTGTTCGCGATAGGCCGAGCCCGAGACTTTGCGGTCGACTTCGACGGTTTTCAC‑TACCCGATcagtgcgtccgtacgcctcg                                                                                                                                                                                                                              >  2:107148/1‑232 (MQ=255)
                            aTGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTg                                                                                                                                                                                                                        >  1:170614/1‑250 (MQ=255)
                              gggTGATCCCCGTCATCAGTGCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTGTGCCAGGCGGCGGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCTATAGGCCGAGCCTGAGACTTGGCGGTCGCCTTCGACGTGTTTCAC‑TACCAGATCGTTACGTCCTGACGCTTCGGCTCTGAc                                                                                                                                                                                                                      <  2:197114/250‑1 (MQ=255)
                                 tGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCg                                                                                                                                                                                                                  >  1:335582/1‑251 (MQ=255)
                                            aaCAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGACGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTccg                                                                                                                                                                                                       >  2:143106/1‑249 (MQ=255)
                                                    ggTGAGGGTGGTCACGCGGGAATGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTGTTAACCGGATTGAGATAGGTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCGGTGGTGGAAGCCGCGATCTTTGGTGTCGCGATAGGACGAGTCTGAGACTTTGCGGTCGCCTTCGAAGGTTTTCAC‑TACCAGATCGTTAAGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTGCAgcgaaggga                                                                                                                                                                                               <  2:164450/251‑7 (MQ=255)
                                                            tGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCAGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGCGCACAATc                                                                                                                                                                                      >  1:379605/1‑251 (MQ=255)
                                                             ggTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGGCGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCATAATCg                                                                                                                                                                                     <  1:86638/251‑1 (MQ=255)
                                                                cacccgttaaCGTACCGACGTAATGCCGTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATATTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCTGTCGCCTGCGACGGTTTTCAC‑TACCAGAGCGTTACGTCCGTACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACg                                                                                                                                                                                  <  2:360017/243‑1 (MQ=255)
                                                                  cgcgGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTg                                                                                                                                                                                >  1:346641/1‑251 (MQ=255)
                                                                    cgGGAACGTACCGACGTAATGCCTGTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGga                                                                                                                                                                              <  1:82821/251‑1 (MQ=255)
                                                                           gTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGAGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACggggggagagtt                                                                                                                                                                       >  2:98038/1‑240 (MQ=255)
                                                                             tccGACGTAATGCCGTTCGCCGCCAGGTTCTGTTTATCCGGAGTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCAGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTGCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGAAT‑TGAGCACAATCGACGTGGAGATAGttt                                                                                                                                                                      <  2:308151/249‑1 (MQ=255)
                                                                                 aCGTAATGCCTTTCGCCGCCAGTTTCTTTTTAGCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCg                                                                                                                                                                 <  2:159357/251‑1 (MQ=255)
                                                                                        gCCTTTCGCCGCCAGTTTCTTGTTATCCGGATTGAGAGAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGGAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGTTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCg                                                                                                                                                          <  2:275043/251‑1 (MQ=255)
                                                                                                                         gagaTAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGc                                                                                                                         >  1:96919/1‑251 (MQ=255)
                                                                                                                          agttAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGTCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCCGCCATAAACGAAAAGCAGGTATTGcc                                                                                                                        >  2:2124/4‑251 (MQ=255)
                                                                                                                               ttttGCCAGGCGGCTGTGGTCGGCATTGAGTAGAACGGTGACGTGGTGGAAGCCGCGCTCTTTGGTTTCGCGATAGACCGAGCCTGAGACTTTGCGGACGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCCTCGGCGCTGACGCCGAGCGCGTTCAGCGACT‑TGCGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTACAGAAGGAATTGCCgagc                                                                                                                    >  2:243229/2‑249 (MQ=255)
                                                                                                                                 ttGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATc                                                                                                                  <  1:8499/250‑1 (MQ=255)
                                                                                                                                                 gTCGGCATTGAGTAGCAAGGGGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATCAAGGTAAAGCAGGTATTGCCGAGATCGTGGACCACCGcgccgc                                                                                                 >  2:164283/1‑251 (MQ=255)
                                                                                                                                                    ggCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCAc                                                                                              >  1:200217/1‑251 (MQ=255)
                                                                                                                                                     ggaTTGAGTAGCAGGGTGCCGGGGTGGAAGCCGCGATCTTTGGTTTCGCTATAGGCCGAGCCTGAGTCTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGATCACAATCGACGTGGAGATAGTTTTATCTTGCTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTTGAACACCGCGCCGCCAcc                                                                                             <  2:286890/249‑1 (MQ=255)
                                                                                                                                                               gCAAGGTGCCGGGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TTCCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTgcgc                                                                                   <  1:243229/251‑1 (MQ=255)
                                                                                                                                                                 aaGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TTCCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGAAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTgcgccg                                                                                 >  2:5048/1‑251 (MQ=255)
                                                                                                                                                                 aaGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTgcgccg                                                                                 >  2:89367/1‑251 (MQ=255)
                                                                                                                                                                  aGGTGCAGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGCGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGAATTGCCGAGATCGGGGCACACCGCGCCGCCACCGCTACTgcgccgc                                                                                >  2:154191/1‑251 (MQ=255)
                                                                                                                                                                                  gCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTAt                                                                <  1:164283/251‑1 (MQ=255)
                                                                                                                                                                                    cgcgATCTTTGGTTTCGCTATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGGATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTAtct                                                              <  2:335582/251‑1 (MQ=255)
                                                                                                                                                                                            tgggTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCCC‑TACCAGATCGTTTCGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCGGGCGGAc                                                                     <  1:370503/234‑1 (MQ=255)
                                                                                                                                                                                               gTTTCGCGATAGGCCGAGCCTGAGACTTTGCGTGGGCCGTCGACGGTGTGCAC‑TACCAGATCGTTACGTCCGGACGCTGCGGCGCGGACGCCGAGCGCGTTCCGCGCATGTG‑GCACAATCGACGTGTAGATAGTTTTATCGTACGCCGGTTTGCCAGCGGTAAGGGTAAAGCAGGGATTGCCGAGTTCGTGGAACGCCGCGCCGCCGCCGTTACGGCGCCGCGCCAGGCTGTCGTTATCTTCTTCCATccg                                                   <  2:206809/251‑1 (MQ=255)
                                                                                                                                                                                                   cgcgATAGGCCGAGCCTGAGACTTGGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACTCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTTCCAGCCATAAAGGTAAAGCAGGTATTGCCGATATCGTGGAGCACCGCGCCGCCCCCGCTACTGCGCCGCGCCAGGCTGACTTTATCTTCTTCCATCCGCCgg                                               <  2:382819/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccTGAGACTTTGCGGTCGCCTTCGACGGTGTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTAGTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGTACGTTATCTTCTTCCATCCGCCGGGTATTACACTCttt                                 <  2:170614/251‑1 (MQ=255)
                                                                                                                                                                                                                    gagaCTTTGCGGTCGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTc                                <  1:62113/249‑1 (MQ=255)
                                                                                                                                                                                                                                 cGCCTTCGACGGTTTTCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGGTCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTgcgcg                 <  1:130638/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGAcc            >  1:227380/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gTTTTCAC‑TCCCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGAGCAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGCATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCCACTGCGCCGCGCCAGGCGGACGTAATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCCCGGGTTCTGCCCGCCACCACTtac      >  2:336665/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ttCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATtactac   >  1:68527/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tCAC‑TACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATtactac   <  1:336665/250‑1 (MQ=255)
                                                                                                                                                                                                                                                tCAC‑TACCAGATCGTTACGTCAGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGAAT‑TGAGCACAATAGACGTGGCGATAGTTTTATCGTACTCAGGCTTGCCAGCCATAAAGGTAAAGCAGGCATTGACGAGATCGTGGAAAACCGCGCCGCCAACGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCTCCGGGGATTACACTCTTTCCACGGGTTCTGCGCGCGAACAATTACTAcc  >  2:75781/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
CCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCAC‑CACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCAT‑TGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACC  >  NC_000913/4623499‑4623989

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