Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,221,689 G→A 78.3% R101C (CGC→TGC)  bglX ← beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,221,6890GA78.3% 246.0 / 69.1 133R101C (CGC→TGC) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (16/13);  new base A (54/50);  total (70/63)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 138 total aligned reads displayed.

ATCGAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGG  >  NC_000913/2221591‑2221786
                                                                                                  |                                                                                                 
aTCGAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAggcg                                                                                                   <  2:2658602/99‑1 (MQ=255)
 tCGAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAggcggc                                                                                                 >  2:2057245/1‑100 (MQ=255)
 tCGAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAggcagc                                                                                                 <  1:2521754/100‑1 (MQ=255)
   gAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCtt                                                                                               >  2:592544/1‑100 (MQ=255)
   gAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACAACGTGAGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCtt                                                                                               <  2:3529508/100‑1 (MQ=255)
        tAATAGACGCGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAACAGGCATTTTCATGCAGCTTCATTc                                                                                          >  2:1456166/1‑100 (MQ=255)
        tAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTc                                                                                          >  2:868279/1‑100 (MQ=255)
        tAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTc                                                                                          >  1:3549627/1‑100 (MQ=255)
        tAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTc                                                                                          >  1:3549310/1‑100 (MQ=255)
         aaaaGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGAAAAGAAAAGAGGAATTTTCAGGCGGCTTAAttac                                                                                         >  2:813735/1‑98 (MQ=255)
         aaaaGACGAGGCCAGACCGAGGCTAATCGGAAAAACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTcc                                                                                         >  2:2739322/1‑100 (MQ=255)
          aaaGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCa                                                                                        >  2:1148499/1‑100 (MQ=255)
          aaaGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCa                                                                                        <  1:2964706/100‑1 (MQ=255)
           aaGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCAt                                                                                       <  2:1784685/100‑1 (MQ=255)
             gACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCa                                                                                     <  1:2870090/100‑1 (MQ=255)
               cGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCAcc                                                                                   <  2:4005776/100‑1 (MQ=255)
               cGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCAcc                                                                                   <  1:2647856/100‑1 (MQ=255)
                gAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCt                                                                                  <  1:1304901/100‑1 (MQ=255)
                 aGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGGAAAGAGGAATTTTCTGGCAGCTTAATTCCATCACCTg                                                                                 >  2:2806950/1‑100 (MQ=255)
                    ccAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATc                                                                              >  1:3288151/1‑100 (MQ=255)
                    ccAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATc                                                                              >  1:1388879/1‑100 (MQ=255)
                    ccAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATc                                                                              >  1:2066926/1‑100 (MQ=255)
                     cAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATcc                                                                             <  1:4010663/100‑1 (MQ=255)
                     cAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATcc                                                                             <  2:743379/100‑1 (MQ=255)
                         ccGAGGCTAATCGGGAACACCGTGCGCTGACCTTTCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCa                                                                         <  2:2812525/100‑1 (MQ=255)
                          ctagtcTAATCGGGAACACCGGTCGCTGCCCGTGCATCACGTCGTAAGCAAAGAAAAGAGGCATTTTCAGGCAGCTTAATTCCATCACCGGAttc                                                                             <  2:445813/90‑3 (MQ=255)
                          cGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCAt                                                                        <  2:2118625/100‑1 (MQ=255)
                           gAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATg                                                                       >  2:370134/1‑100 (MQ=255)
                            aGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATgg                                                                      <  2:1445069/100‑1 (MQ=255)
                            aGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATgg                                                                      <  2:3859081/100‑1 (MQ=255)
                             ggCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGgc                                                                     >  2:2876643/1‑100 (MQ=255)
                             ggCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGgc                                                                     >  1:3159946/1‑100 (MQ=255)
                             ggCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGgc                                                                     >  1:1699221/1‑100 (MQ=255)
                               cTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGgc                                                                     >  2:3546771/1‑98 (MQ=255)
                                 aaTCGGGTACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCgg                                                                 <  2:2235946/100‑1 (MQ=255)
                                 aaTCGGGAACACCGTGCGCTGGCCGTGCAGCACGTCGTAAGCTAAGAAATGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCgg                                                                 <  1:1694461/100‑1 (MQ=255)
                                 aaTCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCgg                                                                 <  1:3935561/100‑1 (MQ=255)
                                   tCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGat                                                               >  2:3582599/1‑100 (MQ=255)
                                    cGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGata                                                              <  1:3102560/100‑1 (MQ=255)
                                     gggAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGatat                                                             >  2:951047/1‑100 (MQ=255)
                                     gggAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAAATCCATCACCTGATCCTGCATGGCGCGGatat                                                             >  2:623220/1‑100 (MQ=255)
                                      ggAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATc                                                            >  1:138225/1‑100 (MQ=255)
                                        aaCACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCt                                                          >  1:1479205/1‑100 (MQ=255)
                                        aaCACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCt                                                          <  1:1811113/100‑1 (MQ=255)
                                        aaCACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCt                                                          <  1:3766712/100‑1 (MQ=255)
                                             cGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACgg                                                     >  2:448507/1‑100 (MQ=255)
                                             cGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACgg                                                     <  2:4088731/100‑1 (MQ=255)
                                              gTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACggg                                                    <  1:3272946/100‑1 (MQ=255)
                                              gTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACggg                                                    <  1:3064423/100‑1 (MQ=255)
                                               tGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGt                                                   >  2:228825/1‑100 (MQ=255)
                                               tGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGt                                                   <  2:163316/100‑1 (MQ=255)
                                                gcgcTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGtt                                                  <  2:618043/100‑1 (MQ=255)
                                                 cgcTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTa                                                 >  1:865768/1‑100 (MQ=255)
                                                  gcTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTAc                                                <  2:1454910/100‑1 (MQ=255)
                                                    tGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACgg                                              <  1:2596351/100‑1 (MQ=255)
                                                     gACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGgt                                             >  1:3302696/1‑100 (MQ=255)
                                                      aCCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGgtg                                            <  1:1064748/100‑1 (MQ=255)
                                                       ccGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGgtgt                                           >  2:1694877/1‑100 (MQ=255)
                                                        cGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGtt                                          <  1:1089468/100‑1 (MQ=255)
                                                         gTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTg                                         <  2:2406963/100‑1 (MQ=255)
                                                           gcagcaCGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGCCGGGTTACGGTGTTGaa                                       >  1:1646134/1‑100 (MQ=255)
                                                            cagcaCGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGaaa                                      >  1:2416860/1‑100 (MQ=255)
                                                              ggaCGTCGTAAGCGAAGAAAAGAGGAATTTTGAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAAt                                    >  1:2805910/3‑100 (MQ=255)
                                                              gcaCGTCGTAAGCAAAGAAAAGAGGAATTTTCAACCAGCTTAATTCCATCACATGATCCTGCATGGCGCGGATATCATGACGGTTTACAATGTTTaaa                                      >  2:97572/1‑98 (MQ=255)
                                                               caCGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATc                                   <  2:880742/100‑1 (MQ=255)
                                                                  gtcgtAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGcc                                >  2:1698494/1‑100 (MQ=255)
                                                                  gtcgtAAGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGcc                                >  2:305104/1‑100 (MQ=255)
                                                                   gcgtAAGCAAAGCAAAGAGAAAACTTCAGGGAGCGTAATTCCAGGACCTGATGCCGCAGGGCGCGGATATCCTGACGGGTTACGGTCTTGAAAAgcgcn                                <  2:1157688/98‑2 (MQ=255)
                                                                        ngCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAAcc                          >  2:727148/2‑100 (MQ=255)
                                                                        aGCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAAcc                          <  2:1100157/100‑1 (MQ=255)
                                                                         gCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAATATCGCCCCAACCt                         >  2:246437/1‑100 (MQ=255)
                                                                         gCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCt                         >  2:1151320/1‑100 (MQ=255)
                                                                         gCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCt                         >  1:3199350/1‑100 (MQ=255)
                                                                         gCAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCt                         >  1:3821163/1‑100 (MQ=255)
                                                                         gCAAAGAAAAGAGAAATTTTCAGGCGGCTTAATTCCATCAACTGATACTGCATGGCGCGGATATCCTGACGGGTTACGGTGCTGAAAATCGCCCCAAcc                          >  2:1604533/1‑99 (MQ=255)
                                                                          cAAAGAAAAGAGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTg                        >  2:415598/1‑100 (MQ=255)
                                                                           aaagaaaagaGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTg                        >  2:4008890/1‑99 (MQ=255)
                                                                           aaagaaaagaGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGa                       >  2:3347383/1‑100 (MQ=255)
                                                                            aagaaaagaGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGAc                      <  2:1235998/100‑1 (MQ=255)
                                                                            aagaaaagaGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGAc                      >  1:2976238/1‑100 (MQ=255)
                                                                             agaaaagaGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGAcc                     >  1:1102693/1‑100 (MQ=255)
                                                                               aaaagaGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGt                   >  1:1741912/1‑100 (MQ=255)
                                                                                aaagaGGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTc                  >  1:54038/1‑100 (MQ=255)
                                                                                   gagGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCAAGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCttt               <  2:2922982/100‑1 (MQ=255)
                                                                                    agGAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTg              <  1:3347409/100‑1 (MQ=255)
                                                                                     ggAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGTATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGa             >  2:2572020/1‑100 (MQ=255)
                                                                                     ggAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGa             >  2:1894027/1‑100 (MQ=255)
                                                                                     ggAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGa             >  2:1318193/1‑100 (MQ=255)
                                                                                      gAATTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGat            <  1:3934893/100‑1 (MQ=255)
                                                                                       aaTTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGatc           <  1:868284/100‑1 (MQ=255)
                                                                                       aaTTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGTCCCAACCTGACCGTCTTTGatc           <  1:3111257/100‑1 (MQ=255)
                                                                                       aaTTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGatc           >  1:1467167/1‑100 (MQ=255)
                                                                                       aaTTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGatc           <  2:2690439/100‑1 (MQ=255)
                                                                                        aTTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGatca          >  2:310782/1‑100 (MQ=255)
                                                                                        aTTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGatca          <  2:3369304/100‑1 (MQ=255)
                                                                                        aTTTTCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGatca          <  2:1131983/100‑1 (MQ=255)
                                                                                           ttCAGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCAtct       >  2:4122124/1‑100 (MQ=255)
                                                                                              aGGCAGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTcgc    <  2:2356569/100‑1 (MQ=255)
                                                                                               ggcagcTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTcgcg   >  1:2759443/1‑100 (MQ=255)
                                                                                                gcggcTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCgg  >  1:2242243/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
ATCGAGGTTAAAAGACGAGGCCAGACCGAGGCTAATCGGGAACACCGTGCGCTGACCGTGCAGCACGTCGTAAGCAAAGAAAAGAGGAATTTTCAGGCGGCTTAATTCCATCACCTGATCCTGCATGGCGCGGATATCCTGACGGGTTACGGTGTTGAAAATCGCCCCAACCTGACCGTCTTTGATCATCTCGCGG  >  NC_000913/2221591‑2221786

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