Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,600,288 C→T 100% S128L (TCG→TTG)  gcvR → transcriptional repressor, regulatory protein accessory to GcvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,600,2880CT98.8% 562.9 / ‑2.8 170S128L (TCG→TTG) gcvRtranscriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (83/85);  total (84/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 171 total aligned reads displayed.

TCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600192‑2600383
                                                                                                |                                                                                               
 ctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTgcgc                                                                                             >  2:1915670/1‑100 (MQ=255)
 ctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTgcgc                                                                                             <  2:2928956/100‑1 (MQ=255)
  tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                            <  2:1816434/100‑1 (MQ=255)
  tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                            <  1:277467/100‑1 (MQ=255)
   gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                           >  1:3337734/1‑100 (MQ=255)
          ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                     <  1:3269837/99‑1 (MQ=255)
          ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACCTTTCAACGCGCAGCATATGCACATTGCGGAGCTGGTGTTGCGCACGCGAcc                                                                                     >  2:792945/1‑99 (MQ=255)
           gTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcacc                                                                                   <  2:603955/100‑4 (MQ=255)
           gTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctg                                                                                   >  1:790306/1‑100 (MQ=255)
           gTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctg                                                                                   <  2:2680477/100‑1 (MQ=255)
           gTCGATGTGGCAGACTCCCCGCATTTAATTGAAAGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                    >  2:2475886/1‑99 (MQ=255)
            tCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcacca                                                                                  <  2:1224325/100‑5 (MQ=255)
            tCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcacca                                                                                  <  1:3872545/100‑5 (MQ=255)
            tCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgc                                                                                  >  2:4008840/1‑100 (MQ=255)
             cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccag                                                                                 <  2:3326370/100‑6 (MQ=255)
             cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcacca                                                                                  <  1:3485769/99‑5 (MQ=255)
              gATGTGGCAGACTCCCCGCATTTAATTGCACGCTTCACAGCAGTTTTCGACGCGCATCATATGAACATTGCGGGGCTGGTGTTGCGTGCGCAAcaccagc                                                                                <  2:2032303/100‑7 (MQ=255)
              gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccagc                                                                                <  1:1041320/100‑7 (MQ=255)
              gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccagc                                                                                <  2:842002/100‑7 (MQ=255)
              gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccagc                                                                                <  2:3528031/100‑7 (MQ=255)
              gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgctg                                                                                >  1:1525440/1‑100 (MQ=255)
                tgtgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaa                                                                              <  2:3982201/100‑1 (MQ=255)
                   ggCCGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                            >  2:229664/1‑99 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  2:96750/1‑100 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  1:2948961/1‑100 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           <  2:2237379/100‑1 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  2:2443166/1‑100 (MQ=255)
                    gCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGTGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATg                                                                          <  1:345170/100‑1 (MQ=255)
                       gACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAAc                                                                       <  1:2987818/100‑1 (MQ=255)
                       gACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAAc                                                                       <  1:536784/100‑1 (MQ=255)
                           ccccGCATTTAATTGCACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCGGAAAATGAACGGGc                                                                   <  2:763531/100‑1 (MQ=255)
                           ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATAATATGAACATTGCGGAGCTGGTGTTGCGCACGCCCCGTGCTGAAAATGAACGAGc                                                                   >  2:424496/1‑100 (MQ=255)
                            cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGTTGTTGCGCACGCAACCTGCTGAAAGTGAACGGGCt                                                                  >  2:2435217/1‑100 (MQ=255)
                            cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                  >  1:2207244/1‑100 (MQ=255)
                             ccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGATCGGGCTg                                                                 <  1:2168601/100‑1 (MQ=255)
                              cGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTgc                                                                <  1:249886/100‑1 (MQ=255)
                              cGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTgc                                                                >  2:2103882/1‑100 (MQ=255)
                              cGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTgc                                                                <  1:2251599/100‑1 (MQ=255)
                               gCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTgcg                                                               >  1:273089/1‑100 (MQ=255)
                                 aTTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAACCAGTGTTGCGCACGCAACGTGGTGAACAGGAGCGGGc                                                                   >  2:3422866/1‑94 (MQ=255)
                                  tttAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAg                                                            >  1:3510936/1‑100 (MQ=255)
                                    tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                          >  1:3467562/1‑100 (MQ=255)
                                     aaTTGAACGCTCCACAGCACTTTTCGACGCGCACCCTATGAACATTGCGGAGCTGGTGTTCCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         <  2:1825400/100‑1 (MQ=255)
                                      aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                        >  2:1525853/1‑100 (MQ=255)
                                      aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                        >  2:3687961/1‑100 (MQ=255)
                                       ttGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCa                                                       >  2:3397231/1‑100 (MQ=255)
                                        tGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCat                                                      <  1:5871/100‑1 (MQ=255)
                                         gAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCata                                                     >  1:1326597/1‑100 (MQ=255)
                                            cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                  <  1:565787/100‑1 (MQ=255)
                                            cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                  >  1:1456056/1‑100 (MQ=255)
                                             gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  2:517401/100‑1 (MQ=255)
                                             gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:3895786/100‑1 (MQ=255)
                                                 cacaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATa                                             >  2:3572787/1‑100 (MQ=255)
                                                 cacaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATa                                             >  2:284286/1‑100 (MQ=255)
                                                  acaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATaa                                            >  2:2626039/1‑100 (MQ=255)
                                                   caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATa                                             <  2:3079329/98‑1 (MQ=255)
                                                   caGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCAGCGCAGTTGCATATTCAGATAAc                                           >  1:3575318/1‑100 (MQ=255)
                                                     gCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCg                                         <  1:1042536/100‑1 (MQ=255)
                                                     gCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCg                                         >  1:961097/1‑100 (MQ=255)
                                                      cACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGc                                        >  1:982615/1‑100 (MQ=255)
                                                        cTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                      <  2:3672732/100‑1 (MQ=255)
                                                        cTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                      >  1:3671498/1‑100 (MQ=255)
                                                        cTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGAGCAGTTGCATATTCAGATAACCGcc                                       >  1:32223/1‑99 (MQ=255)
                                                        cTTTTCGACGCGCATCAGATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                      <  1:1339676/100‑1 (MQ=255)
                                                        cTTCTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                      <  2:1072108/100‑1 (MQ=255)
                                                         ttttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCca                                     <  1:3777181/100‑1 (MQ=255)
                                                          tttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcac                                    <  1:1840633/100‑1 (MQ=255)
                                                             cGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGc                                 >  1:1713558/1‑100 (MQ=255)
                                                              gACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGcc                                >  1:3085232/1‑100 (MQ=255)
                                                               aCGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGccc                               <  2:122950/100‑1 (MQ=255)
                                                                 gcgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                             <  2:2720185/100‑1 (MQ=255)
                                                                   gcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCa                           >  2:919963/1‑100 (MQ=255)
                                                                     atcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATc                         <  1:618746/100‑1 (MQ=255)
                                                                      tcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATc                         <  2:1417952/99‑1 (MQ=255)
                                                                         tttGAACATTGCTGAGCTGGTGTTGCGCACGCATCCTGCTGAAAATGAACGGGCTGCGCAGTTTCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     <  2:3645274/98‑1 (MQ=255)
                                                                         tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     <  1:2282583/100‑1 (MQ=255)
                                                                         tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCAACAGCCCCGCATCTGCg                     <  2:3231954/100‑1 (MQ=255)
                                                                          atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                    <  2:383687/100‑1 (MQ=255)
                                                                           tGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGa                   >  2:871802/1‑100 (MQ=255)
                                                                           tGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGa                   >  1:159838/1‑100 (MQ=255)
                                                                            gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  <  2:571700/100‑1 (MQ=255)
                                                                            gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  >  1:2018602/1‑100 (MQ=255)
                                                                            gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  <  1:3014538/100‑1 (MQ=255)
                                                                             aaCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACg                 >  2:2959466/1‑100 (MQ=255)
                                                                              aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc                >  1:500585/1‑100 (MQ=255)
                                                                              aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc                <  2:45646/100‑1 (MQ=255)
                                                                               cATTGCGGAGCTGGTGTTGCGCACGCAACCTACTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgca               <  1:1059529/100‑1 (MQ=255)
                                                                                   gCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaa           >  1:1401586/1‑100 (MQ=255)
                                                                                   gCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaa           >  1:1691374/1‑100 (MQ=255)
                                                                                    cGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaaa          >  1:2376061/1‑100 (MQ=255)
                                                                                     ggAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaaa          <  2:1752901/99‑1 (MQ=255)
                                                                                        gCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAtt      >  1:2694761/1‑100 (MQ=255)
                                                                                         cTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCATAGCCCCGCATCTGCGGACGCAGCAAATATTg     <  1:2103895/100‑1 (MQ=255)
                                                                                         cTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTg     <  2:2326604/100‑1 (MQ=255)
                                                                                          tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    <  1:3125198/100‑1 (MQ=255)
                                                                                          tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    >  2:1721834/1‑100 (MQ=255)
                                                                                           ggTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCACATATTGa    >  2:788411/1‑99 (MQ=255)
                                                                                           ggTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAg   <  1:441107/100‑1 (MQ=255)
                                                                                            gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  1:3767906/100‑1 (MQ=255)
                                                                                            gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGATGCGCAGTTGCATATTCAAATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    <  2:3964754/98‑1 (MQ=255)
                                                                                                |                                                                                               
TCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600192‑2600383

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