Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,797 T→G 16.1% N581H (AAC→CAC)  yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,7970TG16.1% 578.5 / 108.6 245N581H (AAC→CAC) yicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (93/111);  new base G (32/7);  total (127/118)
Fisher's exact test for biased strand distribution p-value = 3.66e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

Only 100 of 249 total aligned reads displayed.

CGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGG  >  NC_000913/3832698‑3832896
                                                                                                   |                                                                                                   
 gTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTc                                                                                                    <  1:2005577/100‑1 (MQ=255)
  ttGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCg                                                                                                   <  2:1090625/100‑1 (MQ=255)
  ttGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCg                                                                                                   <  1:4139603/100‑1 (MQ=255)
      cGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcg                                                                                               >  1:4083674/1‑100 (MQ=255)
      cGCCTAACATGTATTCACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGACCTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAc                                                                                                <  2:3885427/99‑1 (MQ=255)
       gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgc                                                                                              >  2:2350821/1‑100 (MQ=255)
       gCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGGTCTGAcgc                                                                                              >  1:3734642/1‑100 (MQ=255)
        ccTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCAcgcg                                                                                             >  1:1444425/1‑100 (MQ=255)
         cTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACgcgc                                                                                            <  1:2152317/100‑1 (MQ=255)
         cTAACATGTATTCACAGCCAAGGTAATCCCAAGCCGGATCGTCCGCGAACTCCATCATCATGGCCCGCATCATCGGCGTACACCGCGCCTTCGCAc                                                                                                >  1:1618102/1‑96 (MQ=255)
            aCATGTATTGGCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcggg                                                                                         <  1:4068128/100‑4 (MQ=255)
            aCATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGc                                                                                         >  2:3621063/1‑100 (MQ=255)
             cATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgggg                                                                                        <  2:1822960/100‑5 (MQ=255)
             cATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACGCGATCAGCATGGCCCGCATCATCGGCGTACCCCGCGCGTGCGAACgcgcgg                                                                                          <  1:968505/98‑3 (MQ=255)
                 tATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCAc                                                                                    <  1:4221403/100‑1 (MQ=255)
                  aTTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACg                                                                                   >  2:3222542/1‑100 (MQ=255)
                    tGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGAt                                                                                 <  2:343520/100‑1 (MQ=255)
                     gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                >  2:1591013/1‑100 (MQ=255)
                     gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                <  2:146694/100‑1 (MQ=255)
                     gACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATa                                                                                >  1:340249/1‑100 (MQ=255)
                      aCGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATAc                                                                               <  2:817712/100‑1 (MQ=255)
                       cGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACa                                                                              <  2:2663346/100‑1 (MQ=255)
                       cGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACa                                                                              >  2:2181315/1‑100 (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)
                             aGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACg                                                                        <  1:2612938/100‑1 (MQ=255)
                              ggTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACgg                                                                       >  2:382266/1‑100 (MQ=255)
                                tAATCACGAGCAGGATCGTCCGGGAACTCCATCATCATGACCCGCATCATCGGCGTACCCCGCGCATCCTCACGCGCAGCTTCACGATACAGATACGGc                                                                      >  1:2310351/1‑99 (MQ=255)
                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACaagatcggaa                                                                     >  2:3300299/1‑91 (MQ=255)
                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                     >  1:2944892/1‑100 (MQ=255)
                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                     <  2:3872409/100‑1 (MQ=255)
                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                     >  1:369959/1‑100 (MQ=255)
                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGca                                                                     >  1:3964435/1‑100 (MQ=255)
                                tAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCAGCGCGTCCGCACGCGCAGCTTCACGATACAGATACGGca                                                                     >  1:3094257/1‑100 (MQ=255)
                                   tCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatca                                                                  >  1:820017/1‑100 (MQ=255)
                                    cacaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatcat                                                                 <  2:4069365/100‑1 (MQ=255)
                                    cacaAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGcatcat                                                                 >  1:1894866/1‑100 (MQ=255)
                                        aGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCgg                                                             <  2:2083699/100‑1 (MQ=255)
                                        aGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCgg                                                             >  2:3997487/1‑100 (MQ=255)
                                         gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGc                                                            >  1:3829450/1‑100 (MQ=255)
                                         gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGc                                                            >  1:3589977/1‑100 (MQ=255)
                                         gCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCCTAACGATACAGATACGGCATCATGCgg                                                             >  1:1449886/1‑99 (MQ=255)
                                          ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCa                                                           <  1:774484/100‑1 (MQ=255)
                                          ccGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCa                                                           >  2:3830025/1‑100 (MQ=255)
                                             gATCGTCCGGGAACTCCATCATCATGGCCCGCATCATGGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACACATACGGCATCATGCGGCAttt                                                        >  2:561091/1‑100 (MQ=255)
                                              aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                       <  2:2594985/100‑1 (MQ=255)
                                              aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                       <  1:1317413/100‑1 (MQ=255)
                                              aTCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTc                                                       <  2:1400681/100‑1 (MQ=255)
                                                 gcccGGGAACTCCTTCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCATCTTCACGATACAGATACGGCATCATGCGGCATTTCAGt                                                    >  2:2991767/3‑100 (MQ=255)
                                                 gTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGt                                                    >  2:1128796/1‑100 (MQ=255)
                                                  tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                   >  2:1715981/1‑100 (MQ=255)
                                                  tCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGtt                                                   <  1:1953011/100‑1 (MQ=255)
                                                   ccGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTg                                                  <  1:4295556/100‑1 (MQ=255)
                                                    cGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGc                                                 >  1:2503296/1‑100 (MQ=255)
                                                     gtgAACTCCATCATCATGGCCCTCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCg                                                <  1:3220134/98‑1 (MQ=255)
                                                       gAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTg                                              <  2:4015879/100‑1 (MQ=255)
                                                          cTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaag                                           >  1:3929699/1‑100 (MQ=255)
                                                          cTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTgaa                                            <  1:3016421/99‑1 (MQ=255)
                                                             catcatcatGGCCCGCTTCCTCTGCGTACCCCGCGCGTTCGCACGCGTAGCTTCGCGATACAGATACGGCACCATGCGGCTTTTCAGTTGCGTgaag                                           <  1:504353/97‑1 (MQ=255)
                                                                 atcatGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGca                                    >  1:1065730/1‑100 (MQ=255)
                                                                     tGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGcaccac                                <  2:1203621/100‑1 (MQ=255)
                                                                      ggCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACcaca                               <  2:21537/100‑1 (MQ=255)
                                                                        cccGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATc                             >  1:279242/1‑100 (MQ=255)
                                                                             atcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAgg                        <  2:2933239/100‑1 (MQ=255)
                                                                             atcatcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAgg                        >  2:1261586/1‑100 (MQ=255)
                                                                               natcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAc                      >  2:3518738/2‑100 (MQ=255)
                                                                               catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAc                      <  2:3256833/100‑1 (MQ=255)
                                                                               catcGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAc                      <  2:13776/100‑1 (MQ=255)
                                                                                  cGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtca                   <  2:653886/100‑1 (MQ=255)
                                                                                  cGGCGTAACCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGAc                      <  1:685672/97‑1 (MQ=255)
                                                                                   ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcat                  <  1:2332190/100‑1 (MQ=255)
                                                                                   ggCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcat                  >  2:4098793/1‑100 (MQ=255)
                                                                                    gggTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatc                 <  2:1994148/98‑1 (MQ=255)
                                                                                    cgtgaagctgcgcgTGCGAACGCGCGGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatc                 >  1:634919/10‑100 (MQ=38)
                                                                                     gtgaagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatca                >  2:2495945/9‑100 (MQ=255)
                                                                                     gtgaagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACtcatca                >  2:2326399/9‑100 (MQ=255)
                                                                                      tgaagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCAGCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTcatcat               >  1:2442936/8‑100 (MQ=255)
                                                                                       tACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCatcatc              <  1:1395863/100‑1 (MQ=255)
                                                                                        aagctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCg             >  1:1033047/6‑100 (MQ=255)
                                                                                         ccccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCg             <  1:72651/99‑1 (MQ=255)
                                                                                         ccccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt            <  1:3561076/100‑1 (MQ=255)
                                                                                         agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt            >  1:1684461/5‑100 (MQ=255)
                                                                                         agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt            >  1:1281525/5‑100 (MQ=255)
                                                                                         agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt            >  2:2037516/5‑100 (MQ=255)
                                                                                         agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt            >  1:3252458/5‑100 (MQ=255)
                                                                                         agctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGt            >  1:4358841/5‑100 (MQ=255)
                                                                                          gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATGCGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa           >  2:3098866/4‑100 (MQ=255)
                                                                                          gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa           <  2:4358872/97‑1 (MQ=255)
                                                                                          gctgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa           >  2:2192455/4‑100 (MQ=255)
                                                                                          cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa           >  2:115313/1‑100 (MQ=255)
                                                                                          cccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTa           >  2:1502578/1‑100 (MQ=255)
                                                                                           ccGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAg          >  2:3690554/1‑100 (MQ=255)
                                                                                            tgcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAgg         >  2:3356881/2‑100 (MQ=255)
                                                                                            cgcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAgg         <  1:205449/100‑1 (MQ=255)
                                                                                             gcgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGc        <  2:1276501/100‑1 (MQ=255)
                                                                                             gcgcgTGCGAACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGc        >  2:2406186/1‑100 (MQ=255)
                                                                                              cgcgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGcc       <  2:1525854/100‑1 (MQ=255)
                                                                                                cgTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGccc      >  2:4134371/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
CGTTGTCGCCTAACATGTATTGACGGTCAAGGTAATCACAAGCCGGATCGTCCGGGAACTCCATCATCATGGCCCGCATCATCGGCGTACCCCGCGCGTTCGCACGCGCAGCTTCACGATACAGATACGGCATCATGCGGCATTTCAGTTGCGTGAAGAAGCGCACCACATCACAGGACTCATCATCGTAGGCCCACGG  >  NC_000913/3832698‑3832896

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: