Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,579,656 A→C 17.4% V62G (GTG→GGG)  symE ← toxic peptide regulated by antisense sRNA symR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,6560AC17.4% 155.6 / 2.0 75V62G (GTG→GGG) symEtoxic peptide regulated by antisense sRNA symR
Reads supporting (aligned to +/- strand):  ref base A (50/7);  new base C (6/10);  total (56/19)
Fisher's exact test for biased strand distribution p-value = 1.29e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGGTGTG  >  NC_000913/4579420‑4579888
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
cgttgaaaTTTTGACTGGGGGGGGAAGTTTCTTTTTGGTCAGGGAAAGTCTCCCGGTGGTTAAGCGGGTATCAACGCCGTTACGCTACTCTTTGTTCACCGGCAATCACTCCAATAAACGGCTGCCCCTGCTTTTGTTTACGCGGCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACt                                                                                                                                                                                                                                           <  2:890758/234‑1 (MQ=255)
   ggaaaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGCGAGGCACATACAGCCTCCCATCACctt                                                                                                                                                                                                                                         >  1:1006938/3‑238 (MQ=255)
   gTAAATTTTTTGTGGGGGGGAATTTTTTTTTTGGGTAGGGGGAAGGCCCCCGGTGGTTAACCGGGAGTAAACGCTTTTACGCGTCTTTCTTTTTACCGGCAATCACTCCAATAAACGCCTGCCCCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTg                                                                                                                                                                                                                                        <  2:926349/239‑1 (MQ=255)
         atTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATc                                                                                                                                                                                                                                  >  1:1007879/1‑238 (MQ=255)
             aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTCACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGCCAGTTTGCACACCGGCCGTAGCGACTGCATCCGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGGGGACCATACAGCCTTCCCTCCCTTTCACATccccc                                                                                                                                                                                                                              >  2:832531/1‑234 (MQ=255)
             aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTAcc                                                                                                                                                                                                                              >  1:54533/1‑238 (MQ=255)
             aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGGTGGGCGGTGAGGACAAAACAGCCTTCCATCACTTTGACATCTAcc                                                                                                                                                                                                                              >  2:378317/1‑238 (MQ=255)
                  gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGGGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGt                                                                                                                                                                                                                         >  1:486857/1‑238 (MQ=255)
                   ggggAAAGTTTCTATTGAGTCAGTGGAATGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGGGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTg                                                                                                                                                                                                                        >  1:935128/1‑238 (MQ=255)
                   ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGTCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCCCTTTGACACCTACCGCCGTg                                                                                                                                                                                                                        >  1:1075774/1‑238 (MQ=255)
                   ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTTACATCTACCGCCGTg                                                                                                                                                                                                                        >  1:709749/1‑238 (MQ=255)
                   ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGCGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTg                                                                                                                                                                                                                        >  1:122408/1‑238 (MQ=255)
                   ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCCCTCCTCGGCGGCGGGGGGTTGGGCGGTGAGGACAATACAGCCTTCCACCACTTTGACTTCTACCGCCGTg                                                                                                                                                                                                                        >  1:511187/1‑238 (MQ=255)
                   ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCAACCCTGGCGTAGCGACTCCATCATTTCGCCCTCCTCGGCGGCGGGTGGTTGGGCGGGGAGGGCAATACAGCCTTCCCTACCTTTGACACTCAccgcccgg                                                                                                                                                                                                                        >  2:265671/1‑235 (MQ=255)
                           tttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCGGTGGCa                                                                                                                                                                                                               >  1:693854/1‑239 (MQ=255)
                                   tcgtgtggggATGGTTCGCGGTGTTTACCCGTGAGGAAACGGTGTTACGGGTCTTTTTGTTTAGCGGCAATCATTTCAATGAACGGCTACCCCTGCTTTTGTTGACGCGCCGACGGTGGGCAGCCCTGGCGGAGCGACTGCATCAGTTCGCCCTCCCCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTCCCACCTCTTTGACATCTACCGCCGTGCCAGTGGCAAAACcggc                                                                                                                                                                                                       <  2:42875/231‑1 (MQ=255)
                                     gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGGGGCAAAACcggcg                                                                                                                                                                                                      >  1:694159/1‑238 (MQ=255)
                                      tCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACCGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTGGGGCGGTGCGGACAATACAGCCTTCCATCCCTTTTGACTCTACCGCCGTGCCAGGGGCAAAACccgcg                                                                                                                                                                                                      >  2:857130/1‑237 (MQ=255)
                                         gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGAtt                                                                                                                                                                                                  >  1:470786/1‑238 (MQ=255)
                                         gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGGGAGGACAAACCAGCCTTCCATCACTTTGACATCTCCCGCCGGGCCAGGGGCAAAACCGGGGGCtc                                                                                                                                                                                                  >  2:643024/1‑237 (MQ=255)
                                           ggAAGGCCTCCGGGTGTTTACCCGGGATAACACGCTTTTCGCCGTTTTTTTTTTTCCCGGCAATCCCCCCAATAACCGCCTGCACCGGTTTTTTTTTACCCGCCGACGTTTGGCCCCCCTGGCGTAGCGACTGCATCAGTTCCTTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAAAATCCACCCTCCCATCCTTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTcc                                                                                                                                                                                                <  2:873719/239‑1 (MQ=255)
                                                  cTCCCGTGGTTTAACGGGGTGTAAACTCGTTTCCGCGTTTTTTTTTGTACCGGCAGCCCCCCCAATAACCGCCGGCACGTGTTTTTGTTTACGCGCGGCACGGTGGCACCCCTGGCGTGGCGACTGCATCAGTCCGCCCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGAAAATACAGCCTTCCACCCCTTTGCCATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACt                                                                                                                                                                                         <  2:529107/239‑1 (MQ=255)
                                                   tCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGCCAATACAGCCTTCCCTCCCTTTTACATTTCCAACCCTTCCACTGCCCAAACCGGCACCTTCCCCCCAcac                                                                                                                                                                                        >  2:108347/1‑237 (MQ=255)
                                                   ccccGGTGGTTAACCGGGAGAAACCGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACCCCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTg                                                                                                                                                                                        <  1:378317/238‑1 (MQ=255)
                                                     ccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGCGGTTGGGCGGTGAGGACACTACAGCCTCCCATCACTTTGCCATCTACCGCCGGGCCAGTGGCAAAACCGGGGGCTTCCAccccccgc                                                                                                                                                                                       >  2:736146/1‑233 (MQ=255)
                                                     ccGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTACACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGa                                                                                                                                                                                       >  1:119183/1‑238 (MQ=255)
                                                               aCCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGt                                                                                                                                                                            >  1:958565/1‑239 (MQ=255)
                                                                             cttgtccccgATTTTCTGTTTACCGGCAATCCCCCCAATAACCGCCGGCCCCTGCTTTTGTTTACGCGCCGACAGTTTGCACCCCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATAc                                                                                                                                                              <  1:721855/231‑1 (MQ=255)
                                                                              tGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGCAATCTACCGCCGTGCCAGTGGCAAAACCGCGGGCCTCCCACCGCAGACCCTTCCGGGGGATGGCggggaaaag                                                                                                                                                             >  2:665616/1‑233 (MQ=255)
                                                                              tGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCATCGGCGGCGGGGGGTTGGGCGGTGAGGCCAAATCAGCCTTCCATCACTTTGACATCAACCGCCGTGCCAGTGGCAAAACCGCGGGCTTCCCGCCCCAGACATTTCCGGGGGGAGGGGGGAAAACg                                                                                                                                                             >  2:84054/1‑239 (MQ=255)
                                                                               gTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGGGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACg                                                                                                                                                             >  1:462987/1‑238 (MQ=255)
                                                                                   cgcgACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGGGAGGCCAATACAGCCTTCCATCCCTTTTACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGaaaagggggga                                                                                                                                                        >  2:1082988/1‑230 (MQ=255)
                                                                                        cTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGCCCTTTCAGGGTGATGGCGGGAATACGGCTGTAATcc                                                                                                                                                   >  1:339286/1‑239 (MQ=255)
                                                                                        cTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAgg                                                                                                                                                                              >  1:64213/1‑212 (MQ=255)
                                                                                        cTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAgg                                                                                                                                                                              >  1:256053/1‑212 (MQ=255)
                                                                                        cTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGACTACGGCTGTAATcc                                                                                                                                                   >  1:193940/1‑239 (MQ=255)
                                                                                        cTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATcc                                                                                                                                                   >  1:959403/1‑239 (MQ=255)
                                                                                          ttCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGCCCTTTCCGGGGGATGGCGGGAATACGGGGGTAAACCgg                                                                                                                                                 >  2:1075518/1‑239 (MQ=255)
                                                                                                  aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGTCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCCGCCCCTGCCCTTTCAGGGTGATGGCGGGAAAACCGGCGTAAACCCGGTGGcaa                                                                                                                                         >  2:635313/1‑236 (MQ=255)
                                                                                                  aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAgg                                                                                                                                                                              >  1:809830/1‑202 (MQ=255)
                                                                                                  aCCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCATGGTGATGGCGGGGATACGGCTGTAATCCGGGTAGCGa                                                                                                                                         >  1:880391/1‑238 (MQ=255)
                                                                                                       cAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCa                                                                                                                                   >  1:464867/1‑239 (MQ=255)
                                                                                                         aTCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCCGGGGGATGGCGGGAAAACGGCTGTACACCCGGTAGAGAcacacaaa                                                                                                                                 >  2:171985/1‑232 (MQ=255)
                                                                                                              tCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACCCGCATAAcct                                                                                                                             >  1:365709/1‑237 (MQ=255)
                                                                                                                        gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGCCGGGGGCATg                                                                                                                  >  1:1000026/1‑239 (MQ=255)
                                                                                                                        gCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGGAGCGACTCGCATAACCGCCGGTGACATg                                                                                                                  >  1:937637/1‑239 (MQ=255)
                                                                                                                          cTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCCAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGCATACGGCTGCAATCCGGGTAGCGCCTCGCCTAACCGACgggtccttgac                                                                                                                >  1:1620/1‑230 (MQ=255)
                                                                                                                                          ttccgcgCCGACAGTTGCCCCCCCGGGCGTGGCTACTGCATCATTTCGCTCTCCTCGGCGGGGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTGGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGga                                                                                                <  1:857130/236‑1 (MQ=255)
                                                                                                                                          ttccgccgcGACGGTTTGCACCCCTGGCGGAGCGCCTCACACGGTTCGCCCTCCCGGGCGGCGGGTGGTTGGTGGGTGAGGAAAATGCAGCCTTCCATCCCTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATTCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGga                                                                                                <  1:84054/231‑1 (MQ=255)
                                                                                                                                                        tttGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGCCGGGGACATGACGGTTATTTGCCCGGGAGACTTCTGCTTCgg                                                                                 >  1:778265/1‑238 (MQ=255)
                                                                                                                                                              cacCGGGCGTAGCGACGGACTAAGTTCGCTCTCCGCGGTGGGGGGTGGTTGGGCGGTGAGGACAATGCAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGtt                                                                           <  1:189184/239‑1 (MQ=255)
                                                                                                                                                                  tGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGaa                                                                                                                                                                 <  1:123540/151‑1 (MQ=255)
                                                                                                                                                                  tGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGaa                                                                                                                                                                 >  2:123540/1‑151 (MQ=255)
                                                                                                                                                                              tGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGGGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACCTCTGCTTCGAACGGTTGTGCAATacctggc                                                             >  1:861926/1‑231 (MQ=255)
                                                                                                                                                                               gCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGCGATGGCGGGAAAACGGCTGGAATCCggg                                                                                                                                                >  1:24066/1‑155 (MQ=255)
                                                                                                                                                                                cATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGGTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAAAATGCgtg                                                          >  1:156695/1‑238 (MQ=255)
                                                                                                                                                                                  tcggTTCTCTCCTCTCGGGGGGGGGTGGTGGGGGGGGGGGAAACACACGGCTTCCCACCCTTTTGACTTCCACCGGCGGGCCAGGGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTgtca                                                       <  2:1000026/236‑1 (MQ=255)
                                                                                                                                                                                       tCGCTCTCCGCGGCGGGGGGTGGTTGGGCGGTGAGGACAAGCCGGCCTTCCATCCCTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCAt                                                  <  1:573500/239‑1 (MQ=255)
                                                                                                                                                                                       cgattcccctgGGCGGCGGGTGGTTGGGGGGGAGGACCAGTGTCGCCTCCACACCTTTGCACCCCTCCCCCAGTGCCAGTGGCAAACCCGGGGGCTTCCAGCCACTGGCCTTGCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCAt                                                  <  2:1620/232‑1 (MQ=255)
                                                                                                                                                                                          ctctCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGGTGTAATCCGGGTAGCGACTCGGCTAACAGACGGGGACCAGACGGTTAATTTCCGGGGGGGCTCTCTGGTCCAAACGTTGGGCAATAACAAGCgtgtagcgtaaaa                                                >  2:209094/1‑229 (MQ=255)
                                                                                                                                                                                               ctcGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGGACGGTTGTGCAATAGAATGGGGGTCAGTCATAACTGCTa                                          >  1:464952/1‑239 (MQ=255)
                                                                                                                                                                                                tgggcgggggTGGGTTGGGCGGGGAGGACAATACCGCCTTCCATCTCTTTGCCACTGCCCGCCGGGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTAt                                         <  1:67347/237‑1 (MQ=255)
                                                                                                                                                                                                                 ggCGGTGAGGGCAATACAGCCTTCCATCACTTTGACCTCTACCGCCGTGCCAGGGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCCGGGGGATGGCGGGAATACGGCTTGTATCCGGGGACGCGCCCCCCTAAACGCACGGGGACAGGCGGTTATTTGCCGGGGGGGCTTTGTCTTCGTACGGGTGTGAAATAGAATGCGGGTCCGTCCACACTTCCATTTTCCAGGGAAAgggat                        >  2:154815/1‑239 (MQ=255)
                                                                                                                                                                                                                    gggtggagaaaaTCCACCCTCCCAACCCTTGGACACTGACCGGCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtgatt                       <  1:271191/228‑1 (MQ=255)
                                                                                                                                                                                                                      tGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATa                                                                   >  1:1052788/1‑191 (MQ=255)
                                                                                                                                                                                                                        aGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGCCGGTGACATGACGGGTATTTGCCGGGGGGACTTCTTCTTCGAACGGTTGTGCAATAGAATGCGGGTCAGCATTAACTGCTATTTTCCCGGAATAgggattgggttt                 >  2:795559/1‑234 (MQ=255)
                                                                                                                                                                                                                        aGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGGAGCGACTCGCATAACCGGCGGGGGCATGACGGTTATTTGCCGGGGAGAATTCTGCTTCGAACGGGTGTGCCATAGAATTGGGGTCAGCCATAACTGCTATTTCCCCGGAAAAgggatttgggtt                 >  2:756400/1‑233 (MQ=255)
                                                                                                                                                                                                                        aGGAAAATACCGCCTCCCATCCTTTTACCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAgtgattgtgat                  <  1:795559/238‑1 (MQ=255)
                                                                                                                                                                                                                                tACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATa                                                                   >  1:807645/1‑182 (MQ=255)
                                                                                                                                                                                                                                tACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTTTCCAGGAATAGTGATTGTGATTAGCGAtc          >  1:299321/1‑238 (MQ=255)
                                                                                                                                                                                                                                 acacCCTACCATCCTTTTGCCTTTCACCGCCGTGCCAGTGGCAAACCCGGGGGCTTCCAGCCACTGCACTTGCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGc         <  1:385485/239‑1 (MQ=255)
                                                                                                                                                                                                                                   cgCCTTCCATCCTTTTGACATCTGCCGGCGTGCCGGTGGCAAAACCGGGGGCTTCCAGCCAATGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCgg       <  1:108347/238‑1 (MQ=255)
                                                                                                                                                                                                                                     ccTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACgg                                                                                                                         >  1:682475/1‑124 (MQ=255)
                                                                                                                                                                                                                                        tCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGCACGGTTGTGCAATAGAATGCGGGTCAGTCATAACTGCTATTTTCCAGGAATAGGGATTGTGATTAGCGATGCggggggg  >  1:384197/1‑235 (MQ=255)
                                                                                                                                                                                                                                        tCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCCGTCATAACTGCTATTCTCCAGGAAAATGGATTGGGATTACCGATGCGGgggtg  >  1:506020/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
CTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAATAGAATGCGTGTCAGTCATAACTGCTATTCTCCAGGAATAGTGATTGTGATTAGCGATGCGGGTGTG  >  NC_000913/4579420‑4579888

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