Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,894,839 C→T intergenic (+68/+80) ygcP → / ← ygcQ putative anti‑terminator regulatory protein/putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,894,8390CT98.4% 192.6 / ‑4.5 64intergenic (+68/+80)ygcP/ygcQputative anti‑terminator regulatory protein/putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (36/25);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

CTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTA  >  NC_000913/2894595‑2895070
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
cTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGGAAATCCGCCAACGGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGCAATTCAGCGACTGa                                                                                                                                                                                                                                     >  2:176367/1‑249 (MQ=255)
  gCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCACCCGCTGATTTCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGCCGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGc                                                                                                                                                                                                                                  <  1:555256/250‑1 (MQ=255)
     aGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACg                                                                                                                                                                                                                              >  1:273931/1‑251 (MQ=255)
          gtaggCCCAAAGTGCTGGGCTGGGGGACAGAGAAAAGCCGCCAACCGTTGATGGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCTTGGCGCTTTCCACCACGAATACCGGGGTCGGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑TTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTg                                                                                                                                                                                                                          <  1:330831/246‑1 (MQ=255)
          gTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCtt                                                                                                                                                                                                                           <  1:53912/249‑1 (MQ=255)
                ccAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAAACAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCAACGCTTGACGCGTc                                                                                                                                                                                                                   >  2:677557/1‑251 (MQ=255)
                 cAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGCCGCGTCt                                                                                                                                                                                                                  >  1:441604/1‑251 (MQ=255)
                     gTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCACGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGAGTCTGGACATTAGCGAAAAAATTACTTTGACTGGATAATC‑GTAACCAATTGCATTTGGATTGATTTGTAGGCCGCACGCCACATCCGACATTCAGATTCTGATACGACGCTTGCAGCGTCTAATc                                                                                                                                                                                                              >  2:174211/1‑251 (MQ=255)
                         tgggctgggTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCt                                                                                                                                                                                                                            >  1:323845/1‑233 (MQ=255)
                         tgggcggggTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGTCGCTTCCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCt                                                                                                                                                                                                                            <  2:323845/233‑1 (MQ=255)
                          gggctgggTTACGGAGAAAACCCGCCATCCGCTGATTGCCGGTGGGCTGGGTTGCGATGAAGAATAGGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATTCCGGGGTCTGGACGTTAGCGAAAAAATTACTTGGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                                                                                                                                          <  2:642807/250‑1 (MQ=255)
                              tgggGGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                                                                                                                                      >  2:239807/1‑250 (MQ=255)
                                  tGACAGAGAAAATCCGCAACCCGCTGATTGCCGATGGGCTGGTGTGCGATGAAGCAGATGCGCGTAATACGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGACACAGCCGACATTCAGCGTCTGATGCGCCTCTTGCCGCGTCTTATCAGGCCTACAAGTc                                                                                                                                                                                                 >  2:434666/1‑251 (MQ=255)
                                             aaccGCCAACCGGGGAGTGCTGGTGGGCTGGTGTGCGATGAAGAAGATGCGTGTAAGGGGATTAACGCGGGTGTAGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATGGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGa                                                                                                                                                                                      <  2:273931/249‑1 (MQ=255)
                                                  ccAACCTCGGATTGCCGGTGGGTTGGGGGGCGTTGAAGAATATGCGTGTAAGGCGATTACCGCGGGTGTCGTGGCGCTTTCCACCACGAATTCCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTGGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTGTCAGGCCTACAAGTCCCCTGGCCAGAACCGt                                                                                                                                                                                 <  2:641133/251‑1 (MQ=255)
                                                           gatggCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGTACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGat                                                                                                                                                                        <  1:676790/247‑1 (MQ=255)
                                                                  ggggggCTGGTGTGCGATTAAGAAAATGCGTGTAAGTCTATTAACGCGGGTTTCGTGTCTCTGTCCACCACTAAGACCGGGGTCTGGACTTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTGGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTAGCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCgg                                                                                                                                                                 <  2:413308/248‑1 (MQ=255)
                                                                    tGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGCCCCCAGGCCAGAACCGTCGGCCGGATCATCCGGc                                                                                                                                                                >  2:332837/1‑250 (MQ=255)
                                                                              ggCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTCCTTTGACGGGATAATC‑TTAACCAATTGAATTTGGTTTGATGTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGAGGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTAtt                                                                                                                                                     <  1:354517/250‑1 (MQ=255)
                                                                                         agaTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCg                                                                                                                                          <  1:174211/251‑1 (MQ=255)
                                                                                           aTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                                                                                        >  1:173976/1‑251 (MQ=255)
                                                                                                    aTGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATACTC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCCGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGCTCATCCGCCAACACGTTATTGGCAATCTGCTTGAATGTTTGt                                                                                                                               >  2:74235/1‑251 (MQ=255)
                                                                                                        gATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                                                                                           >  1:573931/1‑251 (MQ=255)
                                                                                                         aTTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                                                                                           >  1:167409/1‑250 (MQ=255)
                                                                                                          ttAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGACAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                                                                                           >  1:96527/1‑249 (MQ=255)
                                                                                                            aaggcgGGTGTCGTGGCGCTTCCCTCCACGAATACCGGGGTCGGGACGTGAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTGGGTTTGATTTGTAGGCCGTACGCCAAATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTg                                                                                                                       <  1:190/248‑1 (MQ=255)
                                                                                                                       cGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGcac                                                                                                             <  1:332837/250‑1 (MQ=255)
                                                                                                                                     ccgcgaatctcGGGGTTTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATGTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtcatc                                                                                              <  1:49968/241‑1 (MQ=255)
                                                                                                                                      cacGAATACCGGGGTCTGGACGTTAGAGTAAAAATTACTTTGCCGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTCCGTGAATGTTTGTCACCAGCTCCTCCCGCACCACTTTCCAGtcctc                                                                                              >  2:108194/1‑247 (MQ=255)
                                                                                                                                         gcaTACCGGAGTCTGGACGATAGCGAAAAAATGACTTTGACGGGATAATC‑GAAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGCCATTCAGAGTCTGATGCAACGCATGACGCCTCTTATCAGGCCTACACATCCCCAGGCCAGAACCGTAGGCATGATCATCCGGAAACACGTTATGGGCAATCTGCGTGAAAGTTTGTCACAAGTGCATCCCGCACCACTTTCAAGTCATCAAcc                                                                                          >  2:497583/3‑251 (MQ=255)
                                                                                                                                         gAATACGGGGGTTTGGACGTTAGCGAAAAAATTACTTGGACGGGATAATC‑GTAACCAATTGAATTGGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAAGCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTGCCAGTCATCAAcc                                                                                          <  1:426587/251‑1 (MQ=255)
                                                                                                                                                ggggttttgacgttatttaaaaaaTTATTTTGACGG‑ATAACCTGAATCCAGTTGGATTTGGTTTGTTTTTTAGGTCCCACGCCACATCCGACATTCAGCGTCTTATGCGACGCTTGACGCGTCTTATCAGGCCTATAAGTCCCCAGTCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCTTGAATGTTTGTCACCAGTGCCTCCTGCACCCCTTTCCAGTCATCAACCACGctaa                                                                                   <  2:411153/233‑4 (MQ=255)
                                                                                                                                                 gggTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAAc                                                                                  >  1:346152/1‑251 (MQ=255)
                                                                                                                                                 gggTATGGATGTGAGCGAAAAAATTACTTAGACGGGATAAAT‑GTAACCAATAGAATTTGGTTTGATTTGTAGGCCGCACGCCACATACGACATTCAGCGGCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCAAGAACCGTAGGCAGGAGCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTGGTCAACAGTGGCGCTAGCACAACTTTCCAGTCATCAACCACGCCAAc                                                                                  <  2:71490/251‑1 (MQ=255)
                                                                                                                                                   gTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCAAGCACCACTTTCCAGTCATAAACCACGCCAACAt                                                                                >  2:29273/1‑251 (MQ=255)
                                                                                                                                                    tCTGGACGTTAGCGCAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCACCAGCCCCACTTTCCATCCATCAACCACGCCACCATc                                                                               >  2:316792/1‑251 (MQ=255)
                                                                                                                                                      tGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg                                                                             >  1:568706/1‑251 (MQ=255)
                                                                                                                                                             tAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtgtg                                                                       >  1:364969/1‑250 (MQ=255)
                                                                                                                                                              aGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATGTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGaa                                                                     <  1:380660/251‑1 (MQ=255)
                                                                                                                                                                          ttttgaaggggATAA‑C‑GTAACCAATTGAATTTGGTTTGATTTGTAGCCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAgcg                                                         <  2:111883/242‑1 (MQ=255)
                                                                                                                                                                               gACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGACCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTgg                                                     <  1:176367/250‑1 (MQ=255)
                                                                                                                                                                                  gggATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCCGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCACCGCTGGCGTc                                                 >  2:417542/1‑251 (MQ=255)
                                                                                                                                                                                      tAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATg                                              >  1:494204/1‑250 (MQ=255)
                                                                                                                                                                                        aTC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGAtt                                           <  1:453048/251‑1 (MQ=255)
                                                                                                                                                                                              aaCCAATTGATTTGGGTTTGATCTGTAGGCAGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATc                                       <  2:285984/250‑1 (MQ=255)
                                                                                                                                                                                              aaCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCCACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATcc                                      >  2:156503/1‑250 (MQ=255)
                                                                                                                                                                                              aaCCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACAACGCCAAAATCTGCCTGTGAAAATCCCGCAGCGCTGGCGTCAAGATTAATCg                                      >  2:470017/1‑251 (MQ=255)
                                                                                                                                                                                                  aaTTGAATTTGGTTTGATTTGTAGAACGAACGCAAACTCCCAAATTCAGCGTATGATGCGACGCATGTCTCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTGGGCAGGATCATCCCGCAACACGTTATTGGCAATCTGCGCGAATGTTTAACCCCAATGCATCCAGCAACCCTTTCCAGCCCTCAACCACGCCACCACCTGCCTGTGAACATACCGCACCCATGTCGTCATGCTTCATCGccac                                  >  2:339041/1‑251 (MQ=255)
                                                                                                                                                                                                   aTTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGccacc                                 >  1:233263/1‑251 (MQ=255)
                                                                                                                                                                                                    ttGAATTTGGTTTGATATGAAGGCCGCAAGCCACATCCGACATTAAGCGTCGGATGCCAAGCTTGACGCCGCTTATCAGGCCTACCAGTCCCCAGGCCAGACCCGTAGGCAGGATCATCCGGCAACACGTTATTGGAAATATGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAACACCGCAGCTATGGCGTCATGATTCATCGccacc                                 >  2:237807/1‑250 (MQ=255)
                                                                                                                                                                                                     tGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACAGCAGCGCTGGCGTCATGATTAATCGCcaccac                               >  1:560133/1‑251 (MQ=255)
                                                                                                                                                                                                            cccgccgctcctccgatctaCGCCACATCCGACATTCAGCGTCGGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTg                        <  2:356847/232‑1 (MQ=255)
                                                                                                                                                                                                               ttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATc                                                                                                                                              >  2:108641/1‑130 (MQ=255)
                                                                                                                                                                                                               ttGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATc                                                                                                                                              <  1:108641/130‑1 (MQ=255)
                                                                                                                                                                                                                    ttGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGc                 <  1:434666/250‑1 (MQ=255)
                                                                                                                                                                                                                     tGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGt               >  1:164067/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCAACCACAAACTTGCTATTGCGt               >  1:469352/1‑251 (MQ=255)
                                                                                                                                                                                                                      gTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTa              <  1:239807/251‑1 (MQ=255)
                                                                                                                                                                                                                        aGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTACAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcg            >  2:413147/1‑251 (MQ=255)
                                                                                                                                                                                                                            cGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcgcccg        >  1:603504/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcgcccgc       >  2:336656/1‑251 (MQ=255)
                                                                                                                                                                                                                               aCGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTg                        >  1:356847/1‑232 (MQ=255)
                                                                                                                                                                                                                                 gCCACATCCGACATTCAGCGTCTGATGCGAAGCTTGACGCGTCTTATCAGGCCTACAAGTCCCAAGGCAAGAACAGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGAGCCTCCAGCACCACTTTCCAGTAATCAACCACGCCAACATCTGCCTGTGAAAATACCGAAGCGATGGCGTCATGGTTAATCGCCACCACATACTTGCTATTGCGTACGCCCGCCAtt   >  1:458702/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cacaTCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTa  >  1:209536/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
CTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATC‑GTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTA  >  NC_000913/2894595‑2895070

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: ‑