Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,221,860 T→G K44Q (AAA→CAA)  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,8600TG99.3% 516.3 / ‑4.7 145K44Q (AAA→CAA) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (67/77);  total (68/77)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 146 total aligned reads displayed.

GACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACT  >  NC_000913/2221764‑2221956
                                                                                                |                                                                                                
gACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCtt                                                                                               <  2:1296509/100‑1 (MQ=255)
gACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCtt                                                                                               >  1:256407/1‑100 (MQ=255)
gACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCtt                                                                                               >  1:338153/1‑100 (MQ=255)
 aCCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTa                                                                                              <  1:905873/100‑1 (MQ=255)
  ccGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTaa                                                                                             >  2:3571358/1‑100 (MQ=255)
  ccGTCTTTGATCATCTCGCGGATCGCATCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTaa                                                                                             >  2:1553940/1‑100 (MQ=255)
   cGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAg                                                                                            <  2:3599350/100‑1 (MQ=255)
   cGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAg                                                                                            <  2:3668540/100‑1 (MQ=255)
    gTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGc                                                                                           >  2:661663/1‑100 (MQ=255)
    gTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGc                                                                                           <  2:1770203/100‑1 (MQ=255)
    gTCTTTGATCAGCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGc                                                                                           <  2:3905759/100‑1 (MQ=255)
       tttGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGt                                                                                        >  1:2153764/1‑100 (MQ=255)
          gATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCg                                                                                     >  2:2476643/1‑100 (MQ=255)
           atcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCgg                                                                                    <  2:1123200/100‑1 (MQ=255)
           atcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCgg                                                                                    <  2:2338719/100‑1 (MQ=255)
           atcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCgg                                                                                    <  2:1123221/100‑1 (MQ=255)
           atcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCgg                                                                                    <  1:2227469/100‑1 (MQ=255)
            tcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGt                                                                                   <  1:3696944/100‑1 (MQ=255)
            tcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGt                                                                                   <  2:2103320/100‑1 (MQ=255)
              atcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGa                                                                                 <  1:520396/100‑1 (MQ=255)
              atcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGa                                                                                 >  2:907041/1‑100 (MQ=255)
              atcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGa                                                                                 <  2:675736/100‑1 (MQ=255)
               nctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGAc                                                                                >  1:562383/2‑100 (MQ=255)
                ctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACg                                                                               >  1:3681165/1‑100 (MQ=255)
                ctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACg                                                                               >  2:1053242/1‑100 (MQ=255)
                   gcgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAAc                                                                            <  2:3113254/100‑1 (MQ=255)
                       aTCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCAt                                                                        >  1:2799559/1‑100 (MQ=255)
                       aTCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCAt                                                                        >  1:1423960/1‑100 (MQ=255)
                       aTCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCAt                                                                        <  2:305970/100‑1 (MQ=255)
                         cGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATcc                                                                      >  2:188841/1‑100 (MQ=255)
                           ccTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCcg                                                                    >  2:3268024/1‑100 (MQ=255)
                            ctctTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCcgc                                                                   <  2:1411446/100‑1 (MQ=255)
                              ctTTCGGGTTATCCGGGCCGACGCTGAATAAGCCCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCgcgc                                                                 >  2:3101917/1‑100 (MQ=255)
                               tttCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCt                                                                >  1:1309738/1‑100 (MQ=255)
                               tttCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCt                                                                <  2:1953756/100‑1 (MQ=255)
                                 tCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTc                                                              <  2:1087993/100‑1 (MQ=255)
                                   gggTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCgg                                                            >  1:2812223/1‑100 (MQ=255)
                                    ngTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                           >  2:439574/2‑100 (MQ=255)
                                    ggTTATCCGGGCCGACGCTGCTTAAGCCCAGCTGACCAATTTTCTCATCAACTGTCATTTGGTTAAGCAGATGGGTGACGAACGCATCCCGCGCTCCgag                                                           >  2:3248982/1‑98 (MQ=255)
                                    ggTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                           <  2:2493244/100‑1 (MQ=255)
                                    ggTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                           >  1:358272/1‑100 (MQ=255)
                                    ggTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                           <  2:2156421/100‑1 (MQ=255)
                                     gTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGc                                                          >  1:1741967/1‑100 (MQ=255)
                                     gTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGc                                                          <  1:2658463/100‑1 (MQ=255)
                                     gTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGc                                                          >  1:191620/1‑100 (MQ=255)
                                       tATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGt                                                        >  2:3665219/1‑100 (MQ=255)
                                        aTCCGGGCCGACGCTGCCTAAGCGCAGCCAACCAATTTTCTCATAAACTGTCCATTGCTTAACGCGTTCAATGAAGACCGCAACCCGCGCCTCCGGCGtt                                                       <  2:2613377/100‑1 (MQ=255)
                                        aTCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGtt                                                       >  2:606578/1‑100 (MQ=255)
                                        aTCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGtt                                                       >  2:1785771/1‑100 (MQ=255)
                                        aTCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGtt                                                       <  1:2202453/100‑1 (MQ=255)
                                         tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTCAGCAGTTCGGTGACGAACGCATCCCGCGCTTAGGGCGTTa                                                      >  2:3084923/1‑100 (MQ=255)
                                         tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATTCCGCGCGTCGGGCGtga                                                      >  2:386136/1‑98 (MQ=255)
                                         tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTa                                                      <  2:779242/100‑1 (MQ=255)
                                         tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTa                                                      >  2:921725/1‑100 (MQ=255)
                                           cGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAat                                                    <  1:1646246/100‑1 (MQ=255)
                                            gggCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTACGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatg                                                   >  1:1864054/1‑100 (MQ=255)
                                            gggCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatg                                                   >  1:1023739/1‑100 (MQ=255)
                                             ggCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatgg                                                  >  1:922425/1‑100 (MQ=255)
                                                    gCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTg                                           <  2:2832398/100‑1 (MQ=255)
                                                       gATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCg                                        <  1:3649455/100‑1 (MQ=255)
                                                       gATTAAGCGCAGCTGACCAATATTCTCACCAACTGACATTTGATTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCg                                        <  1:187379/100‑1 (MQ=255)
                                                        aTTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGa                                       >  2:632017/1‑100 (MQ=255)
                                                        aTTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGa                                       >  1:2483116/1‑100 (MQ=255)
                                                        aTTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGa                                       >  1:2705624/1‑100 (MQ=255)
                                                         ttAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGaa                                      <  1:2126067/100‑1 (MQ=255)
                                                           aaGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATa                                    <  1:1322389/100‑1 (MQ=255)
                                                            aGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATaa                                   <  1:1717088/100‑1 (MQ=255)
                                                             gcgcAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATaaa                                  <  1:3918798/100‑1 (MQ=255)
                                                              cgcAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAat                                 <  1:105774/100‑1 (MQ=255)
                                                              cgcAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAat                                 <  1:1665992/100‑1 (MQ=255)
                                                              cgcAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAat                                 >  2:3469887/1‑100 (MQ=255)
                                                              cgcAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAat                                 <  2:1707385/100‑1 (MQ=255)
                                                               gcAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAatc                                >  2:457824/1‑100 (MQ=255)
                                                                atcttGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAatca                               <  1:2926680/96‑1 (MQ=255)
                                                                  gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAatcatc                             >  1:2564674/1‑100 (MQ=255)
                                                                   cTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATcc                            <  1:1542344/100‑1 (MQ=255)
                                                                    tGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAatcaagat                           >  2:2882814/1‑96 (MQ=255)
                                                                    tGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCg                           <  1:2065673/100‑1 (MQ=255)
                                                                     gACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGc                          >  2:2306501/1‑100 (MQ=255)
                                                                       ccAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCa                        >  2:3968541/1‑100 (MQ=255)
                                                                       ccAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTACCGAATAAATCATCCGCCa                        >  1:2038254/1‑100 (MQ=255)
                                                                        cAATTTTCTCGTGAGCTGGCAGGTGCTTAAGCAGGTCGGTGGCGAACGCAGCCCGCGCTTCGGGCGTTAATGGGTGGTTGCCGAATAAATCATCCGCCAg                       <  2:1545787/100‑1 (MQ=255)
                                                                         aaTTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGt                      >  2:1214926/1‑100 (MQ=255)
                                                                              tctcATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCTCTTCGGGCGTTAATGGATTGTTGCCGAATAAATCATCCGCCAGTgctgg                 >  1:2041492/1‑100 (MQ=255)
                                                                               ctcATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggc                <  1:1719850/100‑1 (MQ=255)
                                                                               ctcATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggc                <  2:252274/100‑1 (MQ=255)
                                                                                 catcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggctg              >  1:844692/1‑100 (MQ=255)
                                                                                 catcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAACGGATGGTTGCCGAATAAATCATCCGCCAGTgctggctg              >  2:3502067/1‑100 (MQ=255)
                                                                                   tcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCa            >  1:2426833/1‑100 (MQ=255)
                                                                                   tcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCa            >  1:1443194/1‑100 (MQ=255)
                                                                                    caACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAg           >  2:1411467/1‑100 (MQ=255)
                                                                                    caACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAg           >  1:1826267/1‑100 (MQ=255)
                                                                                       cTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGGGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGc        >  1:3878507/1‑100 (MQ=255)
                                                                                        tGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGcc       >  1:2941428/1‑100 (MQ=255)
                                                                                        tGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGcc       >  1:1178215/1‑100 (MQ=255)
                                                                                        tGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGcc       <  2:856078/100‑1 (MQ=255)
                                                                                          tCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAg     <  1:1546354/100‑1 (MQ=255)
                                                                                           cATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGa    >  1:2457331/1‑100 (MQ=255)
                                                                                            aTTTGCTTAAGCAGTCCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGAc   >  2:3037595/1‑100 (MQ=255)
                                                                                             tttgCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACt  >  1:2061613/1‑100 (MQ=255)
                                                                                                |                                                                                                
GACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACT  >  NC_000913/2221764‑2221956

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: