Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,776,715 A→T S377C (AGT→TGT)  ydiF → putative acetate‑CoA transferase YdiF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,776,7150AT92.9% 99.9 / ‑5.1 44S377C (AGT→TGT) ydiFputative acetate‑CoA transferase YdiF
Reads supporting (aligned to +/- strand):  ref base A (1/1);  major base T (19/20);  minor base G (3/0);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAAATCGCTC  >  NC_000913/1776477‑1776958
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
cGAAATGCGTAAAGGCGCGGTGGGGACTGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGGGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGGACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGAGGTCTGAATGTTTGGTATTTGGGTTTTggtgt                                                                                                                                                                                                                                           >  2:322331/1‑245 (MQ=255)
     ggCGTAAAGGCGCGGGGGGTAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGGGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGat                                                                                                                                                                                                                                    <  2:227717/250‑2 (MQ=255)
           agggCGCGGGGGGGAATGGCGGCGGCGGTTTTGCTGACGTGATTTGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATTGCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAATGTGAATACCCGTGCCGTTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGAGGCGGTCTTGATGTTTGTTATTTGTGTTTTGCGGAAGTCGACCTGCAc                                                                                                                                                                                                                              <  2:389826/249‑1 (MQ=255)
               cgcgGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGt                                                                                                                                                                                                                           >  1:281187/1‑250 (MQ=255)
                 cgGTGGGGCATGTCGGGGTCGGTATTTATGACGGCATTTGCCTGGTCTGCCGAGAAAAAGGTTGTGCGTATGACTGTATGCTGACGGGAGAAACATGTCCGATTGGCTGAATTTCTTCAAAGGGGATCGCCTATGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGTCGTCCCAGTTTTATTTTTATCACGTTGGCGGTCTGGATGCTTGTTATTGGTGTTTTGCTTAAGTCTACCAGCCCGGTAAc                                                                                                                                                                                                                        <  2:314828/251‑1 (MQ=255)
                                     ggttgtgcgtactgccttGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCa                                                                                                                                                                                                     <  2:326658/234‑1 (MQ=255)
                                     ggaaTTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGCACGTGAATACCCGTGCCATTCTGGATCTGACGTCCCAGTTTGGTTTTTATCACGGTGGCGGTCTGGATGTTTGTTTTTTGTGTTTTGCTGAAGTCGACCAGCACGGTTACGTCGGCGTGCTTAAATTcc                                                                                                                                                                                                     >  2:136240/4‑249 (MQ=255)
                                       tATTGCTGACGGCATTGGCCTGGTCGCCCGAGTAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGGTATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGGTGTTTGATATTTGGTTTTTGCTGAAGTCGACCAGCACGGACACGTCAGCGTTGATAAATTCact                                                                                                                                                                                                   >  2:98144/1‑248 (MQ=255)
                                                    aTTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATggg                                                                                                                                                                                     >  1:129473/1‑251 (MQ=255)
                                                                gccagagaCGAAGGATGTGCTGATGACTTTATTCTGACGGGAGAAACAGGTCCGCTTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCATATTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGCAGTCGACCAGCACGGTAACGTAGGAGTGCATAAATTCAATTGGAAAATCAAGGGCACCGGTgggtt                                                                                                                                                                         >  2:251642/5‑248 (MQ=255)
                                                                 cccGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTg                                                                                                                                                                                                                                               <  1:76497/180‑1 (MQ=255)
                                                                 cccGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTg                                                                                                                                                                                                                                               >  2:76497/1‑180 (MQ=255)
                                                                    gagaAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCa                                                                                                                                                                                                             <  2:432776/211‑1 (MQ=255)
                                                                    gagaAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCa                                                                                                                                                                                                             >  1:432776/1‑211 (MQ=255)
                                                                         gaagGTTGCGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGCTTGGCGGAATTACTTCATACTGGATAGCCTTTGCCGCGAACGAGACTACCCGTGCCATTCTGGGTATGTCGTCCCAGATTGATTTTTATCACGGTGGCTGTCTGGATGTTTGTTATTTGTGTTTTGCTGAATGCGACCAGCCCGGTAACGTCGGAGGGCATAATTTCAATGGTAAAATCCTGGGCACCGGGGGTTTTATTGatat                                                                                                                                                                >  2:272269/1‑251 (MQ=255)
                                                                          aagGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATc                                                                                                                                                               <  1:175600/251‑1 (MQ=255)
                                                                             gTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAACTTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCatg                                                                                                                                                            >  1:404395/1‑249 (MQ=255)
                                                                                gtgCTGATGACTTTCTTCTGACGGTATAAACAGGTCCGACTGGCGGAATTACTTCACAGGGGAGCGCCTTTAGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGGCTGGATGGTTGTTATGTGTGGTTTGCTGAAGTCGGCCAGCACGGTAACGGCGGCGAGCATAAATTCAATGGTAAAATCAGGGGCCACGGGGGATTTGTTGTTATCAgggac                                                                                                                                                         >  2:52686/1‑247 (MQ=255)
                                                                                     gatgaCTTTATTCTGACGGTAGAAACAGGTCCGATTGTCGGAATTACTGCACAGGGGTGCGCCTTTGGCGCGAACTTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACtt                                                                                                                                                     <  2:390992/250‑1 (MQ=255)
                                                                                                   gACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCAt                                                                                                                                        <  1:264464/249‑1 (MQ=255)
                                                                                                    aCGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCATTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGCAGTCGACCAGCACGGTAACGTCGGAGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTAGAAGAAAATCTtttt                                                                                                                                     >  2:466331/1‑251 (MQ=255)
                                                                                                               aGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTTCTGAAGTCGACCAGCa                                                                                                                                                                                                                               >  2:112331/1‑150 (MQ=255)
                                                                                                               aGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCa                                                                                                                                                                                                                               <  1:112331/150‑1 (MQ=255)
                                                                                                                 gTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTa                                                                                                                        <  1:136240/251‑1 (MQ=255)
                                                                                                                  tcagATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGAGGCGGTCTGGATGTTTGTTATTTGGGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATTGTAAAATCATGGGCACCGGTGGATTTATTGATATCAATGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTa                                                                                                                        >  2:367680/4‑250 (MQ=255)
                                                                                                                          gCGGAATTACTTCACAGTGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGTGTTTGTTGATATCAGTGCCACTTCGACGTAAATCATTTTCTGCGGCACATTAACTGCggg                                                                                                                >  2:115012/1‑250 (MQ=255)
                                                                                                                                ttACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGCTGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCAACGGTGGATTTATTGATCTCAGTGCTACTTCGACGACAATCATTTTCTGCGGCCCCTTAACTGCGGGCAGttt                                                                                                          >  2:120237/1‑250 (MQ=255)
                                                                                                                                  aCTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTTAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATGCAATGGTAAAATCATGGGCACCGGTGGATTTATTCATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTaa                                                                                                        <  1:52686/250‑1 (MQ=255)
                                                                                                                                                  ccTTTGGCGCGAACGTGAATACCCTTGCCATTCGGGATATGGCGTGCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTGTTGCGGAAGACGACCAGCCCGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGTTGGATTTATTGATATCAGTTCCACTTCGAAGAAAATCATTTTCTGCGTCACATTAACTGCGGGCAGTTTAAAAACATAAATTACCGAc                                                                                       <  2:449657/251‑1 (MQ=255)
                                                                                                                                                             aaCGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTaa                                                                             >  2:165443/1‑250 (MQ=255)
                                                                                                                                                              aCGTGAATTCCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTTGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAat                                                                           <  1:290486/251‑1 (MQ=255)
                                                                                                                                                                       cctaggCCAGTCTGGAGATTACGTCCTAGTTTTATTTTTATCACTGTGGCGGTCTTGATGGTTGTTATTTGTTGTTTGCTGAAGTCGACCAGCACGGTAACGTCGTCGTGCATAAATTCAATGGTAAAAGCATGGGCACCGGTGGAGTTATTTATTTCAGTGCCACGTCGAAGAAAAGCATTTTCGTCGGCACATTAACCGCGGGCAGTTTAAAAGCAGAAATGACCGACGGCAATTTAAATATCGTCCag                                                                  <  2:242394/246‑1 (MQ=255)
                                                                                                                                                                       cccGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCag                                                                  <  1:165443/251‑1 (MQ=255)
                                                                                                                                                                       cccGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCag                                                                  <  1:120237/251‑1 (MQ=255)
                                                                                                                                                                                ttCTGGATATGACGTCCCAGTTTGATTTTTATCACGATGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAATGACgg                                                         >  1:410013/1‑251 (MQ=255)
                                                                                                                                                                                    ggATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTTTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGTCggg                                                        <  1:453855/248‑1 (MQ=255)
                                                                                                                                                                                                  aGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCgg                                       >  1:363063/1‑251 (MQ=255)
                                                                                                                                                                                                     ttGATTTTGATCACGGAGGCGGTCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGAAAAATCATGGGCACCGGTGGATTTATTGATATCAGTACCACTTCGAAGAAAATCATTTTCTGCGTCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAACTTAAGTATCGTCAAGGAAGGACGGGTGAAGAAATTTATTCGGGa                                     >  2:120176/1‑250 (MQ=255)
                                                                                                                                                                                                          ttttATCACGGTGGCGGTCTGGATTTTTGTTATTTGTGTGTTGCTGAAGTCGACCAGCACGGTAACGTCGTCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTTGATTTATTGATATCAGGGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCTGGCAGTTTAAAAACAGAAATTACCGATGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTat                                <  2:332050/250‑2 (MQ=255)
                                                                                                                                                                                                            ttATCCCGGGGGCGGGCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTACTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATGTAGTCGGGAACTACCg                              <  2:278362/250‑1 (MQ=255)
                                                                                                                                                                                                               tCACGGTGGCGGTCTGGGTGTTTGTTATTTGTGTTTTGCTGAAGTCGAACAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACCTTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCggca                           >  2:446169/1‑248 (MQ=255)
                                                                                                                                                                                                                     tGGCGGGCTGGATGTTTGTTATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTACAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACtt                    >  2:371541/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gtttgttATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGaaaaa        <  2:363063/250‑1 (MQ=255)
                                                                                                                                                                                                                                       ttATTTGTGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGCAATTACTTTCAGCGGAAAAATCGctc  >  1:345060/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
CGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAAAATCATGGGCACCGGTGGATTTATTGATATCAGTGCCACTTCGAAGAAAATCATTTTCTGCGGCACATTAACTGCGGGCAGTTTAAAAACAGAAATTACCGACGGCAAATTAAATATCGTCCAGGAAGGACGGGTGAAGAAATTTATTCGGGAACTACCGGAAATTACTTTCAGCGGAAAAATCGCTC  >  NC_000913/1776477‑1776958

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