Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,334,010 C→T intergenic (‑4213/+326) yfaQ ← / ← yfaT tandem DUF2300 domain‑containing protein YfaQ/DUF1175 domain‑containing protein YfaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,334,0100CT81.0% 72.1 / 13.4 43intergenic (‑4213/+326)yfaQ/yfaTtandem DUF2300 domain‑containing protein YfaQ/DUF1175 domain‑containing protein YfaT
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (21/13);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

CTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGC  >  NC_000913/2333785‑2334252
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
ctgcGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGGTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAAAGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGt                                                                                                                                                                                                                           >  1:190607/1‑251 (MQ=255)
    ggCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGTGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGcc                                                                                                                                                                                                                       <  2:522073/251‑1 (MQ=255)
                              tGGATAACGGAACTGTTTCACCAGGGGCTATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAAAGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCAACAGTCAGGTTAGCGTATTGTTCAGCCCGTAGCCCAGATATTGCGGTTGCACCACAATGc                                                                                                                                                                                             >  2:567673/1‑251 (MQ=255)
                                 aTAACGGAACGGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATATCGGGAAGGTTTAATAATGGCAGTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTTAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTGAGCGTATCGTTTAGCCCGGCGCCCAGATAGTGCGGTTGCACCACAATGCGAt                                                                                                                                                                                          <  2:211808/251‑1 (MQ=255)
                                  tAACGGAACTGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATg                                                                                                                                                                                         >  1:163382/1‑251 (MQ=255)
                                     cGGAACTGTTTCCCCAGGGGATATTTTTTCAGCGGGGAAAACTGGCTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAg                                                                                                                                                                                      <  1:230283/251‑1 (MQ=255)
                                                                  cAGCGGGGAAAACTGGTTGTTCTGTACATAGCGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTTCACCACAATGCGATGCATGTTTTTCTGCTGTCGCAAAAATGCCATc                                                                                                                                                          <  2:402297/250‑1 (MQ=255)
                                                                     cGGGGAAAACTGGTTGTTCTGTACATAACGGGACGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTAAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGCTATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCCTCGGGt                                                                                                                                                      >  2:426532/1‑251 (MQ=255)
                                                                         gAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAgg                                                                                                                                                  >  1:187942/1‑251 (MQ=255)
                                                                                 gttgttCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGAt                                                                                                                                          <  1:396827/251‑1 (MQ=255)
                                                                                   tgttCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATa                                                                                                                                         <  2:36270/250‑1 (MQ=255)
                                                                                    gttCTGTACATAAAGGGAAGGTTTAATAATGGCATTGCCGAGCTGTACTTCGGGTAACGCCTGAGTCACATTCTGCCGTGCCTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAACTGCCATCGGGTCAGGAATACGATAc                                                                                                                                        >  1:609263/1‑250 (MQ=255)
                                                                                          tACATGACGGGAAGGTTTAATAATGGCATTGCCGAGCTGAAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGGCTGAGAAGAGAAAGGACGCAGCATCACGCGGCGAGATTTGCCGAACCAGTTATCCGACAGTCAGGTTAGAGTATTGTTCAGCCCGGCTCCCAGATATTGCGATTGCACCACAATGCGACGCCGGTTTTTCTGCTGGCGATAGATTGCCATCGGGTCAGGAAGACGATACAGGCGa                                                                                                                                  >  1:106833/1‑250 (MQ=255)
                                                                                                    gAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGCGTCACATTCTGCCGTGACTGAGAAGAGCAAGTACGCTGCCTCACGCGGCGAGATTTGCCGTACCAGTTATCCCACCGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTTGCGCCAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCAcgcc                                                                                                                       >  1:330044/1‑251 (MQ=255)
                                                                                                    gAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCAcgcc                                                                                                                       >  1:397796/1‑251 (MQ=255)
                                                                                                    gAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCAcgcc                                                                                                                       >  2:449850/1‑251 (MQ=255)
                                                                                                          ttAATAAGGGCATTGCCTAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTGGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGAGATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGTCGCAAAAATGCCATCGGGTCATGAATGCGATACAGGCGAACATCCACGCCGCCGt                                                                                                                  <  2:609263/250‑1 (MQ=255)
                                                                                                                    cATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATc                                                                                                       <  1:243517/251‑1 (MQ=255)
                                                                                                                        gCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGa                                                                                                    <  1:190455/250‑1 (MQ=255)
                                                                                                                                   aTTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACATGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGAtc                                                                                         >  2:504937/1‑249 (MQ=255)
                                                                                                                                      cggggaccGCCTGAGTCACATTCGGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATGGCGGTTGCACCACAAGGGGATGCAGGTTTTCCTGCGGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGGGAACTTCCACGCCGCCGTACGCTTCCAGCTGATAGCGCCGTTAATCa                                                                                     <  2:77793/244‑1 (MQ=255)
                                                                                                                                      cGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCAGCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCa                                                                                     <  1:394135/251‑1 (MQ=255)
                                                                                                                                       gggTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCAc                                                                                    <  1:487402/251‑1 (MQ=255)
                                                                                                                                         atAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGCATACGATACAGGCGAACGTCCACGCCGCCGTAATCTTCCATCTGATACCGCCGATAATCACGc                                                                                  >  2:353646/2‑251 (MQ=255)
                                                                                                                                                     tCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTcc                                                                      >  1:159274/1‑251 (MQ=255)
                                                                                                                                                     tCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTcc                                                                      >  1:277753/1‑251 (MQ=255)
                                                                                                                                                         aTTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGCACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTc                                                                  >  2:503908/1‑251 (MQ=255)
                                                                                                                                                              gCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGAAAGGATAGCGTATTGTTCAGACCGTCGCCCAGATATTGCGGTTGCACCACAATACGATGCAGGTTTCTCTGATGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATAACACCCCGGCGCTTCCAGTCGCACt                                                             >  1:328835/1‑251 (MQ=255)
                                                                                                                                                                 gTGACTGAGAAGAGAAAGTCCGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTc                                                          <  2:544841/251‑1 (MQ=255)
                                                                                                                                                                          aagagaaAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTc                                                  >  1:277592/1‑250 (MQ=255)
                                                                                                                                                                           agagaaAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCAc                                                >  1:142803/1‑251 (MQ=255)
                                                                                                                                                                             agaaAGTACGCGGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCActg                                              <  1:254292/251‑1 (MQ=255)
                                                                                                                                                                               aaaGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCActgct                                            <  2:20204/251‑1 (MQ=255)
                                                                                                                                                                                    aCGCTGCATCACGCGGCGAGATTTGCCGTCCCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGACATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATAAGGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGagag                                                                                                 <  2:164552/193‑4 (MQ=255)
                                                                                                                                                                                    aCGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGGATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAg                                                                                                 >  1:164552/1‑193 (MQ=255)
                                                                                                                                                                                        tGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCCGCCCGTCCCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGTGCGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTacaaac                                   >  2:230268/1‑246 (MQ=255)
                                                                                                                                                                                            tCACGCGGCGAGATTTGCCGTACCAGTTATCCCGCAGTCAGGTGAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCCCGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAActgct                                <  2:337620/250‑1 (MQ=255)
                                                                                                                                                                                            tCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGCTACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAActgctg                               >  2:292390/1‑251 (MQ=255)
                                                                                                                                                                                                          ttGCCGTACCAGTTATCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAACTGCTGTCAGCAAGCaaaa                  >  1:385342/1‑250 (MQ=255)
                                                                                                                                                                                                                tACCAATTATACAACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCACATATTGCGGATGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACGGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCCCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTCAAACTGCTGTCAGCAAGCAAAAAGAATGt           >  2:260837/1‑251 (MQ=255)
                                                                                                                                                                                                                 aCCAGGTATCCCACAGTCAGGATAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAAATGCTGTCAGCAAGCAAAAAGAATGtt          >  1:205847/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTCCCACAGTCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGGTCCCCCgg   >  1:127516/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccACAGTCAGGTTAGCGCATTGTTCAGCCCGTAGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAACTGCTGTCAGCAAGCAACAAGAATGTTCCCCCGGc  >  1:527783/1‑249 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
CTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCACCAGGGGATATTTTTTCAGCGGGGAAAACTGGTTGTTCTGTACATAACGGGAAGGTTTAATAATGGCATTGCCGAGCTGTAATTCGGGTAACGCCTGAGTCACATTCTGCCGTGACTGAGAAGAGAAAGTACGCTGCATCACGCGGCGAGATTTGCCGTACCAGTTATCCCACAGCCAGGTTAGCGTATTGTTCAGCCCGTCGCCCAGATATTGCGGTTGCACCACAATGCGATGCAGGTTTTTCTGCTGGCGCAAAAATGCCATCGGGTCAGGAATACGATACAGGCGAACGTCCACGCCGCCGTACTCTTCCATCTGATAGCGCCGATAATCACGCCCCGGCGCTTCCAGTCGCACTTTCGCCTCTTCACTGCTGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGC  >  NC_000913/2333785‑2334252

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: