Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,971,075 2 bp→GC coding (942‑943/1068 nt) rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,0750CG99.8% 1294.5 / NA 400C314W (TGC→TGGrffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (195/204);  minor base T (0/1);  total (195/205)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,971,0760GC99.8% 1292.6 / NA 400V315L (GTG→CTG) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (195/204);  minor base T (0/1);  total (195/205)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 410 total aligned reads displayed.

TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162
                                                                                       ||                                                                                      
tGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGcc                                                                                    >  1:4135441/1‑93 (MQ=255)
  gccTGATCGCCTTTGTCGCTGACCGTCCGGGTCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGc                                                                                  <  1:6023406/92‑1 (MQ=255)
  aCCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGc                                                                                  <  2:6202560/93‑1 (MQ=255)
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 <  1:2396242/93‑1 (MQ=255)
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 >  1:2224002/1‑93 (MQ=255)
   ccTGATCACCTTTGTCGCTGAACGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 >  2:6449243/1‑93 (MQ=255)
    cTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAg                                                                                >  2:2062697/1‑93 (MQ=255)
    cTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAg                                                                                >  1:8633788/1‑93 (MQ=255)
      gATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGa                                                                              >  2:7997969/1‑93 (MQ=255)
        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGaaa                                                                            >  2:3977202/1‑93 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           >  2:1238611/1‑93 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           <  2:5124523/93‑1 (MQ=255)
         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           <  1:8120448/93‑1 (MQ=255)
          aCCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAcc                                                                          >  2:5360238/1‑93 (MQ=255)
          aCCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAcc                                                                          <  1:3850575/93‑1 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  2:6633686/93‑1 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  1:9981773/93‑1 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  1:1360191/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCTTGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  1:5289305/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  1:6476708/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  2:1642904/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  2:3896217/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTACCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        >  1:763645/1‑93 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  1:4586551/1‑93 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  2:6270161/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  1:7235169/1‑93 (MQ=255)
              ttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTg                                                                      >  1:3010069/1‑93 (MQ=255)
               tGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGa                                                                     >  1:3549410/1‑93 (MQ=255)
                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   <  1:1120/93‑1 (MQ=255)
                  cGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAg                                                                  <  2:5421198/93‑1 (MQ=255)
                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 >  1:8379417/1‑93 (MQ=255)
                    cTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                >  2:1710128/1‑93 (MQ=255)
                    cTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                <  1:6759266/93‑1 (MQ=255)
                    cTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                <  1:9768952/93‑1 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  1:5032497/1‑93 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  1:2715600/1‑93 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  2:6332791/1‑93 (MQ=255)
                       gccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:9872859/92‑1 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:9028264/93‑1 (MQ=255)
                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            <  2:5658036/93‑1 (MQ=255)
                         cGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATg                                                           <  1:905907/93‑1 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  2:2485257/93‑1 (MQ=255)
                           tCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCg                                                         >  2:720718/1‑93 (MQ=255)
                           tCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCg                                                         >  1:6259292/1‑93 (MQ=255)
                             cGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTa                                                       <  1:4448792/93‑1 (MQ=255)
                              ggggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaa                                                      >  2:9535869/1‑93 (MQ=255)
                               gggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaa                                                     <  1:4761761/93‑1 (MQ=255)
                                ggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaaa                                                    >  2:8874450/1‑93 (MQ=255)
                                ggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaaa                                                    >  2:9750808/1‑93 (MQ=255)
                                   aTGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGTCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACgg                                                 <  1:3246367/93‑1 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAACCCTTTGAAAGTGGAATGCGTAAAACGGt                                                <  2:7583257/93‑1 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:8582827/1‑93 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:302475/1‑93 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:8386060/1‑93 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:5988767/1‑93 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:492899/1‑93 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:7422083/1‑93 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGtt                                               <  1:3159091/93‑1 (MQ=255)
                                      aTCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTc                                              >  2:1421486/1‑93 (MQ=255)
                                      aTCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTc                                              <  2:9205217/93‑1 (MQ=255)
                                        cTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAg                                            <  2:1565671/93‑1 (MQ=255)
                                        cTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAg                                            >  1:9461537/1‑93 (MQ=255)
                                        cTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGCGGAATGCGTAAAACGGTTCAg                                            >  1:9230878/1‑93 (MQ=255)
                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAACACGGTTTAGTgga                                        >  1:6886054/1‑92 (MQ=255)
                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGt                                        >  1:202031/1‑93 (MQ=255)
                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      <  1:1045105/93‑1 (MQ=255)
                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      <  2:8481126/93‑1 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     >  2:2934617/1‑93 (MQ=255)
                                                  ccATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTgg                                  <  1:1227098/93‑1 (MQ=255)
                                                  ccATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTgg                                  <  2:830190/93‑1 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  1:6981078/1‑93 (MQ=255)
                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  2:3261701/1‑93 (MQ=255)
                                                     ttGATGCTTCGAAAATTGACCGTGAACTTGGATGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTa                               >  2:6388006/1‑93 (MQ=255)
                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             >  1:8881715/1‑93 (MQ=255)
                                                        aTGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATg                            <  2:9071935/93‑1 (MQ=255)
                                                          gCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaa                          <  1:8917595/93‑1 (MQ=255)
                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        <  2:3469356/93‑1 (MQ=255)
                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       >  1:3926374/1‑93 (MQ=255)
                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       <  2:1061769/93‑1 (MQ=255)
                                                              cGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCt                      <  2:4390519/93‑1 (MQ=255)
                                                                aaaaGTGCCCGTGAACTTGGCGGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    <  2:2295756/93‑1 (MQ=255)
                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  1:7868682/1‑93 (MQ=255)
                                                                    ttGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGa                <  2:3696961/93‑1 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               >  1:1756498/1‑93 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               >  2:2691627/1‑93 (MQ=255)
                                                                       cccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGc             <  1:8271776/93‑1 (MQ=255)
                                                                        ccGTGAGCTTTGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGAAACAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  1:5591914/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            >  1:4293454/1‑93 (MQ=255)
                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAg           <  2:4739934/93‑1 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  1:2910882/93‑1 (MQ=255)
                                                                           tGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGt         >  1:4977699/1‑93 (MQ=255)
                                                                           tGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGt         <  1:7084772/93‑1 (MQ=255)
                                                                            gAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTg        >  1:7475849/1‑93 (MQ=255)
                                                                              aCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCa      <  2:4578335/93‑1 (MQ=255)
                                                                               cTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAg     <  1:1456186/93‑1 (MQ=255)
                                                                               cTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAg     >  1:1869320/1‑93 (MQ=255)
                                                                               cTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAg     >  2:5799645/1‑93 (MQ=255)
                                                                               cTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAg     <  2:4159912/93‑1 (MQ=255)
                                                                                 tggctggcTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGa   <  1:5670360/93‑1 (MQ=255)
                                                                                  ggctggctGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAc  >  2:6828875/1‑93 (MQ=255)
                                                                                       ||                                                                                      
TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162

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