Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,800,517 T→A 5.5% 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,5170TA 5.5% 904.6 / 29.7 254intergenic (+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 (116/124);  new base A (8/6);  total (124/130)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-01

Only 100 of 259 total aligned reads displayed.

CTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCACGCCACCAAACGGG  >  NC_000913/2800247‑2800780
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                       
                                 ggaaagcaagaaaGCTCGGGCATTTCGGAAAGAAACCAGGGCGGGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaa                                                                                                                                                                                                                                     <  4:422547/264‑1 (MQ=255)
                                                  aGGCATTGCTGAAAGAAACCAGGGGTCGAATGCATGGTAGTACTCGCGTCCAGCAAGCCGGGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGGAAGACCCTGGTGTGGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCa                                                                                                                                                                                                                   <  4:242103/275‑1 (MQ=255)
                                                  aGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATca                                                                                                                                                                                                                                         >  2:422551/1‑253 (MQ=255)
                                                      aTTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGt                                                                                                                                                                                                                                                              <  2:212722/228‑1 (MQ=255)
                                                              aagaaaCCAGGGCGCGAATGCAGGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATc                                                                                                                                                                                                         <  2:499856/273‑1 (MQ=255)
                                                                       gggCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCNGCGCACTGTTAAGCATGCGAAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGCAGCNTGTCTTCGGGGTTTCTTTTTTCCTGGTGAATCACAAAAATCCCCCNANCCCGTTCCGCTCCTATCACGGGTAAtattga                                                                                                                                                                                           >  3:79903/1‑273 (MQ=255)
                                                                        ggCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  2:128296/225‑1 (MQ=255)
                                                                                         gTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGAc                                                                                                                                                                                          >  2:369778/1‑261 (MQ=255)
                                                                                                                  cgATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCa                                                                                                                                                   <  2:307192/275‑1 (MQ=255)
                                                                                                                      gCCGGTGGCGGTGCGGATTCTTTTGTTCGTTAACGACCCTGTTCTGGCTTGTGTACAGCATCTGCGGTAGGTATGTTTATCGGCGCACTGTTAAGCATGCGCAAATAGTAGTGCAAAAAGTATAATAAATACGCGTCTTTGACACCGAAGCCTGTCTTCGGGGTTTCTTTTTTCCTGGGGAAGCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATGTTTTGTTTGGTTGAAACTTTTTTCaa                                                                                                                                                  <  4:599097/272‑1 (MQ=255)
                                                                                                                            ggCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGt                                                                                                                                                                                                                                                                       <  5:1313805/149‑1 (MQ=255)
                                                                                                                               tGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCtt                                                                                                                                                                                                                                                                    >  8:1057427/1‑149 (MQ=255)
                                                                                                                                 ggcgGGTTCTTTGGGTTCTGTAAGACCCCGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTAACGGCGAACTGTAAAGCAGACGCAAAACGTAGTGCAAAAATGATAATAAAGACGCGTCGTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTAGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTAtc                                                                                                                                    <  4:199616/274‑1 (MQ=255)
                                                                                                                                 cgcgGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCg                                                                                                                                                                                                                                                                  <  6:841278/149‑1 (MQ=255)
                                                                                                                                    ggATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGACAGGCTTCgggg                                                                                                                                                                                                                                                               <  8:827204/149‑1 (MQ=255)
                                                                                                                                       ttCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    <  1:116295/221‑1 (MQ=255)
                                                                                                                                          ttttGTTGGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGAAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGAGGTTTCtt                                                                                                                                                                                                                                                         >  5:1514080/1‑149 (MQ=255)
                                                                                                                                          ttttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtt                                                                                                                                                                                                                                                         <  6:26085/149‑1 (MQ=255)
                                                                                                                                            ttGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCtttt                                                                                                                                                                                                                                                       <  6:980155/149‑1 (MQ=255)
                                                                                                                                                tCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATcaca                                                                                                                                                                                                                                       >  4:83479/1‑161 (MQ=255)
                                                                                                                                                  gTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTg                                                                                                                                                                                                                                                 <  8:404851/149‑1 (MQ=255)
                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               >  8:1603342/1‑149 (MQ=255)
                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  8:681656/149‑1 (MQ=255)
                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  8:752160/149‑1 (MQ=255)
                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                                                                                                                                                                                                   <  2:83479/161‑1 (MQ=255)
                                                                                                                                                    gAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcc                                                                                                                                                                                                                         >  2:64509/1‑171 (MQ=255)
                                                                                                                                                     aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                                                                                                                              >  8:166649/1‑149 (MQ=255)
                                                                                                                                                     aaaGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTg                                                                                                                                                                                                                                              <  5:41242/149‑1 (MQ=255)
                                                                                                                                                             tGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTGTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaa                                                                                                                                                                                                                                      <  7:1125253/149‑1 (MQ=255)
                                                                                                                                                                tgtAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaaaa                                                                                                                                                                                                                                   <  7:279824/149‑1 (MQ=255)
                                                                                                                                                                    gCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCTCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAAt                                                                                                                                                            >  1:179302/1‑216 (MQ=255)
                                                                                                                                                                     cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                                                                                                                                                                                                              <  8:1556822/149‑1 (MQ=255)
                                                                                                                                                                     cGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATcccc                                                                                                                                                                                                                              <  8:180215/149‑1 (MQ=255)
                                                                                                                                                                         ggTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAc                                                                                                                                                                                                                          <  6:499560/149‑1 (MQ=255)
                                                                                                                                                                           tACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAccc                                                                                                                                                                                                                        >  5:1393411/1‑149 (MQ=255)
                                                                                                                                                                            aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                                                                                                                                                                                       <  5:1009021/149‑1 (MQ=255)
                                                                                                                                                                            aCAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTAcccc                                                                                                                                                                                                                       >  8:1329585/1‑149 (MQ=255)
                                                                                                                                                                             cagcagCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    <  2:158730/183‑1 (MQ=255)
                                                                                                                                                                              agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGACTCACAAAAATCCCCCTACCCCGTCACGCTCa                                                                                                                                                                                                             >  4:497590/1‑157 (MQ=255)
                                                                                                                                                                              agcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGt                                                                                                                                                                                                                     >  5:772727/1‑149 (MQ=255)
                                                                                                                                                                               gcagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTc                                                                                                                                                                                                                    >  5:1065393/1‑149 (MQ=255)
                                                                                                                                                                                cagcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGa                                                                                                                                                                                                                                             >  7:1531328/1‑123 (MQ=255)
                                                                                                                                                                                 agcTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCAc                                                                                                                                                                                                                  <  6:1612916/149‑1 (MQ=255)
                                                                                                                                                                                   cTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACaaa                                                                                                                                                                                                                                     >  3:149311/1‑128 (MQ=255)
                                                                                                                                                                                    tGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCt                                                                                                                                                                                                               <  7:615952/149‑1 (MQ=255)
                                                                                                                                                                                     gCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTAt                                                                                                                                                                                    >  1:39395/1‑175 (MQ=255)
                                                                                                                                                                                      cGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               <  8:599683/115‑1 (MQ=255)
                                                                                                                                                                                         taggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCatat                                                                                                                                                                                                          >  8:1472819/1‑149 (MQ=255)
                                                                                                                                                                                          aggtaTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGt                                                                                                                                                                                                                                               >  7:368330/1‑111 (MQ=255)
                                                                                                                                                                                             taTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAg                                                                                                                                                                                                      >  5:365457/1‑149 (MQ=255)
                                                                                                                                                                                              aTTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAgg                                                                                                                                                                                                     <  8:738576/149‑1 (MQ=255)
                                                                                                                                                                                                  ttATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGaaa                                                                    <  2:280703/274‑1 (MQ=255)
                                                                                                                                                                                                       ggCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCg                                                                                                                                                                                            <  6:1570753/149‑1 (MQ=255)
                                                                                                                                                                                                           cACTGTTAAGCATGAGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATc                                                            <  2:448012/273‑1 (MQ=255)
                                                                                                                                                                                                            aCTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCAc                                                                                                                                                                                       >  7:1357501/1‑149 (MQ=255)
                                                                                                                                                                                                                ttAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCAAAAAAATCCCACTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCAATAtt                                                                                                                                                                                   >  7:1580197/1‑149 (MQ=255)
                                                                                                                                                                                                                    gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               >  8:381444/1‑149 (MQ=255)
                                                                                                                                                                                                                    gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               >  7:1205580/1‑149 (MQ=255)
                                                                                                                                                                                                                    gCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               >  5:941667/1‑149 (MQ=255)
                                                                                                                                                                                                                         cgcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCtatata                                                                                                                                                                          <  6:1065393/149‑1 (MQ=255)
                                                                                                                                                                                                                          gcAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAAtt                                                                                                                                                                                               >  8:1000221/1‑127 (MQ=255)
                                                                                                                                                                                                                            aaaTCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGa                                                                      <  2:478538/246‑1 (MQ=255)
                                                                                                                                                                                                                               tCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgt                                                                                                                                                                    >  7:1314527/1‑149 (MQ=255)
                                                                                                                                                                                                                                cGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACTCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtgtgtg                                                                                                                                                                   <  6:1459986/149‑1 (MQ=255)
                                                                                                                                                                                                                                   aGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTCTGTGGtt                                                                                                                                                                <  6:649773/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCt                                                                                                                                                                               <  4:40605/130‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACTCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGtt                                                                                                                                                                >  5:299361/1‑144 (MQ=255)
                                                                                                                                                                                                                                         aaaaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCt                                                                                                                                                          >  6:1212845/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                                                                                                                                                        <  5:1642502/149‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                                                                                                                                                        >  7:700601/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aTGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttt                                                                                                                                                       <  7:775174/149‑1 (MQ=255)
                                                                                                                                                                                                                                              gATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttttt                                                                                                                                                     <  7:435386/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 aataaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                                                                                                                                                  <  6:1167512/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   taaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACt                                                                                                                                                >  8:104282/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATca                                                                                                                                                                                                                                         <  8:1143787/59‑1 (MQ=255)
                                                                                                                                                                                                                                                    aaataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATca                                                                                                                                                                                                                                         >  7:1143787/1‑59 (MQ=255)
                                                                                                                                                                                                                                                      ataCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACa                                                                                                                                             <  6:173324/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       taCGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACAt                                                                                                                                            <  5:21939/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           cgTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAg                                                                                                                                        <  8:272138/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            gTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGc                                                                                                                                                                                                                >  5:1218887/1‑76 (MQ=255)
                                                                                                                                                                                                                                                               tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                                                                                                                                                                                                 >  1:244397/1‑72 (MQ=255)
                                                                                                                                                                                                                                                               tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACg                                                                                                                                                                                                                 <  3:244397/72‑1 (MQ=255)
                                                                                                                                                                                                                                                               tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCttt                                                                                                                                                        >  6:1256568/1‑129 (MQ=255)
                                                                                                                                                                                                                                                               tttGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATAACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTCCGACCACTATTTGCTATATCTTGTGTGGTTGAATCTTTGTTCAACTCCATCTAGTATCCCTGTATCAACAGCGAGACAACCCGACGCGTACCATCGCGCCGTATCTTCATTTTAAACGGAACTACGAATCATGCGCATTACTCTTTACACTCGTAACGATTACGTTCAATGCCACGCAACCaaa      >  4:149575/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 tGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTatatat                                                                                                                                                                         >  4:676573/1‑110 (MQ=255)
                                                                                                                                                                                                                                                                 tGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTAACATCGTAACGATTGCGTTTAGTgcctgtgcctataaca    >  4:506255/1‑262 (MQ=255)
                                                                                                                                                                                                                                                                  gACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATctct                                                                                                                                 <  5:856451/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAACTCCCCCTCCCCCGTCACGCTCTGATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCACGCCACCAAACggg  <  2:257990/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCCCCGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTg                                                                                                                                >  7:561019/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ccccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt                                                                                                                               <  6:390296/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGt                                                                                                                               >  6:914095/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttt                                                                                                                                                       >  8:483832/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                     cccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGAtt                             <  2:506259/246‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATtg                                                                                                                                                                       <  7:829203/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCtttt                                                                                                                                                       >  7:1650737/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                     cccGAAGACAGGCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCaa                                                                                                                                                  >  5:160538/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                      ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTAt                                                                                                                             <  7:171505/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTAt                                                                                                                             >  6:1582270/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           gCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACa                                                                                                                        >  6:761942/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            ccTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACag                                                                                                                       >  8:1377643/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                       
CTGGGGCAGCGACGCCAAAGGGGAAGCTGAAGCCGCACGCAGTAAAGCTCAGGCATTGCTGAAAGAAACCAGGGCGCGAATGCATGGTCGTACTCGCGTCCAGCAAGCCGCGCGCGATGCCGTTGGCTGCGCGGATTCTTTTGTTCGTGAAAGACCCTGGTGTAGCGTGGGTACAGCAGCTGCGGTAGGTATTTTTATCGGCGCACTGTTAAGCATGCGCAAATCGTAGTGCAAAAATGATAATAAATACGCGTCTTTGACCCCGAAGCCTGTCTTCGGGGTTTCTTTTTGCCTGGTGAATCACAAAAATCCCCCTACCCCGTCACGCTCATATCCAGGGTAATTTCGACCACTATTTGCTATATATTGTGTGGTTGAATCTTTTTTCAACTACATCTAGTATCTCTGTATCAACAGAGAGACAACCCGACGCGTATCATCGCGCCGTATCTTCATTTTAAACGGAAATACGAATCATGCGCATTACTATTTACACTCGTAACGATTGCGTTCAGTGCCACGCCACCAAACGGG  >  NC_000913/2800247‑2800780

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