Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,756 T→C T196A (ACT→GCT)  gatZ ← putative tagatose‑6‑phosphate kinase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,7560TC99.3% 460.3 / ‑4.6 145T196A (ACT→GCT) gatZputative tagatose‑6‑phosphate kinase 1
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (64/80);  total (64/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

Only 100 of 147 total aligned reads displayed.

CCACACCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGTATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACAT  >  NC_000913/2173666‑2173845
                                                                                          |                                                                                          
ccACACCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCAt                                                                                          <  2:2992360/93‑1 (MQ=255)
 caccCCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCAtt                                                                                         <  2:3496436/93‑1 (MQ=255)
 cacaCCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGTAtt                                                                                         <  1:1446312/93‑1 (MQ=255)
 cacaCCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCAtt                                                                                         <  1:2894862/93‑1 (MQ=255)
 cacaCCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCAtt                                                                                         <  1:2937618/93‑1 (MQ=255)
 cacaCCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCAtt                                                                                         >  2:2422353/1‑93 (MQ=255)
 cacaCCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCAtt                                                                                         <  2:1962177/93‑1 (MQ=255)
  ccccccGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTg                                                                                        <  2:584055/90‑1 (MQ=255)
   cacCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTgg                                                                                       <  2:165718/93‑1 (MQ=255)
     cccGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTggcg                                                                                     <  1:683945/93‑1 (MQ=255)
      ccGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTggcgg                                                                                    >  1:3163267/1‑93 (MQ=255)
       cGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTggcggc                                                                                   <  1:2272446/93‑1 (MQ=255)
       cGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTggcggc                                                                                   >  2:2560080/1‑93 (MQ=255)
         gCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGTGGCAt                                                                                 <  1:4483611/93‑1 (MQ=255)
         gCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCAt                                                                                 >  1:2006996/1‑93 (MQ=255)
         gCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCAt                                                                                 <  1:497501/93‑1 (MQ=255)
          cTGCGCCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATc                                                                                <  2:988945/93‑1 (MQ=255)
          cTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATc                                                                                <  1:1066580/93‑1 (MQ=255)
           tGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCt                                                                              <  2:3933233/93‑1 (MQ=255)
           tGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCt                                                                              >  2:2857813/1‑93 (MQ=255)
            gCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCtt                                                                             >  2:1512885/1‑93 (MQ=255)
             caccacGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTc                                                                            <  2:1080902/93‑1 (MQ=255)
              accacGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCa                                                                           >  1:890368/1‑93 (MQ=255)
               ccacGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCaa                                                                          >  2:2608878/1‑93 (MQ=255)
                cacGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAAc                                                                         <  1:4337445/93‑1 (MQ=255)
                 ccGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACa                                                                        <  2:2648683/92‑1 (MQ=255)
                 acGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACa                                                                        <  1:1226831/93‑1 (MQ=255)
                  cGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACAt                                                                       <  1:4020198/93‑1 (MQ=255)
                  cGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACAt                                                                       >  2:1521197/1‑93 (MQ=255)
                   gATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATg                                                                      <  1:525804/93‑1 (MQ=255)
                   gATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATg                                                                      <  2:54232/93‑1 (MQ=255)
                    agggCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATgg                                                                     >  2:3579740/3‑93 (MQ=255)
                    aTGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATgg                                                                     <  1:434731/93‑1 (MQ=255)
                     tGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATggg                                                                    <  2:974656/93‑1 (MQ=255)
                      ggCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGt                                                                   >  2:2681331/1‑93 (MQ=255)
                         aaTCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGAt                                                                <  2:2461784/93‑1 (MQ=255)
                         aaTCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGAt                                                                <  1:3947506/93‑1 (MQ=255)
                         aaTCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGAt                                                                >  2:4290119/1‑93 (MQ=255)
                           tCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGAtgt                                                              >  2:2569181/1‑93 (MQ=255)
                             acacGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATgtgt                                                            <  1:3765904/93‑1 (MQ=255)
                             acacGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATgtgt                                                            >  1:954420/1‑93 (MQ=255)
                             acacGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATgtgt                                                            >  2:663033/1‑93 (MQ=255)
                               acGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTac                                                          >  1:791528/1‑93 (MQ=255)
                                cGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTact                                                         <  1:2835553/93‑1 (MQ=255)
                                   gtTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTactgac                                                      <  2:777732/93‑1 (MQ=255)
                                    ttAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTactgact                                                     >  1:2799062/1‑93 (MQ=255)
                                       aCGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGaa                                                  >  2:1591374/1‑93 (MQ=255)
                                       aCGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGaa                                                  <  1:527106/93‑1 (MQ=255)
                                       aCGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGaa                                                  >  2:3964345/1‑93 (MQ=255)
                                       aCGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGGACTGACTGaa                                                  >  2:1408131/1‑93 (MQ=255)
                                        cGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAAt                                                 <  2:451739/93‑1 (MQ=255)
                                        cGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAAt                                                 <  2:3290/93‑1 (MQ=255)
                                         gCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATg                                                >  2:1586186/1‑93 (MQ=255)
                                          ccTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATgg                                               >  2:307577/1‑93 (MQ=255)
                                          ccTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATgg                                               <  1:1745508/93‑1 (MQ=255)
                                          ccTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATgg                                               <  1:2971274/93‑1 (MQ=255)
                                          ccTCTGTCAGCCCACGGCAATAAAAGACTTTTGATGCGGACCTAAAACAATGGCGGCATCTTTCAACTGGGTGATTTGTAATGGATTAATggg                                              >  2:372340/1‑92 (MQ=255)
                                           ctctGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGgc                                              >  1:3474274/1‑93 (MQ=255)
                                             ttGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTGCAACATGGGTGATGTGTACTGACTGAATGgcgc                                            <  2:3502180/92‑1 (MQ=255)
                                               gTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTg                                          >  1:188988/1‑93 (MQ=255)
                                                tCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTgg                                         <  1:858697/93‑1 (MQ=255)
                                                  aGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGcc                                       <  2:3719106/93‑1 (MQ=255)
                                                   gCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCt                                      >  1:2923856/1‑93 (MQ=255)
                                                   gCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTAATGACTGAATGGCGCTGGGCTc                                     >  1:2491825/1‑93 (MQ=255)
                                                      cACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCAc                                   <  1:4573445/93‑1 (MQ=255)
                                                       aCGGGCAATAAAGGCCTTTTGATGCGTACGTAACGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCAcc                                  <  1:1076208/93‑1 (MQ=255)
                                                       aCGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCAcc                                  <  2:3284778/93‑1 (MQ=255)
                                                            cAATAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccg                             >  1:877314/1‑93 (MQ=255)
                                                             aaTAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccgg                            >  1:4011297/1‑93 (MQ=255)
                                                             aaTAAAGGCCTTTTGATGCGTACGTAAAGCATTGGAGGCATCTTCAAAATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccgg                            >  2:4599626/1‑93 (MQ=255)
                                                              aTAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccgga                           >  1:3191691/1‑93 (MQ=255)
                                                               tAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaa                          >  1:961457/1‑93 (MQ=255)
                                                               tAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaa                          <  2:3544324/93‑1 (MQ=255)
                                                               tAAAGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaa                          >  2:3701675/1‑93 (MQ=255)
                                                                aaaGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaac                         >  1:3979937/1‑93 (MQ=255)
                                                                 aaGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaacc                        >  1:1399854/1‑93 (MQ=255)
                                                                 aaGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaacc                        <  2:232864/93‑1 (MQ=255)
                                                                  aGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGGCTGAATGGCGCTGGCCTCACCGCccggaaccg                       >  1:1833491/1‑93 (MQ=255)
                                                                  aGGCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaaccg                       >  1:2269426/1‑93 (MQ=255)
                                                                   ggCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaaccgg                      >  2:659399/1‑93 (MQ=255)
                                                                    gCCTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaaccgga                     >  1:1035783/1‑93 (MQ=255)
                                                                     ccTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGTCTGAATGGCGCTGGCCTCACAGCccggaaccggac                    >  2:4205843/1‑92 (MQ=255)
                                                                     ccTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCccggaaccggaa                    >  1:1532529/1‑93 (MQ=255)
                                                                      cTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCcggaaccggaac                   <  1:2765801/93‑1 (MQ=255)
                                                                      cTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCcggaaccggaac                   >  1:3464065/1‑93 (MQ=255)
                                                                      cTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCcggaaccggaac                   <  1:3740814/93‑1 (MQ=255)
                                                                      cTTTTGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCcggaaccggaac                   <  2:1995407/93‑1 (MQ=255)
                                                                       ttttGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACt                  >  1:4454714/1‑93 (MQ=255)
                                                                       ttttGATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACt                  >  2:778432/1‑93 (MQ=255)
                                                                           gATGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCgg              >  1:3056886/1‑93 (MQ=255)
                                                                            aTGCGTACGTAAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGt             <  1:2110390/93‑1 (MQ=255)
                                                                                 tacgtaAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCaa        <  2:1103310/93‑1 (MQ=255)
                                                                                 tacgtaAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCaa        <  1:2819350/93‑1 (MQ=255)
                                                                                    gtaAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGa     >  1:4299737/1‑93 (MQ=255)
                                                                                    gtaAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGa     >  1:889260/1‑93 (MQ=255)
                                                                                     taAAGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGAc    >  1:678926/1‑93 (MQ=255)
                                                                                      aaaGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACa   <  2:1206463/93‑1 (MQ=255)
                                                                                       aaGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACAt  <  1:1168626/93‑1 (MQ=255)
                                                                                       aaGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACAt  <  1:3083080/93‑1 (MQ=255)
                                                                                       aaGCATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACAt  <  2:178310/93‑1 (MQ=255)
                                                                                          |                                                                                          
CCACACCCGGCTGCACCACGATGGCAATCACACGTGTTAACGCCTCTGTCAGCCCACGGGCAATAAAGGCCTTTTGATGCGTACGTAAAGTATTGGCGGCATCTTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACAT  >  NC_000913/2173666‑2173845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: