Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,146,389 (T)6→7 coding (132/879 nt) pflC → putative pyruvate formate‑lyase 2 activating enzyme PflC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,146,3891.T94.1% 172.3 / 2.4 51coding (132/879 nt)pflCputative pyruvate formate‑lyase 2 activating enzyme PflC
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (22/26);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑T‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGG  >  NC_000913/4146147‑4146631
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
aCGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCt                                                                                                                                                                                                                                             <  2:16541/251‑1 (MQ=255)
           ttGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCtt                                                                                                                                                                                                                                  >  1:261962/1‑251 (MQ=255)
              cGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGc                                                                                                                                                                                                                               <  2:62494/251‑1 (MQ=255)
              cGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGc                                                                                                                                                                                                                               <  2:139347/251‑1 (MQ=255)
                  aaGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCgt                                                                                                                                                                                                                           <  2:8142/251‑1 (MQ=255)
                  aaGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATACGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCgg                                                                                                                                                                                                                           >  1:380286/1‑250 (MQ=255)
                      gcAACAGCGCCCACAAGATTATGCCGGGCTGGTGCTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGCTCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCgtggtg                                                                                                                                                                                                                       >  1:339647/1‑251 (MQ=255)
                              gcCCACAATATTATGCCGGGCTGGTGGTGCGCGTTGCCTGATACAGCGCCTGCTTTGTCGAACGGTCGAAGGAGAGCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACTCGCCGCGATGATGTGGCGCGCATTTTCAACATTGATCGTTAGTCACTGAATGGCGGCGTGGGTATTCGTACGGGGGTGTTTTTT‑T‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATc                                                                                                                                                                                                               <  1:38163/251‑1 (MQ=255)
                               cccACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATcc                                                                                                                                                                                                              >  1:31175/1‑251 (MQ=255)
                                ccACAAGAGTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCg                                                                                                                                                                                                             <  2:120448/251‑1 (MQ=255)
                                        ttATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGAt                                                                                                                                                                                                     >  1:4999/1‑251 (MQ=255)
                                         tATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGAtc                                                                                                                                                                                                    >  1:98224/1‑251 (MQ=255)
                                                       ggtgCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCa                                                                                                                                                                                      <  1:286718/251‑1 (MQ=255)
                                                           cgcgTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACg                                                                                                                                                                                  <  2:303065/251‑1 (MQ=255)
                                                                 gCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGc                                                                                                                                                                            >  2:14022/1‑251 (MQ=255)
                                                                 gCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGc                                                                                                                                                                            >  2:502/1‑251 (MQ=255)
                                                                    ggATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCaga                                                                                                                                                                         <  1:54808/251‑1 (MQ=255)
                                                                     gATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAAATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCag                                                                                                                                                                          >  1:164456/1‑249 (MQ=255)
                                                                        aCAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAgg                                                                                                                                                                     >  2:311024/1‑251 (MQ=255)
                                                                          aGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGTTGGT‑CTTTTTGT‑AAAGGCTGTCCGCAGCTTTGCTGGTGGTGGGCTAATCCGGAGTCGATATCCGGCAAAAGCCGGTCCATACGCagagagtcg                                                                                                                                                                   <  2:156556/251‑4 (MQ=255)
                                                                                   tttGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTg                                                                                                                                                          >  1:209903/1‑251 (MQ=255)
                                                                                          aaCTGTCGAAGGAGATCCAGGATAACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACtgtg                                                                                                                                                   <  1:215535/251‑1 (MQ=255)
                                                                                           actgtcgaaagcgataatgcaTGACATCATCCGCCGGACCGCGCATCAGATGTCACGTTGTGGAAACGCGCCGAGATCATGTGGCGCGAATTTTCAACATTCAGCGTTCTTCACTGAATGACGGCGAGGACTTTCGTACGGTGGTATTTTTT‑T‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGc                                                                                                                                                  >  2:376232/21‑251 (MQ=255)
                                                                                            cTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATAATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCg                                                                                                                                                 >  2:185067/1‑251 (MQ=255)
                                                                                                   aGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGtt                                                                                                                                          <  2:115402/251‑1 (MQ=255)
                                                                                                          ccAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTg                                                                                                                                    <  2:219243/250‑1 (MQ=255)
                                                                                                                             cgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGAGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTcc                                                                                                                <  2:112301/251‑1 (MQ=255)
                                                                                                                              ggACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCg                                                                                                               >  1:109954/1‑251 (MQ=255)
                                                                                                                              ggACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCg                                                                                                               >  1:150859/1‑251 (MQ=255)
                                                                                                                                aCAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCggg                                                                                                             <  1:201075/251‑1 (MQ=255)
                                                                                                                                              gTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCgtg                                                                                                                                                                                                                          >  2:196598/1‑128 (MQ=255)
                                                                                                                                              gTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCgtg                                                                                                                                                                                                                          <  1:196598/128‑1 (MQ=255)
                                                                                                                                                      gtgGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTcgcg                                                                                       <  2:31175/251‑1 (MQ=255)
                                                                                                                                                            aCGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCa                                                                                 <  2:333188/251‑1 (MQ=255)
                                                                                                                                                            aCGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAAGGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGGCGCGATATCa                                                                                 <  2:232400/251‑1 (MQ=255)
                                                                                                                                                               cgccgcGATGATGTGGCGCTCATTGTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGcc                                                                              <  2:361733/251‑1 (MQ=255)
                                                                                                                                                                gccgcGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCAATGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCt                                                                             >  1:20241/1‑251 (MQ=255)
                                                                                                                                                                              gcgcgcATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAAc                                                               <  1:160512/251‑1 (MQ=255)
                                                                                                                                                                                    aTTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAACGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGa                                                           >  2:95252/1‑249 (MQ=255)
                                                                                                                                                                                    aTTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAg                                                         <  1:326830/251‑1 (MQ=255)
                                                                                                                                                                                    aTTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAg                                                         <  2:19056/251‑1 (MQ=255)
                                                                                                                                                                                         cAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCGTGACGCTCTGGAACGGGAAGtgatg                                                    <  1:292138/251‑1 (MQ=255)
                                                                                                                                                                                           aCATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGaa                                                  <  2:261962/251‑1 (MQ=255)
                                                                                                                                                                                                  gCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGAc                                           <  2:98224/251‑1 (MQ=255)
                                                                                                                                                                                                         tCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATtttct                                    >  1:251483/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTcc                         <  1:286601/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTcc                         <  1:101981/251‑1 (MQ=255)
                                                                                                                                                                                                                      gCGAGGGCATTCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCcgg                       >  1:123807/1‑251 (MQ=255)
                                                                                                                                                                                                                               ttCGTACGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGt              >  1:32884/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aCGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGc         >  2:254211/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cGGTGGT‑CTTTTT‑TTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCt        >  2:255285/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ggtctttttt‑t‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCt    >  1:298472/6‑251 (MQ=255)
                                                                                                                                                                                                                                              ttttt‑t‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTgg  <  2:54728/248‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
ACGCCGACACGTTGCGGGAAGCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTT‑T‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCGCGATATCAGCCTTGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTCTTTCGCACGTCCGGCGGCGGCGTCACGCTTTCTGG  >  NC_000913/4146147‑4146631

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