Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,517 T→A 12.0% intergenic (+42/‑206) ygaM → / → nrdH putative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,800,5170TA12.0% 208.1 / 21.1 75intergenic (+42/‑206)ygaM/nrdHputative membrane‑anchored DUF883 family ribosome‑binding protein/hydrogen donor for NrdEF electron transport system; glutaredoxin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (4/5);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-02

AGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCA  >  NC_000913/2800243‑2800766
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
aGTCCTGGGGAAGAGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCGCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCt                                                                                                                                                                                                                                                           <  1:154507/275‑1 (MQ=255)
    cTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGACTCTTTTGTTCGTGAAAGACCCTGGTTTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCAGGCGCAAATCGTAGTCCAAAAATGATAATAAATACGATTCTTTGACCCCGAAGCCTGTCt                                                                                                                                                                                                                                                       >  4:271468/1‑275 (MQ=255)
            gCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGAGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATAAGCGTCTTTTACCCCTAAGCCTGTCTGCGGGGtt                                                                                                                                                                                                                                               >  4:526811/1‑275 (MQ=255)
                                                      aGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTACTGCAAAAATGATAATAAATATGCGTCTTTGACCCCGAAGCCTGTCTTCGGGTTTTCTTTTTGCCTGGTGAATCACAACAATCCCCCTACACCGTCa                                                                                                                                                                                                     >  2:288583/1‑275 (MQ=255)
                                                        gCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTAACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTTTCTGGTNAATAANNAAAATccccccccaccctca                                                                                                                                                                                                     >  3:15322/1‑263 (MQ=255)
                                                           ttGCTGAAAGAAACCAGGGCGCGAAGGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                                                                                                                                                                                 <  2:526812/274‑1 (MQ=255)
                                                            tGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTAGTCGCGGCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                                                                                                                                     <  1:15322/269‑1 (MQ=255)
                                                             gCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATCCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAANNNCAAAAATCCCCCTACCCCCTcccccccat                                                                                                                                                                                              >  3:152750/1‑267 (MQ=255)
                                                             gCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCa                                                                                                                                                                                               >  2:12902/1‑274 (MQ=255)
                                                                  aagaaaCCCGCGCTCGAATGCATGGTCGTACCCGCGTCCAGCAAGCCGCGCGCGATGCCGGTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGTGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                                                                                                                                                                           <  1:72428/273‑1 (MQ=255)
                                                                      aaCCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGACAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGt                                                                                                                                                                                     >  4:271065/1‑275 (MQ=255)
                                                                           gggCGCGAATGCATGGTCGTACTCGCGTCCAGCCAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACACTGGTGTAGCGTGGGTACAGCAGCTGCGGTCGGTATTTTTATCGGCGCACTGTTAAGCATGCTCAAATCGTAGTGCAACCATGATAATAACAAACAGTCTTTGACCCCAAAGCCTTTCTTCGGTGTTTGTTTTTTCCTGGTGAATCCCAAAAATCCCCCTACCCCGTCCCCCTCCTACCCAggggtatt                                                                                                                                                                                 >  4:515630/1‑269 (MQ=255)
                                                                           gggCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGCAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTCCATGGTGAATCAAAAAAATCCCCCTACCCCGTCACCCCCATATCCGGGGTAAttt                                                                                                                                                                                >  4:421102/1‑275 (MQ=255)
                                                                                     gCATGTTCCTACTCCGGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGTGCTGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                      <  4:122332/275‑1 (MQ=255)
                                                                                     gCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTGGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                      <  2:515631/275‑1 (MQ=255)
                                                                                     gCATGGTCGTACTCGCGGCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATGTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                      <  2:152470/275‑1 (MQ=255)
                                                                                          gTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGc                                                                                                                                                                  <  2:271066/274‑1 (MQ=255)
                                                                                              tACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTACGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCAAAAAACTCCCCCTACCCCGTCACGCTCATATCTAGGGTAATTCCGACCa                                                                                                                                                                          >  4:57925/1‑262 (MQ=255)
                                                                                                  cgcgTCGAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                      <  2:57926/262‑1 (MQ=255)
                                                                                                     gTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt                                                                                                                                                      >  4:425809/1‑275 (MQ=255)
                                                                                                                    cgcgATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTCTTTGCTATATATTGTGTGGTTGAATCtttttt                                                                                                                                       >  4:152916/1‑275 (MQ=255)
                                                                                                                          gCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAAc                                                                                                                                   <  2:44079/273‑1 (MQ=255)
                                                                                                                                   tGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttttt                                                                                                                                       >  4:338295/1‑260 (MQ=255)
                                                                                                                                       cgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAAc                                                                                                                                   <  2:338296/260‑1 (MQ=255)
                                                                                                                                       cgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtctc                                                                                                                    <  1:152750/275‑1 (MQ=255)
                                                                                                                                       cgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgg                                                                                                                                                                                                                                                   <  5:274283/148‑1 (MQ=255)
                                                                                                                                         gATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAcc                                                                                                                                                                           <  1:104069/218‑1 (MQ=255)
                                                                                                                                               tttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCggg                                                                                                                                                                                                                                                  <  5:331988/141‑1 (MQ=255)
                                                                                                                                               tttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCggg                                                                                                                                                                                                                                                  >  6:331988/1‑141 (MQ=255)
                                                                                                                                                    tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACagag                                                                                                       >  1:33988/1‑275 (MQ=255)
                                                                                                                                                    tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCggg                                                                                                                                                                                                                                                  <  7:607197/136‑1 (MQ=255)
                                                                                                                                                    tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCggg                                                                                                                                                                                                                                                  >  8:607197/1‑136 (MQ=255)
                                                                                                                                                              cccTGGTGTAGCGTGGGTACAGCAGCTTCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCg                                                                                             <  2:152917/275‑1 (MQ=255)
                                                                                                                                                                cggggttaggctGGGTAACGCAGCTGCGGTATGTATTTTTATCGACGCACTGTTAAGCATGGGCAAATCCTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCTACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGAc                                                                                           <  4:39706/264‑1 (MQ=255)
                                                                                                                                                                         cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGAc                                                                                                   >  2:61687/1‑258 (MQ=255)
                                                                                                                                                                          gTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcaca                                                                                                                                                                                                                         >  3:145820/1‑139 (MQ=255)
                                                                                                                                                                           tGGGTTCGGCAGCTGCGGTAGGTATTTTTGTCGGCGCACTGTTAAACCGGCGCAAGTTGTAGTGCAAAAATGATATTAAATACGCGTCTTTGAGCCCGAAGCCTGGCTTCTGGGTTTCTTATTGCCTGATGAATCAGAAAAATCCCCCTGCCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGATTGAATGTTTTTTCAACTACATCTGGTATCTCTGTATCAACAGAGAAACAACCCTACGCGTatcat                                                                                  <  4:288582/273‑1 (MQ=255)
                                                                                                                                                                             ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAAcc                                                                                               <  4:61686/258‑1 (MQ=255)
                                                                                                                                                                              gTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                                                                                                                                                                                     <  1:145820/139‑1 (MQ=255)
                                                                                                                                                                                         gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATCTCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCa                                                                   >  2:65385/1‑274 (MQ=255)
                                                                                                                                                                                            gtaggtaTTTTTATCGGCTCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCa                                                                                                                                     >  2:459650/1‑205 (MQ=255)
                                                                                                                                                                                            gtaggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                                                                                                                                  >  2:491446/1‑192 (MQ=255)
                                                                                                                                                                                             taggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                                                                                                                                                                               <  5:898010/98‑1 (MQ=255)
                                                                                                                                                                                             taggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                                                                                                                                                                               >  6:898010/1‑98 (MQ=255)
                                                                                                                                                                                                gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                                                                                                                              <  4:491445/192‑1 (MQ=255)
                                                                                                                                                                                                    ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                                                                                                                                                                    >  7:774615/1‑86 (MQ=255)
                                                                                                                                                                                                    ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                                                                                                                                                                    <  8:774615/86‑1 (MQ=255)
                                                                                                                                                                                                           ggCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTCTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTCTATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGaa                                                >  4:85113/1‑275 (MQ=255)
                                                                                                                                                                                                             cgcACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTTTCAACAGAGAAACAACCCGACGCGTATCATCCCCCCGTATCTTTATTTTAAACGGAAATaacaacc                                              >  4:441922/1‑269 (MQ=255)
                                                                                                                                                                                                                ccttttAAGCATTCCCAAATCGTTGTTCNANAAGGATACTAAATCCGCGTCTTTGCCCCCGAATCCTCTCTTCGGGNTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACNACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATg                                           <  3:33988/271‑1 (MQ=255)
                                                                                                                                                                                                                  tGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATgcg                                         <  2:425810/275‑1 (MQ=255)
                                                                                                                                                                                                                     tAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCa                                                                                                                                                                                         <  8:8757/128‑1 (MQ=255)
                                                                                                                                                                                                                     tAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCa                                                                                                                                                                                         >  7:8757/1‑128 (MQ=255)
                                                                                                                                                                                                                      aaGCATGCGCAAATGGTAATGCAAAAATGATAATAAATACGCGTTTTTGACCCCGAAGCCTGCCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCGGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCAACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCAtt                                     <  2:208424/275‑1 (MQ=255)
                                                                                                                                                                                                                            gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTTTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACt                                                                                                                                                                        >  4:311700/1‑138 (MQ=255)
                                                                                                                                                                                                                            gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                                                                                                                                                                         >  2:547975/1‑105 (MQ=255)
                                                                                                                                                                                                                            gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTCTTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATATTCATTTTAAACGGAAATACGAATCATgcgc                                        >  2:60905/1‑266 (MQ=255)
                                                                                                                                                                                                                            gcgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGCTTCTTTTTGCCTGGTGAATCACAAAAATCTCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATATAGTATCTCTGTATCTACAGAGAGACAACCCGACGCGCATCATCGCCACGTCTTTTCATTTTAAACGGAAATACGTAACAATCGTattattatt                               >  4:208423/1‑275 (MQ=255)
                                                                                                                                                                                                                                caaTCGTAGTTCAAAAATGATAATAAATACGCGTCTTTTACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTa                                    <  4:60904/265‑1 (MQ=255)
                                                                                                                                                                                                                                aaaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTTTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAttt                                                                                                                                                                    <  2:311701/138‑1 (MQ=255)
                                                                                                                                                                                                                                aaaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                                                                                                                                     <  4:547974/105‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtcaa                                                                                                                                                                                                                                             <  8:258884/61‑4 (MQ=21)
                                                                                                                                                                                                                                    cGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtcaa                                                                                                                                                                                                                                             >  7:258884/1‑58 (MQ=25)
                                                                                                                                                                                                                                     gTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGt                      <  2:271469/275‑1 (MQ=255)
                                                                                                                                                                                                                                        gTGCAAAAATGATAATTAATAAGCGTCTTTGACCCCGAAGCCTTTCTTCGGTGTTTCTTTTTGCCTGGTTAATAACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAGTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTaa                    <  1:89227/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      ataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTg                                                                                                                                                 >  8:544612/1‑135 (MQ=255)
                                                                                                                                                                                                                                                      ataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTg                                                                                                                                                 <  7:544612/135‑1 (MQ=255)
                                                                                                                                                                                                                                                          ataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACAAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTgcca  <  2:441923/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTg                                                                                                                  >  3:107223/1‑153 (MQ=255)
                                                                                                                                                                                                                                                                       aCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATc                                                                                                              <  1:107223/153‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTg                                                                                                                                                 >  7:788385/1‑115 (MQ=255)
                                                                                                                                                                                                                                                                          ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTg                                                                                                                                                 <  8:788385/115‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTa                                                                                                                           <  6:653181/136‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTa                                                                                                                           >  5:653181/1‑136 (MQ=255)
                                                                                                                                                                                                                                                                             aaGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAAt                                               >  2:395875/1‑210 (MQ=255)
                                                                                                                                                                                                                                                                                 cTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATg                                           <  4:395874/210‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATcc                                                                                                                                                                                          >  6:346275/1‑66 (MQ=255)
                                                                                                                                                                                                                                                                                  tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATcc                                                                                                                                                                                          <  5:346275/66‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTa                                                                                                                                                                                    >  6:969693/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                                  tGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGCTa                                                                                                                                                                                    <  5:969693/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         
AGTCCTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCA  >  NC_000913/2800243‑2800766

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