Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 160,462 C→A R48L (CGT→CTT)  dksA ← RNA polymerase‑binding transcription factor DksA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913160,4620CA83.6% 130.0 / 9.3 55R48L (CGT→CTT) dksARNA polymerase‑binding transcription factor DksA
Reads supporting (aligned to +/- strand):  ref base C (6/3);  new base A (20/26);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CGCTTCCAGACGGCGAATACCAATTTCAACACCGCAGGATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGT  >  NC_000913/160221‑160705
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
cgcTTCCAGACGGCGAATACCAATTTCAACACCGCAGGATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCtt                                                                                                                                                                                                                                            <  1:184084/251‑1 (MQ=255)
           ggCGAATACCAATTTCAACACCGCAGGATTCGCAGTAGCCGAAATCGTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATacga                                                                                                                                                                                                                                 <  2:429961/251‑1 (MQ=255)
                         tCAACACCGCAGGATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGcc                                                                                                                                                                                                                     >  2:446546/1‑249 (MQ=255)
                            acacCGCAGGATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGAGCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCGTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTg                                                                                                                                                                                                                <  2:305439/251‑1 (MQ=255)
                               ccGCAGGATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGc                                                                                                                                                                                                             >  2:59909/1‑251 (MQ=255)
                                    ggATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCttcat                                                                                                                                                                                                        >  2:99524/1‑251 (MQ=255)
                                    ggATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCttcat                                                                                                                                                                                                        >  2:135703/1‑251 (MQ=255)
                                     gATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGCCTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCttcatt                                                                                                                                                                                                       >  2:621718/1‑251 (MQ=255)
                                             gTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTGTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACTACGGAAGGGCGCCAGCTGGGCTTCATTCATATAct                                                                                                                                                                                               <  2:269621/251‑1 (MQ=255)
                                                    aaaTCTTCGTCTTCCTCTTTTGTCAGCGTCTTCGCGATCTTTGTGATCAGCTTACGCTCGCGAGCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGGCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCTTATACTCTTCGcc                                                                                                                                                                                        <  1:84015/251‑1 (MQ=255)
                                                          tCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAAATCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCtt                                                                                                                                                                                  >  2:503563/1‑251 (MQ=255)
                                                                       ttttCATCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTGCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGTCTTCTCCTGATATgg                                                                                                                                                                      <  2:388409/250‑1 (MQ=255)
                                                                         ttCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGtt                                                                                                                                                                    <  1:38375/250‑1 (MQ=255)
                                                                           cAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGTCTTCTCCTGATATGGTTcc                                                                                                                                                                  <  2:627981/250‑1 (MQ=255)
                                                                               gTCTTCTCGATCTTTTTGATCAGCTTAAGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCAcccc                                                                                                                                                             >  2:410967/1‑251 (MQ=255)
                                                                                      cGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGa                                                                                                                                                        <  1:538700/249‑1 (MQ=255)
                                                                                      cGATCTGTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCg                                                                                                                                                         <  2:396715/248‑1 (MQ=255)
                                                                                           tttttGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCg                                                                                                                                                   >  1:593063/1‑249 (MQ=255)
                                                                                              ttGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAAt                                                                                                                                              >  2:453970/1‑251 (MQ=255)
                                                                                               tGATCAGCTTACGCTCGCTATCGAGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCgagact                                                                                                                                              >  2:450666/1‑248 (MQ=255)
                                                                                                gATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATAc                                                                                                                                            >  2:180494/1‑251 (MQ=255)
                                                                                                     gCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCg                                                                                                                                                         >  1:25240/1‑233 (MQ=255)
                                                                                                     gCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCg                                                                                                                                                         <  2:25240/233‑1 (MQ=255)
                                                                                                         aCGCTCGCGCTCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGCAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACg                                                                                                                                   >  2:620271/1‑251 (MQ=255)
                                                                                                             tCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAgg                                                                                                                                       >  1:608659/1‑243 (MQ=255)
                                                                                                             tCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAgg                                                                                                                                       <  2:608659/243‑1 (MQ=255)
                                                                                                               gcgATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGttt                                                                                                                             <  1:126794/251‑1 (MQ=255)
                                                                                                                cgATCGAGGGTACGCAGTTCGAGGCTGAACTATTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGCAGTTGGCTGCTTCATCCTGCCTATGTGTAACGGTGCGATCGACTTCATTCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGtttt                                                                                                                            >  2:418567/1‑251 (MQ=255)
                                                                                                                     gcgGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGt                                                                                                                       <  2:87532/251‑1 (MQ=255)
                                                                                                                             gCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCct                                                                                                               <  1:564282/251‑1 (MQ=255)
                                                                                                                               aGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCcttc                                                                                                             >  1:71542/1‑251 (MQ=255)
                                                                                                                                    gAGGCTGATCTCTTCGTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATTCTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCa                                                                                                        <  2:329033/251‑1 (MQ=255)
                                                                                                                                          gAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTg                                                                                                   <  1:453970/250‑1 (MQ=255)
                                                                                                                                                                      aCCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGTAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGGTTTACGGTGTTGCCCTTCTTGCATGTTGCTTCGCCTTAACACTCACTAGCGATc                                                                         <  2:564629/248‑1 (MQ=255)
                                                                                                                                                                       ccGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTTATTACGCCATGCTTCCAGAATACGACGGAAGTGCTCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTACCGATCCCCa                                                                     <  2:285895/251‑1 (MQ=255)
                                                                                                                                                                             ccGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCg                                                               >  1:627600/1‑251 (MQ=255)
                                                                                                                                                                                 gAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCggggg                                                           <  1:563744/251‑1 (MQ=255)
                                                                                                                                                                                 gAAGGTGGCCGCTTCTTCCTGCATATGTGTTACGGTGCGATCGACTGCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGTTTCAGCCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCTCCCCAGCGATGGCGAGAATTCTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCTTGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCggggg                                                           <  2:331668/251‑1 (MQ=255)
                                                                                                                                                                                    gTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGaaa                                                        <  2:337303/251‑1 (MQ=255)
                                                                                                                                                                                    gTTGGATGCTTCATCCTGCATATGTGTACCGGTGCGATCGACTTCCTCCCTGAGTTGATTACGCCATGATTCCCGAATACGACGGAAGTGAGCCAGCTGGGCTTCATTCCTATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCAAGCGTTGGCGAGCATACTCAGGGACGTTGTTTTACGGTTTTGCCCATCTTGCATGTTGCTTCTCCTTACCACGCTCTATCGATCCCCATGTTCGGGGGaa                                                         >  2:290181/1‑250 (MQ=255)
                                                                                                                                                                                       ggCTGCTTCATTCTGCATATGTGTAACGGTGCGATCGACTTCATCCCGGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGTGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGaaaaa                                                      <  2:223425/250‑1 (MQ=255)
                                                                                                                                                                                        gctgctTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATg                                                    >  1:121601/1‑251 (MQ=255)
                                                                                                                                                                                             ttCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGcc                                               >  1:37828/1‑251 (MQ=255)
                                                                                                                                                                                                aTCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCt                                            <  1:59909/251‑1 (MQ=255)
                                                                                                                                                                                                               aaCGGTGCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCATGCTATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTCCGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGc                             >  2:26042/1‑251 (MQ=255)
                                                                                                                                                                                                                     gCGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTcc                       >  1:86770/1‑251 (MQ=255)
                                                                                                                                                                                                                      cGATTGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGCCGGAAGTGCGCCAGCTGGGCGTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTCCGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTGCGGGGGAAAAATGAGGCCGCTATAAATAGCATATGTTTTTCCg                      <  2:87434/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGATCGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCGTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCTATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTGCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCg                      <  2:251988/251‑1 (MQ=255)
                                                                                                                                                                                                                        aTCGACTTCATCCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCgg                     >  1:15446/1‑250 (MQ=255)
                                                                                                                                                                                                                          cGACTTCATCCCTGAGTTGATTAAGTCATGCTTCCAGAATACGACGGAAGTGCGCAAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCTCCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCGTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGAt                   <  1:498/250‑1 (MQ=255)
                                                                                                                                                                                                                          cGACTTCATCCCTGAGTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGAt                   <  2:156054/250‑1 (MQ=255)
                                                                                                                                                                                                                                cATCCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAAtt            >  1:463804/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cccTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCATTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATc         >  2:487458/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gTTGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACg   <  2:490868/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ttGATTAAGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGt  >  1:523785/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
CGCTTCCAGACGGCGAATACCAATTTCAACACCGCAGGATTCGCAGTAGCCGAAATCTTCGTCTTCCACTTTTTTCAGCGTCTTCTCGATCTTTTTGATCAGCTTACGCTCGCGATCGCGGTTACGCAGTTCGAGGCTGAACTCTTCTTCCTGGGCTGCACGGTCTACCGGGTCCGGGAAGTTGGCTGCTTCATCCTGCATATGTGTAACGGTGCGATCGACTTCATCCCTGAGTTGATTACGCCATGCTTCCAGAATACGACGGAAGTGCGCCAGCTGGGCTTCATTCATATACTCTTCGCCCGGCTTCTCCTGATATGGTTCCACCCCAGCGATGGCGAGAATACTCAGGGACGATGTTTTACGGTTTTGCCCTTCTTGCATGTTGCTTCTCCTTAACACGCACTATCGATCCCCATGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGT  >  NC_000913/160221‑160705

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