Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,012,553 Δ1 bp intergenic (+160/+61) yqeA → / ← yqeB putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,5520G.100.0% 315.4 / NA 81intergenic (+159/+62)yqeA/yqeBputative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (40/41);  total (40/41)

TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGT  >  NC_000913/3012419‑3012683
                                                                                                                                     |                                                                                                                                   
tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAAtt                                                                                                                               >  2:890846/1‑139 (MQ=255)
     ccAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCg                                                                                                                          >  1:328890/1‑139 (MQ=255)
     ccAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCg                                                                                                                          >  1:204279/1‑139 (MQ=255)
        aaGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGttt                                                                                                                       >  2:1001381/1‑139 (MQ=255)
             tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:809098/139‑1 (MQ=255)
                        tGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                     <  1:231730/110‑1 (MQ=255)
                        tGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                     >  2:231730/1‑110 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  2:473509/139‑1 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:1371599/139‑1 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:1001381/139‑1 (MQ=255)
                             tttGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggt                                                                                                  <  2:1571124/139‑1 (MQ=255)
                             tttGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggt                                                                                                  >  1:1571124/1‑139 (MQ=255)
                               tgttgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:1593945/135‑1 (MQ=255)
                               tgttgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  2:1593945/1‑135 (MQ=255)
                                 ttgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtgggt                                                                                              >  2:689680/1‑139 (MQ=255)
                                  tgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTa                                                                                             >  1:787898/1‑139 (MQ=255)
                                   gATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTaa                                                                                            <  2:204279/139‑1 (MQ=255)
                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTTAGGCCGGGTGGGTAATa                                                                                          <  2:137954/139‑1 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                      >  1:468312/1‑139 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                      >  2:1091607/1‑139 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                      >  1:1428485/1‑139 (MQ=255)
                                          gATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc                                                                                     <  1:864516/139‑1 (MQ=255)
                                               ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  2:1076120/105‑1 (MQ=255)
                                               ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  >  1:1076120/1‑105 (MQ=255)
                                               ggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGcc                                                                                >  2:423368/1‑139 (MQ=255)
                                                            tATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg                                                                   <  2:1428485/139‑1 (MQ=255)
                                                              tCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                                 <  2:1423994/139‑1 (MQ=255)
                                                                   ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTc                                                            <  1:1091607/139‑1 (MQ=255)
                                                                       cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                                 <  2:1576635/130‑1 (MQ=255)
                                                                       cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                                 >  1:1576635/1‑130 (MQ=255)
                                                                        aaaaCATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGcc                                                                                                      >  1:306605/1‑92 (MQ=255)
                                                                        aaaaCATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGcc                                                                                                      <  2:306605/92‑1 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  >  2:1567579/1‑76 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:1567579/76‑1 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:1374102/76‑1 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  2:764080/1‑90 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:764080/90‑1 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACGTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcgt                                                                                                                  >  2:1374102/1‑74 (MQ=18)
                                                                               cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCttt                                                >  1:856769/1‑139 (MQ=255)
                                                                               cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCttt                                                >  2:816534/1‑139 (MQ=255)
                                                                                    aTTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                                           <  1:704791/139‑1 (MQ=255)
                                                                                          tAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                                     <  1:689680/139‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  2:263703/1‑71 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:263703/71‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCa                                                                >  1:1335034/1‑107 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCa                                                                <  2:1335034/107‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                                >  1:1624818/1‑139 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                                >  1:212432/1‑139 (MQ=255)
                                                                                                    tGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATg                           >  2:164092/1‑139 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:57801/50‑1 (MQ=17)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  >  2:57801/1‑50 (MQ=17)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtc                                                                                                                    >  2:421934/1‑48 (MQ=17)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtc                                                                                                                    <  1:421934/48‑1 (MQ=17)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                   >  2:1354280/1‑129 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                   <  1:1354280/129‑1 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         >  1:528310/1‑139 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         <  2:328890/139‑1 (MQ=255)
                                                                                                          tGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATGCCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa                     <  2:1891063/139‑1 (MQ=255)
                                                                                                           gTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc                    >  2:590075/1‑139 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAcatcat               >  1:1579392/1‑139 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:1328531/37‑1 (MQ=14)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  >  2:1328531/1‑37 (MQ=14)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa            <  1:164092/139‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa            <  2:468312/139‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa            <  1:590075/139‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa            <  2:856769/139‑1 (MQ=255)
                                                                                                                    aTAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAAc           <  2:787898/139‑1 (MQ=255)
                                                                                                                    aTAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAAc           <  2:528310/139‑1 (MQ=255)
                                                                                                                             gcgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGt  >  1:441494/1‑139 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                               >  2:170129/1‑45 (MQ=25)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                               <  1:170129/45‑1 (MQ=25)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:1487855/40‑1 (MQ=25)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  1:1269648/1‑40 (MQ=25)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  2:1487855/1‑40 (MQ=21)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  2:1269648/40‑1 (MQ=25)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCa                                                       >  2:1082017/1‑85 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCa                                                       <  1:1082017/85‑1 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTc                                      <  2:900048/102‑1 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTc                                      >  1:900048/1‑102 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTTACGCCCTGAt                            >  1:668412/1‑112 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGAt                            <  2:668412/112‑1 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCAtt              <  2:1604861/126‑1 (MQ=255)
                                                                                                                              cgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCAtt              >  1:1604861/1‑126 (MQ=255)
                                                                                                                                     |                                                                                                                                   
TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGT  >  NC_000913/3012419‑3012683

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