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

CTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTA  >  NC_000913/2894595‑2895070
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
cTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGGAAATCCGCCAACGGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGCAATTCAGCGACTGa                                                                                                                                                                                                                                     >  2:176367/1‑249 (MQ=255)
  gCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCACCCGCTGATTTCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGCCGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGc                                                                                                                                                                                                                                  <  1:555256/250‑1 (MQ=255)
     aGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACg                                                                                                                                                                                                                              >  1:273931/1‑251 (MQ=255)
          gtaggCCCAAAGTGCTGGGCTGGGGGACAGAGAAAAGCCGCCAACCGTTGATGGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCTTGGCGCTTTCCACCACGAATACCGGGGTCGGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCTTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTg                                                                                                                                                                                                                          <  1:330831/246‑1 (MQ=255)
          gTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCtt                                                                                                                                                                                                                           <  1:53912/249‑1 (MQ=255)
                ccAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAAACAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCAACGCTTGACGCGTc                                                                                                                                                                                                                   >  2:677557/1‑251 (MQ=255)
                 cAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGCCGCGTCt                                                                                                                                                                                                                  >  1:441604/1‑251 (MQ=255)
                     gTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCACGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGAGTCTGGACATTAGCGAAAAAATTACTTTGACTGGATAATCGTAACCAATTGCATTTGGATTGATTTGTAGGCCGCACGCCACATCCGACATTCAGATTCTGATACGACGCTTGCAGCGTCTAATc                                                                                                                                                                                                              >  2:174211/1‑251 (MQ=255)
                         tgggctgggTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCt                                                                                                                                                                                                                            >  1:323845/1‑233 (MQ=255)
                         tgggcggggTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGTCGCTTCCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCt                                                                                                                                                                                                                            <  2:323845/233‑1 (MQ=255)
                          gggctgggTTACGGAGAAAACCCGCCATCCGCTGATTGCCGGTGGGCTGGGTTGCGATGAAGAATAGGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATTCCGGGGTCTGGACGTTAGCGAAAAAATTACTTGGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                                                                                                                                          <  2:642807/250‑1 (MQ=255)
                              tgggGGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                                                                                                                                      >  2:239807/1‑250 (MQ=255)
                                  tGACAGAGAAAATCCGCAACCCGCTGATTGCCGATGGGCTGGTGTGCGATGAAGCAGATGCGCGTAATACGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGACACAGCCGACATTCAGCGTCTGATGCGCCTCTTGCCGCGTCTTATCAGGCCTACAAGTc                                                                                                                                                                                                 >  2:434666/1‑251 (MQ=255)
                                             aaccGCCAACCGGGGAGTGCTGGTGGGCTGGTGTGCGATGAAGAAGATGCGTGTAAGGGGATTAACGCGGGTGTAGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATGGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGa                                                                                                                                                                                      <  2:273931/249‑1 (MQ=255)
                                                  ccAACCTCGGATTGCCGGTGGGTTGGGGGGCGTTGAAGAATATGCGTGTAAGGCGATTACCGCGGGTGTCGTGGCGCTTTCCACCACGAATTCCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTGGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTGTCAGGCCTACAAGTCCCCTGGCCAGAACCGt                                                                                                                                                                                 <  2:641133/251‑1 (MQ=255)
                                                           gatggCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGTACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGat                                                                                                                                                                        <  1:676790/247‑1 (MQ=255)
                                                                  ggggggCTGGTGTGCGATTAAGAAAATGCGTGTAAGTCTATTAACGCGGGTTTCGTGTCTCTGTCCACCACTAAGACCGGGGTCTGGACTTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTGGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTAGCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCgg                                                                                                                                                                 <  2:413308/248‑1 (MQ=255)
                                                                    tGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGCCCCCAGGCCAGAACCGTCGGCCGGATCATCCGGc                                                                                                                                                                >  2:332837/1‑250 (MQ=255)
                                                                              ggCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTCCTTTGACGGGATAATCTTAACCAATTGAATTTGGTTTGATGTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGAGGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTAtt                                                                                                                                                     <  1:354517/250‑1 (MQ=255)
                                                                                         agaTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCg                                                                                                                                          <  1:174211/251‑1 (MQ=255)
                                                                                           aTGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTg                                                                                                                                        >  1:173976/1‑251 (MQ=255)
                                                                                                    aTGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATACTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCCGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGCTCATCCGCCAACACGTTATTGGCAATCTGCTTGAATGTTTGt                                                                                                                               >  2:74235/1‑251 (MQ=255)
                                                                                                        gATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                                                                                           >  1:573931/1‑251 (MQ=255)
                                                                                                         aTTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                                                                                           >  1:167409/1‑250 (MQ=255)
                                                                                                          ttAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGACAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCAcc                                                                                                                           >  1:96527/1‑249 (MQ=255)
                                                                                                            aaggcgGGTGTCGTGGCGCTTCCCTCCACGAATACCGGGGTCGGGACGTGAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTGGGTTTGATTTGTAGGCCGTACGCCAAATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTg                                                                                                                       <  1:190/248‑1 (MQ=255)
                                                                                                                       cGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGcac                                                                                                             <  1:332837/250‑1 (MQ=255)
                                                                                                                                     ccgcgaatctcGGGGTTTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATGTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGtcatc                                                                                              <  1:49968/241‑1 (MQ=255)
                                                                                                                                      cacGAATACCGGGGTCTGGACGTTAGAGTAAAAATTACTTTGCCGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTCCGTGAATGTTTGTCACCAGCTCCTCCCGCACCACTTTCCAGtcctc                                                                                              >  2:108194/1‑247 (MQ=255)
                                                                                                                                         gcaTACCGGAGTCTGGACGATAGCGAAAAAATGACTTTGACGGGATAATCGAAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGCCATTCAGAGTCTGATGCAACGCATGACGCCTCTTATCAGGCCTACACATCCCCAGGCCAGAACCGTAGGCATGATCATCCGGAAACACGTTATGGGCAATCTGCGTGAAAGTTTGTCACAAGTGCATCCCGCACCACTTTCAAGTCATCAAcc                                                                                          >  2:497583/3‑251 (MQ=255)
                                                                                                                                         gAATACGGGGGTTTGGACGTTAGCGAAAAAATTACTTGGACGGGATAATCGTAACCAATTGAATTGGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAAGCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTGCCAGTCATCAAcc                                                                                          <  1:426587/251‑1 (MQ=255)
                                                                                                                                                ggggttttgacgttatttaaaaaaTTATTTTGACGGATAACCTGAATCCAGTTGGATTTGGTTTGTTTTTTAGGTCCCACGCCACATCCGACATTCAGCGTCTTATGCGACGCTTGACGCGTCTTATCAGGCCTATAAGTCCCCAGTCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCTTGAATGTTTGTCACCAGTGCCTCCTGCACCCCTTTCCAGTCATCAACCACGctaa                                                                                   <  2:411153/233‑4 (MQ=255)
                                                                                                                                                 gggTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAAc                                                                                  >  1:346152/1‑251 (MQ=255)
                                                                                                                                                 gggTATGGATGTGAGCGAAAAAATTACTTAGACGGGATAAATGTAACCAATAGAATTTGGTTTGATTTGTAGGCCGCACGCCACATACGACATTCAGCGGCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCAAGAACCGTAGGCAGGAGCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTGGTCAACAGTGGCGCTAGCACAACTTTCCAGTCATCAACCACGCCAAc                                                                                  <  2:71490/251‑1 (MQ=255)
                                                                                                                                                   gTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCAAGCACCACTTTCCAGTCATAAACCACGCCAACAt                                                                                >  2:29273/1‑251 (MQ=255)
                                                                                                                                                    tCTGGACGTTAGCGCAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCACCAGCCCCACTTTCCATCCATCAACCACGCCACCATc                                                                               >  2:316792/1‑251 (MQ=255)
                                                                                                                                                      tGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTg                                                                             >  1:568706/1‑251 (MQ=255)
                                                                                                                                                             tAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCtgtg                                                                       >  1:364969/1‑250 (MQ=255)
                                                                                                                                                              aGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATGTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGaa                                                                     <  1:380660/251‑1 (MQ=255)
                                                                                                                                                                          ttttgaaggggATAACGTAACCAATTGAATTTGGTTTGATTTGTAGCCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAgcg                                                         <  2:111883/242‑1 (MQ=255)
                                                                                                                                                                               gACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGACCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTgg                                                     <  1:176367/250‑1 (MQ=255)
                                                                                                                                                                                  gggATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCCGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCACCGCTGGCGTc                                                 >  2:417542/1‑251 (MQ=255)
                                                                                                                                                                                      tAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATg                                              >  1:494204/1‑250 (MQ=255)
                                                                                                                                                                                        aTCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGTCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGAtt                                           <  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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
CTGCCAGGCTGTATGCCCAAAGTGCTGGGCTGGGTGACAGAGAAAATCCGCCAACCGCTGATTGCCGGTGGGCTGGTGTGCGATGAAGAAGATGCGCGTAATGCGATTAACGCGGGTGTCGTGGCGCTTTCCACCACGAATACCGGGGTCTGGACGTTAGCGAAAAAATTACTTTGACGGGATAATCGTAACCAATTGAATTTGGTTTGATTTGTAGGCCGCACGCCACATCCGACATTCAGCGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCCCCAGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTA  >  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: