Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,918,017 C→T 5.4% intergenic (+204/+28) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,918,0170CT 5.4% 724.7 / 18.5 218intergenic (+204/+28)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/D‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base C (103/102);  new base T (4/8);  total (108/110)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-02

Only 100 of 222 total aligned reads displayed.

TGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAA  >  NC_000913/2917752‑2918276
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                   
tgctgGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATagag                                                                                                                                                                                                                                                             <  1:67032/274‑1 (MQ=255)
                          aaTGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                                                                                                                                  <  2:603473/275‑1 (MQ=255)
                                                                        gTCGAGTGGCAATCGGATTTTTATGACGCTGGTTTCAGATTGGTGACAAGGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTc                                                                                                                                                                                                                                              <  2:305523/217‑1 (MQ=255)
                                                                           gggTAGAATTCGGCTTTTTATGAGCCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                                                                                                                                                 <  4:570207/273‑1 (MQ=255)
                                                                                                                            gCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATaga                                                                                                                                                                                                                                                              >  7:36812/1‑149 (MQ=255)
                                                                                                                            gCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATaga                                                                                                                                                                                                                                                              <  7:449745/149‑1 (MQ=255)
                                                                                                                            gCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATaga                                                                                                                                                                                                                                                              >  5:736300/1‑149 (MQ=255)
                                                                                                                                gtgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCg                                                                                                                                                                                                                                                          >  7:1050490/1‑149 (MQ=255)
                                                                                                                                  gtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggt                                                                                                                                                                                                                                                        <  5:1078528/149‑1 (MQ=255)
                                                                                                                                   tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGGTGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggtt                                                                                                                                                                                                                                                       >  8:331108/1‑149 (MQ=38)
                                                                                                                                   tttATGCCGGAGGCGGCGTTAGCGCCTTATCCAGCCTACAAAACCAATTAGATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtt                                                                                                                                                                                                                                                       >  8:926038/1‑149 (MQ=255)
                                                                                                                                     tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                     >  7:678281/1‑149 (MQ=255)
                                                                                                                                      aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCTATCTGGAGCCggtttga                                                                                                                                                                                                                                                    <  8:1069587/149‑3 (MQ=39)
                                                                                                                                       ggCCGGATAAGGCGTTAACGCCGCATCCGACCTACAAAACCAATTAAATGCAATGACTCAAAAAATGATGTATACCTGATACATAGCACAAAAGGCTATTTTGCGATTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggctttgg                                                                                                                                                                                                                                                   >  8:1060468/2‑149 (MQ=255)
                                                                                                                                        gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttggt                                                                                                                                                                                                                                                  <  5:237631/149‑1 (MQ=255)
                                                                                                                                          cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                                                                                                                                                >  7:194663/1‑149 (MQ=255)
                                                                                                                                           ggATGCGGCGTTAACGCCTTATCCGGCCTACAAGACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                                                                                                                                               >  8:627069/1‑149 (MQ=255)
                                                                                                                                              tGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             <  7:848781/148‑1 (MQ=255)
                                                                                                                                              tGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCt                                                                                                                                                                                                                                             >  8:848781/1‑148 (MQ=255)
                                                                                                                                                 ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCt                                                                                                                                                                                                                                         >  6:526131/1‑149 (MQ=255)
                                                                                                                                                    gTTAACGCCCTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGCTAGAGCCGGTTTTGGTTTTCTGTCTTat                                                                                                                                                                                                                                      >  8:694439/1‑148 (MQ=255)
                                                                                                                                                      tAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTg                                                                                                                                                                                                                                            <  2:716789/141‑1 (MQ=255)
                                                                                                                                                         cGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAc                                                                                                                                                                                                                                 <  5:584246/149‑1 (MQ=255)
                                                                                                                                                         cGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAc                                                                                                                                                                                                                                 <  6:6235/149‑1 (MQ=255)
                                                                                                                                                          gCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAcc                                                                                                                                                                                                                                >  5:535253/1‑149 (MQ=255)
                                                                                                                                                           ccTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCa                                                                                                                                                                                                                               >  6:40818/1‑149 (MQ=255)
                                                                                                                                                            cTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCggtttt                                                                                                                                                                                                                                                     >  7:288483/1‑126 (MQ=255)
                                                                                                                                                             ttATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCtt                                                                                                                                                                                                                                        >  1:45266/1‑138 (MQ=255)
                                                                                                                                                                 ccGGCCTANNAAACCAATTAAATTCAATGAATTNAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACg                                                                                                                                                                                                                                    <  3:45266/138‑1 (MQ=255)
                                                                                                                                                                  cGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAAACCGGCTCTATCTGGAGCCggttt                                                                                                                                                                                                                                                      <  6:567037/120‑1 (MQ=255)
                                                                                                                                                                  cGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                                                                                                                                        <  5:710757/149‑1 (MQ=255)
                                                                                                                                                                   ggCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACggg                                                                                                                                                                                                                       >  5:259805/1‑149 (MQ=255)
                                                                                                                                                                      cTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccag                                                                                                                                                                                                  >  4:707735/1‑167 (MQ=255)
                                                                                                                                                                        aCAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCtt                                                                                                                                                                                                                  >  7:136246/1‑149 (MQ=255)
                                                                                                                                                                          aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGt                                                                                                                                                                                                                >  8:793497/1‑149 (MQ=255)
                                                                                                                                                                             aCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                                                                                                                             <  6:667086/149‑1 (MQ=255)
                                                                                                                                                                              ccAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATc                                                                                                                                                                                                            <  6:736300/149‑1 (MQ=255)
                                                                                                                                                                                 aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGaa                                                                                                                                                                                                         <  7:1062485/149‑1 (MQ=255)
                                                                                                                                                                                 aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                  >  2:446508/1‑204 (MQ=255)
                                                                                                                                                                                  ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        >  7:575957/1‑149 (MQ=255)
                                                                                                                                                                                  ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        <  6:862040/149‑1 (MQ=255)
                                                                                                                                                                                  ttAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        <  6:464984/149‑1 (MQ=255)
                                                                                                                                                                                     aaTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCgt                                                                                                                                              <  4:446508/204‑1 (MQ=255)
                                                                                                                                                                                       ttCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTcca                                                                                                                                                                                                   >  7:143475/1‑149 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGt                                                                                                                                                                                                                                           >  4:197570/1‑106 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                                                                                                                                    <  6:635113/129‑1 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                                                                                                                                             <  6:109203/136‑1 (MQ=255)
                                                                                                                                                                                          aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAc                                                                                                                                                                                                        <  8:454906/141‑1 (MQ=255)
                                                                                                                                                                                           aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                                                                                                                               >  5:507735/1‑149 (MQ=255)
                                                                                                                                                                                           aTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcca                                                                                                                                                                                               <  7:780486/149‑1 (MQ=255)
                                                                                                                                                                                            tGAATTAAAAAATGATGTAGACCTGGTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                                                                                                                              <  8:269644/149‑1 (MQ=255)
                                                                                                                                                                                              aaTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTa                                                                                                                                                                                                                                       <  2:197570/106‑1 (MQ=255)
                                                                                                                                                                                               aTTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGaa                                                                                                                                                                                           <  5:905248/149‑1 (MQ=255)
                                                                                                                                                                                                ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCAGCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                                                                                                                                          >  6:946112/1‑149 (MQ=255)
                                                                                                                                                                                                 tAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         <  6:91091/149‑1 (MQ=255)
                                                                                                                                                                                                                  cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                                                                                                                                        <  6:1021925/149‑1 (MQ=255)
                                                                                                                                                                                                                  cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                                                                                                                                        >  5:552384/1‑149 (MQ=255)
                                                                                                                                                                                                                    gATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcatcg                                                                                                                                                                      >  5:33924/1‑149 (MQ=255)
                                                                                                                                                                                                                     aTACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCCCTTGATc                                                                                                         >  1:192009/1‑209 (MQ=255)
                                                                                                                                                                                                                       aCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGc                                                                                                                                                                                                                            <  5:586070/92‑1 (MQ=255)
                                                                                                                                                                                                                          tAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGt                                                                                                                                                                <  5:418805/149‑1 (MQ=255)
                                                                                                                                                                                                                            gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCa                                                                                                                                                              >  5:858342/1‑149 (MQ=255)
                                                                                                                                                                                                                            gCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGc                                                                                                                                                  >  2:630539/1‑161 (MQ=255)
                                                                                                                                                                                                                             cacaACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGttt                                                                           >  4:100627/1‑231 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATANTGCATTCCCATCGCATCGTCACCNCCGCNAAGCCCCTGTTTCTGATACAGCTCATTGGCTTTCCTCACTTGATCCataaccc                                                                                                 >  4:653163/1‑202 (MQ=255)
                                                                                                                                                                                                                                aacaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                        >  1:152238/1‑167 (MQ=255)
                                                                                                                                                                                                                                 acaGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCt                                                                       <  2:100627/231‑1 (MQ=255)
                                                                                                                                                                                                                                  caGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                                                                                                                                      >  1:108848/1‑71 (MQ=255)
                                                                                                                                                                                                                                  caGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACt                                                                                                              >  4:308439/1‑191 (MQ=255)
                                                                                                                                                                                                                                      cTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCa                                                                                                                                                                                                                                  <  3:108848/71‑1 (MQ=255)
                                                                                                                                                                                                                                      cTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCCTTCCCATCGCATCGTCACGCGCGCCaa                                                                                                                                                    <  8:471408/149‑1 (MQ=255)
                                                                                                                                                                                                                                      cTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGAt                                                                                                          <  2:308439/191‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTTTCCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGGGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGgcgc                    <  2:530188/275‑1 (MQ=255)
                                                                                                                                                                                                                                         ttttGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGcc                                                                                                                                                 >  6:809436/1‑149 (MQ=255)
                                                                                                                                                                                                                                         ttttGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCAc                                                                                                               <  4:727097/183‑1 (MQ=255)
                                                                                                                                                                                                                                               gTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGtt                                                                                                                                           >  6:159648/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  tAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                        <  6:974156/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCt                                                                                                                                        <  5:809436/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTg                                                                                                                                       >  7:316227/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                                                                                                                                      <  8:178773/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         <  7:437194/96‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATc                                                                                                                                                                                         >  8:437194/1‑96 (MQ=255)
                                                                                                                                                                                                                                                       tCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                                                                                                                              <  5:374415/90‑1 (MQ=255)
                                                                                                                                                                                                                                                       tCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagccag                                                                                                                                                                                              >  6:374415/1‑90 (MQ=255)
                                                                                                                                                                                                                                                       tCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAc                                                                                                                                   >  8:28982/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       tCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAc                                                                                                                                   >  5:953833/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        cAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACa                                                                                                                                  >  6:585943/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             tc‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                                                                                                                                          <  5:776899/104‑1 (MQ=255)
                                                                                                                                                                                                                                                             tc‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                                                                                                                                          <  6:714839/104‑1 (MQ=255)
                                                                                                                                                                                                                                                             tc‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                                                                                                                                          >  6:776899/1‑104 (MQ=255)
                                                                                                                                                                                                                                                             tc‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                                                                                                                             >  8:148108/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              caaaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCa                                                                                                                             >  8:729819/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATgg                                                                                                                          <  7:926038/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATgg                                                                                                                          >  8:881149/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACt                                                                                                                                                                                  >  3:174442/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                      gCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCt                                                                                                                       >  4:206460/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                       cTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtc                                                                                                                    <  5:677546/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   <  5:40818/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ccAGATGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtcat                                                                                                                  <  6:473763/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTtca                                                                                                                   <  2:206460/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                   
TGCTGGAACTGTTGGACGAGATGGATAATGTCGCGCGCGAAGCCAGCAAAATTCTCGGGTAATGGATGTCTTGTCGAGTAGCAATCGGCTTTTTATGACGCTGGTTTCAGATTGGTGACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTC‑AAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTTTCATCACTTGATCCATATCGATTTCTACACCCAGCCCCGGTTTTTCTGGCACCTGTACCAGCCCGCCTTTGATCTCAAACGGTTCTTTGGTCAGGCGCTGATTGCCTTCCTGCCAA  >  NC_000913/2917752‑2918276

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