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,8600TG98.9% 661.8 / ‑1.9 188K44Q (AAA→CAA) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (78/108);  total (79/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 189 total aligned reads displayed.

TGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACT  >  NC_000913/2221763‑2221956
                                                                                                 |                                                                                                
tGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCt                                                                                                <  2:796406/100‑1 (MQ=255)
  aCCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTa                                                                                              <  1:2913778/100‑1 (MQ=255)
   ccGTCTGTGATCATCTCGCGGAGCGCCTCTGTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCGCATCAACTGTCATTTGCTTaa                                                                                             <  2:4878009/100‑1 (MQ=255)
    cGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAg                                                                                            >  2:401553/1‑100 (MQ=255)
     gTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGc                                                                                           <  1:1182852/100‑1 (MQ=255)
     gTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGc                                                                                           <  1:3655833/100‑1 (MQ=255)
     gTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGc                                                                                           <  1:2638691/100‑1 (MQ=255)
      tCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCa                                                                                          <  2:328572/100‑1 (MQ=255)
      tCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCa                                                                                          <  1:3471809/100‑1 (MQ=255)
       cTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAg                                                                                         >  1:4342911/1‑100 (MQ=255)
          tGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTc                                                                                      >  1:1552426/1‑100 (MQ=255)
          tGATCATCTCGCGGATCGCCGCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTc                                                                                      >  1:4341403/1‑100 (MQ=255)
           gATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCg                                                                                     >  2:2244298/1‑100 (MQ=255)
            atcatcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCgg                                                                                    <  1:3116228/100‑1 (MQ=255)
              catcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTg                                                                                  <  1:2210885/100‑1 (MQ=255)
               atcTCGCGGTTCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGa                                                                                 <  2:3457747/100‑1 (MQ=255)
               atcTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGa                                                                                 >  1:4571383/1‑100 (MQ=255)
                tctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGAc                                                                                <  1:4090081/100‑1 (MQ=255)
                 ctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCCTCAACTGTCATTTGCTTAAGCAGCTCGGTGACg                                                                               >  2:4250210/1‑100 (MQ=255)
                 ctcGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACg                                                                               >  1:1663774/1‑100 (MQ=255)
                   cgcgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGaa                                                                             <  1:4064627/100‑1 (MQ=255)
                   cgcgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGaa                                                                             <  1:3173449/100‑1 (MQ=255)
                   cgcgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGaa                                                                             >  2:4027766/1‑100 (MQ=255)
                    gcgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAAc                                                                            <  2:468938/100‑1 (MQ=255)
                     cgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACg                                                                           <  1:4859822/100‑1 (MQ=255)
                     cgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACg                                                                           <  2:4222164/100‑1 (MQ=255)
                     cgGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACg                                                                           <  1:1297584/100‑1 (MQ=255)
                      ggATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGc                                                                          >  2:3634298/1‑100 (MQ=255)
                       gATCGCCTCTTTCGTGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCa                                                                         <  1:5042350/100‑1 (MQ=255)
                        aTCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCAt                                                                        <  2:3529828/100‑1 (MQ=255)
                         tCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATc                                                                       >  2:4145871/1‑100 (MQ=255)
                          cGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATcc                                                                      >  2:4224877/1‑100 (MQ=255)
                          cGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATcc                                                                      <  2:4302807/100‑1 (MQ=255)
                          cGCCTCGTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTGTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATcc                                                                      <  2:2887477/100‑1 (MQ=255)
                           gCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATccc                                                                     >  1:3644911/1‑100 (MQ=255)
                           gCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATccc                                                                     >  2:1077988/1‑100 (MQ=255)
                              tctTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCcgcg                                                                  >  2:4855708/1‑100 (MQ=255)
                              tctTTCGGGTTATCCGGGCCGACGCTGATTAAGCGAAGCTGACCAATTTTCTCCTCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCcgcg                                                                  <  1:5298137/100‑1 (MQ=255)
                                  tCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTc                                                              <  1:3399222/100‑1 (MQ=255)
                                   cGGGTTATTCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCg                                                             <  2:267341/100‑1 (MQ=255)
                                   cGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCg                                                             <  1:4623438/100‑1 (MQ=255)
                                    gggTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCgg                                                            >  1:653372/1‑100 (MQ=255)
                                     ggTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCggg                                                           <  2:3342064/100‑1 (MQ=255)
                                       ttATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCg                                                         <  1:991404/100‑1 (MQ=255)
                                        tATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGt                                                        >  2:1810209/1‑100 (MQ=255)
                                        tATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTAATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGt                                                        >  2:669736/1‑100 (MQ=255)
                                         aTCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGtt                                                       <  2:729632/100‑1 (MQ=255)
                                          tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTa                                                      <  2:4588629/100‑1 (MQ=255)
                                          tCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTa                                                      <  1:1222755/100‑1 (MQ=255)
                                            ctggCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAat                                                    <  2:1451543/98‑1 (MQ=255)
                                               gCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatgga                                                 <  1:3053882/100‑1 (MQ=255)
                                               gCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatgga                                                 >  2:3464574/1‑100 (MQ=255)
                                                ccGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatggat                                                <  2:2150259/100‑1 (MQ=255)
                                                  gACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAatggatgg                                              <  1:1840977/100‑1 (MQ=255)
                                                   aCGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGt                                             <  1:3097767/100‑1 (MQ=255)
                                                   aCGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGt                                             <  2:2202122/100‑1 (MQ=255)
                                                    cGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGtt                                            >  1:459595/1‑100 (MQ=255)
                                                     gCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTg                                           <  2:4512732/100‑1 (MQ=255)
                                                     gCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTg                                           <  1:1423053/100‑1 (MQ=255)
                                                       tGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGcc                                         <  1:2422846/100‑1 (MQ=255)
                                                        gATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCg                                        <  2:5442727/100‑1 (MQ=255)
                                                        gATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCg                                        <  2:3935435/100‑1 (MQ=255)
                                                        gATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCg                                        <  2:3706172/100‑1 (MQ=255)
                                                         aTTAAGCGCAGCTGACCAATTTTCTCATCAACTGCCATTTGCTTAAGCAGTTCGCTTACGAACGCATCCCGCGCTTCGGCCGTCAACGTCTGATTGCCGa                                       >  2:4594477/1‑100 (MQ=255)
                                                            aaGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATa                                    <  2:3980831/100‑1 (MQ=255)
                                                             aGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATaa                                   >  1:816310/1‑100 (MQ=255)
                                                               cgcAGCTGACCAATTTTCTCATCAACTGTCATGTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGGGTTAATGGATGGTTGCCGAATAAat                                 <  2:1244316/100‑1 (MQ=255)
                                                                   gCTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAatcatc                             >  1:3715701/1‑100 (MQ=255)
                                                                    cTGACCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATcc                            <  2:4430851/100‑1 (MQ=255)
                                                                       aCCAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGcc                         >  1:4758551/1‑100 (MQ=255)
                                                                        ccAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCa                        >  2:3690693/1‑100 (MQ=255)
                                                                        ccAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCa                        <  1:5236531/100‑1 (MQ=255)
                                                                         cAATTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAg                       <  1:572246/100‑1 (MQ=255)
                                                                          aaTTTTCTCATCACCTGACATTTGCTTAAGCAGTTCGGTGACGAACGCACACCGCGCTTCGTGCGTTAATCTATGCTTGCCGAATAAATCATCCGCCAGt                      >  2:2972041/1‑100 (MQ=255)
                                                                          aaTTTTCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGt                      <  1:2956981/100‑1 (MQ=255)
                                                                              ttCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctg                  <  2:1663069/100‑1 (MQ=255)
                                                                              ttCTCATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctg                  >  1:3127482/1‑100 (MQ=255)
                                                                               tctcATCAACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctgg                 >  2:2731544/1‑100 (MQ=255)
                                                                                ctcATCAACTGTCATTTGCTTGAGCGGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggc                <  2:4468571/100‑1 (MQ=255)
                                                                                  catcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTgctggctg              >  1:1814451/1‑100 (MQ=255)
                                                                                  catcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGGTGCCGAATAAATCATCCGCCAGTgctggctg              >  2:649318/1‑100 (MQ=255)
                                                                                   atcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGc             <  2:497443/100‑1 (MQ=255)
                                                                                    tcgaCTGTGATTTGCTTGAGCAGTGCGGGGACGGACGCAGCCCGCGCCTCGGGCGTTAAAGGATGGTTGCCGAATAAATCATCCCGCAGTGCTGGCTGCa            <  2:4715438/97‑1 (MQ=255)
                                                                                    tcaACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCa            >  2:4544327/1‑100 (MQ=255)
                                                                                     caACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAg           <  2:104900/100‑1 (MQ=255)
                                                                                     caACTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAg           >  2:2446459/1‑100 (MQ=255)
                                                                                     caACTGTCATTAGCTTTAGCAGTACGATGACCACCGCTTCCCGCGCTTCGGGCGTTAATGTATGGTTTCCGAATAAACAATGCTCTAGTTCTGGCTGctg           >  2:298875/1‑98 (MQ=255)
                                                                                      aaCTGTCATTTGCTTAAGCAGTTCGGTGACGAGCGCATCCCGCGCTTCGGGCGTTAATGGATCGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAgg          >  2:4995286/1‑100 (MQ=255)
                                                                                       aCTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAggg         <  2:2559620/100‑1 (MQ=255)
                                                                                       aCTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAggg         <  1:5231313/100‑1 (MQ=255)
                                                                                        cTGTCATTTGCTTAATCATTTCTGTGACGAACGCATCTCGCTCTTCGGTCGTTATTTTATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGc        >  2:1777370/1‑100 (MQ=255)
                                                                                        cTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCGAGTGCTGGCGGCAGGGc        >  1:5123706/1‑100 (MQ=255)
                                                                                        cTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGc        <  2:3837258/100‑1 (MQ=255)
                                                                                        cTGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGACGAATCAATCATCCGCCAGTGCTGGCTGCAGGGc        >  2:4452707/1‑100 (MQ=255)
                                                                                         tGTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGcc       <  2:4039051/100‑1 (MQ=255)
                                                                                          gTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCa      >  2:5137280/1‑100 (MQ=255)
                                                                                          gTCATTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCa      <  1:1966725/100‑1 (MQ=255)
                                                                                             aTTTGCTTAAGCAGTTCGGTGACGAACGTATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGAc   >  1:3222421/1‑100 (MQ=255)
                                                                                             aTTTGCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGAc   <  1:4924834/100‑1 (MQ=255)
                                                                                              tttgCTTAAGCAGCTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACt  >  2:873439/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGACCGTCTTTGATCATCTCGCGGATCGCCTCTTTCGGGTTATCCGGGCCGACGCTGATTAAGCGCAGCTGACCAATTTTCTCATCAACTGTCATTTTCTTAAGCAGTTCGGTGACGAACGCATCCCGCGCTTCGGGCGTTAATGGATGGTTGCCGAATAAATCATCCGCCAGTGCTGGCTGCAGGGCCAGACT  >  NC_000913/2221763‑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: ‑