Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,800 C→A 15.8% A580S (GCG→TCG)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,8000CA15.8% 642.0 / 102.7 248A580S (GCG→TCG) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (97/111);  new base A (32/7);  total (130/118)
Fisher's exact test for biased strand distribution p-value = 4.13e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-02

Only 100 of 251 total aligned reads displayed.

CGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCAC  >  NC_000913/3832704‑3832899
                                                                                                |                                                                                                   
cGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcg                                                                                                  >  1:4083674/1‑100 (MQ=255)
 gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                 >  2:2350821/1‑100 (MQ=255)
 gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                 >  1:962880/1‑100 (MQ=255)
 gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                                 >  2:488451/1‑100 (MQ=255)
  ccTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgcg                                                                                                >  1:1444425/1‑100 (MQ=255)
   cTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACgcgc                                                                                               <  2:4038264/100‑1 (MQ=255)
   cTAACATGTATTCACAGCCAAGGTAATCCCAAGCCGGATCGTCCGCGAACTCCATCATCATGGCCCGCATCATCGGCGTACACCGCGCCTTCGCAc                                                                                                   >  1:1618102/1‑96 (MQ=255)
      aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGc                                                                                            >  1:3180708/1‑100 (MQ=255)
      aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGc                                                                                            >  2:3621063/1‑100 (MQ=255)
      aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                            <  2:676634/100‑4 (MQ=255)
       cATGTATTGCCGGTCAAGGTACTCACAAGCCGGATCGTCCGGGAACTCCACCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCt                                                                                           <  1:2623701/100‑1 (MQ=255)
       cATGTATTGACGGTCAAGGTAATCACAAGCCGGATTGTCCGGGAACTTCATCATCATGGCCCGCATTATCGGCGTACCCCGAGCGATCGCACGCGCa                                                                                              >  1:2022611/1‑97 (MQ=255)
       cATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgggg                                                                                           <  2:1822960/100‑5 (MQ=255)
         tGTATTGACGGTCAAGGTAATCAGATGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgggg                                                                                           <  1:1699751/98‑5 (MQ=255)
           tATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCAc                                                                                       >  1:2061561/1‑100 (MQ=255)
            aTTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACg                                                                                      >  2:3222542/1‑100 (MQ=255)
               gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                   <  2:316388/100‑1 (MQ=255)
               gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                   >  2:1591013/1‑100 (MQ=255)
                aCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATAc                                                                                  <  1:1311601/100‑1 (MQ=255)
                aCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATAc                                                                                  >  1:2986905/1‑100 (MQ=255)
                 cGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATAc                                                                                  >  1:646194/1‑99 (MQ=255)
                  ggTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAg                                                                                >  1:1841904/1‑100 (MQ=255)
                  ggTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAg                                                                                >  2:2823648/1‑100 (MQ=255)
                   gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGa                                                                               >  2:953931/1‑100 (MQ=255)
                   gTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGa                                                                               <  1:2585930/100‑1 (MQ=255)
                    tCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGAt                                                                              >  2:3914663/1‑100 (MQ=255)
                       aGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                           <  1:2612938/100‑1 (MQ=255)
                       aGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                           <  1:3577466/100‑1 (MQ=255)
                        ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACgg                                                                          >  2:2665144/1‑100 (MQ=255)
                        ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                           >  2:3832527/1‑99 (MQ=255)
                          tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACaagatcggaa                                                                        >  2:3300299/1‑91 (MQ=255)
                          tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                        >  1:2944892/1‑100 (MQ=255)
                          tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                        <  2:3875282/100‑1 (MQ=255)
                          tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                        >  1:3964435/1‑100 (MQ=255)
                          tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                        >  1:369959/1‑100 (MQ=255)
                          tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCAGCGCGTCCGCACGCGCAGCTTCACGATACAGATACGGca                                                                        >  1:3094257/1‑100 (MQ=255)
                             tCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatca                                                                     >  2:1127164/1‑100 (MQ=255)
                                   gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGc                                                               >  1:3829450/1‑100 (MQ=255)
                                   gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGc                                                               <  2:1280292/100‑1 (MQ=255)
                                    ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCa                                                              <  1:774484/100‑1 (MQ=255)
                                    ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCa                                                              >  2:3830025/1‑100 (MQ=255)
                                       gATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCAttt                                                           <  1:4279951/100‑1 (MQ=255)
                                        aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                          <  1:1317413/100‑1 (MQ=255)
                                        aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                          <  2:15192/100‑1 (MQ=255)
                                        aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                          <  2:2365603/100‑1 (MQ=255)
                                           gcccGGGAACTCCTTCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCATCTTCACGATACAGATACGGCATCATGCGGCATTTCAGt                                                       >  2:2991767/3‑100 (MQ=255)
                                           gTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCCTGCGGCATTTCAg                                                        >  1:4105049/1‑99 (MQ=255)
                                            tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                      >  1:4003480/1‑100 (MQ=255)
                                            tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                      <  1:1953011/100‑1 (MQ=255)
                                             ccGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTg                                                     <  1:4295556/100‑1 (MQ=255)
                                               gggAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCg                                                   >  1:3725198/1‑100 (MQ=255)
                                                 cttccgatctatcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTg                                                 <  2:1408225/90‑1 (MQ=255)
                                                    cTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaag                                              >  1:3922014/1‑100 (MQ=255)
                                                     tCCATCATCCTCGCCCGCATCACCGGCGTACCCCCCGCGTTCGCACGCGCGGCTTCACGATACAGATACGGCCTCATGCGGCATTTCAGTTGCGt                                                  <  1:4403410/95‑1 (MQ=255)
                                                     tCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCATGATACAGATACGGCATCATGCGGCATTTGAGTTGCGTgaaga                                             <  2:2973216/100‑1 (MQ=255)
                                                      ccATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                            >  2:3386373/1‑100 (MQ=255)
                                                      ccATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaagaa                                            >  2:2215309/1‑100 (MQ=255)
                                                       catcatcatGGCCCGCTTCCTCTGCGTACCCCGCGCGTTCGCACGCGTAGCTTCGCGATACAGATACGGCACCATGCGGCTTTTCAGTTGCGTgaag                                              <  1:504353/97‑1 (MQ=255)
                                                           atcatcGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                       <  1:1522173/100‑1 (MQ=255)
                                                           atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgagatcggaag                                       >  1:1408220/1‑91 (MQ=255)
                                                           atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGCAGAAGCGca                                       >  2:3287718/1‑100 (MQ=255)
                                                           atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                       >  1:1065730/1‑100 (MQ=255)
                                                                ggCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACcaca                                  <  2:21537/100‑1 (MQ=255)
                                                                 gCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACAt                                 >  1:758125/1‑100 (MQ=255)
                                                                    cGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                              <  1:336491/100‑1 (MQ=255)
                                                                    cGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcac                              <  1:2324853/100‑1 (MQ=255)
                                                                     gCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATcaca                             >  2:2592737/1‑100 (MQ=255)
                                                                       atcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAgg                           >  2:1261586/1‑100 (MQ=255)
                                                                         catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAc                         <  2:13776/100‑1 (MQ=255)
                                                                          atcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACt                        <  1:1910641/100‑1 (MQ=255)
                                                                          atcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACt                        >  2:3820713/1‑100 (MQ=255)
                                                                            cGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtca                      <  2:2181137/100‑1 (MQ=255)
                                                                            cGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtca                      <  2:653886/100‑1 (MQ=255)
                                                                            cGGCGTAACCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAc                         <  1:685672/97‑1 (MQ=255)
                                                                             ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcat                     >  2:4098793/1‑100 (MQ=255)
                                                                               gtgaagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatca                   >  2:2495945/9‑100 (MQ=255)
                                                                                tgaagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat                  >  2:3472559/8‑100 (MQ=255)
                                                                                tgaagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCAGCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat                  >  1:2442936/8‑100 (MQ=255)
                                                                                gTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat                  >  2:1677940/1‑100 (MQ=255)
                                                                                  aagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCg                >  1:1033047/6‑100 (MQ=255)
                                                                                   ccccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt               <  2:4121775/100‑1 (MQ=255)
                                                                                   ccccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt               <  1:3561076/100‑1 (MQ=255)
                                                                                   agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGCAGAAGCGCACCACATCACAGGAATCATCATCGt               >  2:4008638/5‑100 (MQ=255)
                                                                                   agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt               >  2:2037516/5‑100 (MQ=255)
                                                                                   agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt               >  2:1980311/5‑100 (MQ=255)
                                                                                   agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt               >  2:1245209/5‑100 (MQ=255)
                                                                                   agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt               >  1:3252458/5‑100 (MQ=255)
                                                                                    gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa              <  2:4358872/97‑1 (MQ=255)
                                                                                    gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa              >  2:2192455/4‑100 (MQ=255)
                                                                                    gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa              >  2:3535174/4‑100 (MQ=255)
                                                                                    cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa              >  2:115313/1‑100 (MQ=255)
                                                                                    cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa              >  2:1502578/1‑100 (MQ=255)
                                                                                       gcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGc           <  2:2066173/100‑1 (MQ=255)
                                                                                       gcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGc           >  1:2503454/1‑100 (MQ=255)
                                                                                       gcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGc           >  1:2234069/1‑100 (MQ=255)
                                                                                              cGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGc    <  2:1074826/100‑1 (MQ=255)
                                                                                              cGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGc    <  2:369960/100‑1 (MQ=255)
                                                                                              cGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGc    <  2:1452235/100‑1 (MQ=255)
                                                                                               gCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGca   >  1:841375/1‑100 (MQ=255)
                                                                                                cACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGca   >  1:2423104/1‑99 (MQ=255)
                                                                                                |                                                                                                   
CGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGGCAC  >  NC_000913/3832704‑3832899

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: