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% 330.8 / NA 82intergenic (+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 . (41/41);  total (41/41)

TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAA  >  NC_000913/3012419‑3012673
                                                                                                                                     |                                                                                                                         
tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:716468/1‑139 (MQ=255)
tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:509199/1‑139 (MQ=255)
tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAtt                                                                                                                     >  1:1589504/1‑139 (MQ=255)
 ccGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAAttt                                                                                                                    >  2:1150934/1‑139 (MQ=255)
             tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTATGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        <  2:1043325/139‑1 (MQ=255)
               tCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcatca                                                                                                      >  1:1649431/1‑139 (MQ=255)
                        tGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                             >  1:805912/1‑139 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:716468/139‑1 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:1589504/139‑1 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  1:1275751/139‑1 (MQ=255)
                             tttGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggt                                                                                        <  2:724521/139‑1 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAGTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                           >  2:1299425/1‑93 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                           <  1:1299425/93‑1 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                             <  2:1656735/122‑1 (MQ=255)
                                         ggATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                             >  1:1656735/1‑122 (MQ=255)
                                           aTGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                     >  2:1134985/1‑128 (MQ=255)
                                           aTGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                     <  1:1134985/128‑1 (MQ=255)
                                             gcggcgCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATTCgg                                                                        <  1:1150934/139‑1 (MQ=255)
                                                gcgcGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        >  1:866746/1‑104 (MQ=255)
                                                gcgcGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        <  2:866746/104‑1 (MQ=255)
                                                gcgcGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          >  2:1302603/1‑118 (MQ=255)
                                                gcgcGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  1:1302603/118‑1 (MQ=255)
                                                 cgcgAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCaggtggg                                                                                     >  2:1555521/1‑122 (MQ=255)
                                                 cgcgAGCGCCTTATCCCTCCTACAAAACATCGCAAATTCAAAAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCaggtggg                                                                                     <  1:1555521/122‑1 (MQ=21)
                                                    gAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAtttt                                                                 >  1:1633549/1‑139 (MQ=255)
                                                       cgcCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt                                                              >  2:1380036/1‑139 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                     >  2:758823/1‑105 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                     <  1:758823/105‑1 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                                                                 >  1:673970/1‑109 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                                                                 <  2:673970/109‑1 (MQ=255)
                                                                   ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTc                                                  <  2:1184648/139‑1 (MQ=255)
                                                                   ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTc                                                  >  1:1184648/1‑139 (MQ=255)
                                                                    cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTcc                                                 >  1:812623/1‑139 (MQ=255)
                                                                      aCAAAACATCGCAAATTCAATAGATTGCAATGCGCTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACtt                                                   <  2:1080606/135‑1 (MQ=255)
                                                                      aCAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACtt                                                   >  1:1080606/1‑135 (MQ=255)
                                                                       cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGGATCA‑GTAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTCAGCCACTTCCAGc                                              >  2:1123605/1‑139 (MQ=38)
                                                                       cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGc                                              <  2:509199/139‑1 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                                                                        >  2:1608866/1‑108 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                                                                        <  1:1608866/108‑1 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtc                                         >  1:272685/1‑139 (MQ=255)
                                                                             aTCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtct                                        >  1:123682/1‑139 (MQ=255)
                                                                              tCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtctt                                       <  2:1649431/139‑1 (MQ=255)
                                                                               cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:1354323/87‑1 (MQ=255)
                                                                               cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          >  1:1354323/1‑87 (MQ=255)
                                                                                gCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCtttt                                     <  2:123682/139‑1 (MQ=255)
                                                                                        aaTAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCttt                             <  2:1548721/139‑1 (MQ=255)
                                                                                         aTAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTc                            <  2:272685/139‑1 (MQ=255)
                                                                                          tAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTc                                  >  2:1345265/1‑132 (MQ=255)
                                                                                          tAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTc                                  <  1:1345265/132‑1 (MQ=255)
                                                                                            gATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                         <  2:1633549/139‑1 (MQ=255)
                                                                                             aTTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                         >  1:455068/1‑138 (MQ=255)
                                                                                             aTTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                         <  2:455068/138‑1 (MQ=255)
                                                                                              ttGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTATTTTGTCGCTTTCACGcc                       <  2:941461/139‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGttt                                                                                                             <  1:726790/52‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGttt                                                                                                             >  2:726790/1‑52 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGtct                                                                                                  <  1:834068/63‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGtct                                                                                                  >  2:834068/1‑63 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                                                                >  2:1299000/1‑81 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                                                                <  1:1299000/81‑1 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          >  1:1400467/1‑64 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:1400467/64‑1 (MQ=255)
                                                                                                              agCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa           >  2:243939/2‑135 (MQ=255)
                                                                                                              agCCTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCaa           <  1:243939/134‑1 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGTACTTCTTTTGTCGCTTTCACGCCCt                     <  1:754346/123‑1 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGTACTTCTTTTGTCGCTTTCACGCCCt                     >  2:754346/1‑123 (MQ=255)
                                                                                                                 cTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  2:1321642/53‑1 (MQ=255)
                                                                                                                 cTGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          >  1:1321642/1‑53 (MQ=255)
                                                                                                                  tGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        >  2:1029041/1‑38 (MQ=14)
                                                                                                                  tGATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGtcat                                                                                                        <  1:1029041/38‑1 (MQ=14)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          >  2:1230915/1‑51 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                          <  1:1230915/51‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGt                                   <  1:1016863/106‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGt                                   >  2:1016863/1‑106 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc          <  2:1163234/131‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc          >  1:1163234/1‑131 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa  >  1:1050440/1‑139 (MQ=255)
                                                                                                                         cGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                                                                       >  2:697829/1‑64 (MQ=255)
                                                                                                                         cGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                                                                       <  1:697829/64‑1 (MQ=255)
                                                                                                                         cGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                           >  2:632113/1‑108 (MQ=255)
                                                                                                                         cGTAGCGCATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                           <  1:632113/108‑1 (MQ=255)
                                                                                                                             gcgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                                                                 >  2:1151086/1‑50 (MQ=37)
                                                                                                                             gcgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAAt                                                                                 <  1:1151086/50‑1 (MQ=37)
                                                                                                                              cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                       >  2:1095024/1‑75 (MQ=255)
                                                                                                                              cgcATCA‑GCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGcc                                                       <  1:1095024/75‑1 (MQ=255)
                                                                                                                                     |                                                                                                                         
TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAA  >  NC_000913/3012419‑3012673

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: ‑