Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,519 T→G 11.8% intergenic (+44/‑204) 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,5190TG11.8% 215.7 / 24.6 81intergenic (+44/‑204)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 (35/35);  new base G (4/6);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-02

CTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCA  >  NC_000913/2800247‑2800766
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
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                                                                                                                                                                                                                                               >  6:898010/1‑98 (MQ=255)
                                                                                                                                                                                         taggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGtt                                                                                                                                                                                                                                               <  5:898010/98‑1 (MQ=255)
                                                                                                                                                                                            gtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                                                                                                                              <  4:491445/192‑1 (MQ=255)
                                                                                                                                                                                                ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                                                                                                                                                                    <  8:774615/86‑1 (MQ=255)
                                                                                                                                                                                                ttttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                                                                                                                                                                    >  7:774615/1‑86 (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                                                                                                                                                                                                                                             >  7:258884/1‑58 (MQ=25)
                                                                                                                                                                                                                                cGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCGGGGtcaa                                                                                                                                                                                                                                             <  8:258884/61‑4 (MQ=21)
                                                                                                                                                                                                                                 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)
                                                                                                                                                                                                                                                                               gTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCa                                                                                                                                     <  6:251654/118‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCa                                                                                                                                     >  5:251654/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
CTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCA  >  NC_000913/2800247‑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: