Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,194 A→T 16.8% intergenic (+28/+25) yicH → / ← yicI putative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,1940AT16.8% 526.6 / 95.6 224intergenic (+28/+25)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (84/100);  new base T (23/14);  total (108/116)
Fisher's exact test for biased strand distribution p-value = 7.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

Only 100 of 224 total aligned reads displayed.

AAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCG  >  NC_000913/3832096‑3832291
                                                                                                  |                                                                                                 
aaGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                  <  2:2572389/100‑1 (MQ=255)
 aGTGACTCAAACCATGAATGCTGGCATCGTTTCCCGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                 <  1:3688935/100‑1 (MQ=255)
 aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                 <  1:3920648/100‑1 (MQ=255)
 aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                 <  1:2721566/100‑1 (MQ=255)
    gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGtgc                                                                                              >  1:588003/1‑98 (MQ=255)
      cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  2:1980567/1‑100 (MQ=255)
      cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  1:682157/1‑100 (MQ=255)
      cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  1:3588229/1‑100 (MQ=255)
      cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  2:4021648/1‑100 (MQ=255)
      cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  1:1091233/1‑100 (MQ=255)
       tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                           <  1:3620508/100‑1 (MQ=255)
       tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                           <  1:2409767/100‑1 (MQ=255)
       tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCg                                                                                           >  1:2288554/1‑100 (MQ=255)
        naaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                          >  1:3288774/2‑100 (MQ=255)
        cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                          <  2:1192454/100‑1 (MQ=39)
        cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                          <  1:1304156/100‑1 (MQ=39)
        cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                          >  1:1397435/1‑100 (MQ=255)
        cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACACgg                                                                                          >  2:733454/1‑100 (MQ=255)
         aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGc                                                                                         <  2:172002/100‑1 (MQ=255)
         aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGc                                                                                         >  2:1446544/1‑100 (MQ=255)
          aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                         <  2:413362/99‑1 (MQ=39)
          aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                        <  1:3047801/100‑1 (MQ=255)
             cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                     <  2:1318305/100‑1 (MQ=255)
                gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                  <  2:279885/100‑1 (MQ=255)
                gAATGCCGGCGTCGTTTCCAGCGGTGGAGTTACAACGACGGAGCCAGTGCGGTAAACATCAATCTCGGCGTGTACGACTCAGACGTACCCGGCATTCCAt                                                                                  <  2:2177784/100‑1 (MQ=255)
                 aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                 <  1:1033591/100‑1 (MQ=255)
                  aTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                >  2:760333/1‑100 (MQ=255)
                   tGCTGGCATCGTTTCCAGCGGTGAATTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCACCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                 <  2:1076253/98‑1 (MQ=255)
                   tGCTGGCATCGTTTCCAGCGGTGAAGTTTCATCGCCGGAGCCGGTGCGGTAAGCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                  <  2:2831812/97‑1 (MQ=255)
                   tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                               >  1:137988/1‑100 (MQ=255)
                   tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                               <  2:2303445/100‑1 (MQ=255)
                   tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                               <  2:1505030/100‑1 (MQ=255)
                    gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                              >  1:42200/1‑100 (MQ=255)
                             gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                     >  2:3432388/1‑100 (MQ=255)
                             gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                     <  1:85638/100‑1 (MQ=255)
                             gTTTCCAGCGGTGAAGTTACATCGACGGAGCCCGGGCGGTAAACGTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                     <  2:3894458/100‑1 (MQ=255)
                              tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                    <  1:3656682/100‑1 (MQ=255)
                               ttCCGGCGGTGCTATTACCTCGACGGGGCCGCGACTGTCAAGATCAATGTCGCCGGGTACGACGCAGCCGTCCCCGGCATGCCATCAATAGATTACAACg                                                                   <  2:4090751/100‑1 (MQ=255)
                               ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                                   <  2:4034426/100‑1 (MQ=255)
                               ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCGCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                                   <  2:2000120/100‑1 (MQ=255)
                                tCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGt                                                                  <  1:725148/100‑1 (MQ=255)
                                  cAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                                >  1:3491457/1‑100 (MQ=255)
                                     cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                             <  1:3907630/100‑1 (MQ=255)
                                       gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                           <  1:152394/100‑1 (MQ=255)
                                         gAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAg                                                         <  1:2097605/100‑1 (MQ=255)
                                          aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                        <  2:137988/100‑1 (MQ=255)
                                             ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCa                                                     <  2:1945765/100‑1 (MQ=255)
                                              tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                    <  2:2928480/100‑1 (MQ=255)
                                                cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                  >  2:3687404/1‑100 (MQ=255)
                                                  tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                                <  2:1171666/100‑1 (MQ=255)
                                                    gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                              >  2:3846938/1‑100 (MQ=255)
                                                    gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                              >  1:1285076/1‑100 (MQ=255)
                                                    gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                              >  2:3586067/1‑100 (MQ=255)
                                                      cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                            <  2:681880/100‑1 (MQ=255)
                                                      cGGAGCCAGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATATATTACAACGTAATTGTCAGCGCATTCCCTTGa                                            <  2:10229/100‑1 (MQ=255)
                                                       ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                           >  1:1769288/1‑100 (MQ=255)
                                                       ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                           <  2:435690/100‑1 (MQ=255)
                                                        gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                          <  1:1182487/100‑1 (MQ=255)
                                                        gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                           <  1:3244854/99‑1 (MQ=255)
                                                           ccGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                       >  1:1974271/1‑100 (MQ=255)
                                                           ccGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                       <  2:3820266/100‑1 (MQ=255)
                                                             ggTGCGGTAAACATCAATCTCGCCGGGTACGACTTAGACGTACACGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCTTTGAGGCTTca                                     <  2:1836270/100‑1 (MQ=255)
                                                             ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                     >  1:1580564/1‑100 (MQ=255)
                                                             ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                     <  2:3000055/100‑1 (MQ=255)
                                                              gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                    >  1:3968302/1‑100 (MQ=255)
                                                              gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                     <  2:42200/99‑1 (MQ=255)
                                                                gCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                                  >  2:1866688/1‑100 (MQ=255)
                                                                  ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                                <  2:1087629/100‑1 (MQ=255)
                                                                   gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                               >  1:2879016/1‑100 (MQ=255)
                                                                     aaaCATCAATCTCGCCGGGTACGACTCAGACGTACCCCGCATTCCATCAATAGCTCACAACGTAATTCTCAGCGCCTTCCCTTGCGGCTTCACCACCAGc                             >  2:1874818/1‑100 (MQ=255)
                                                                      aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                            >  2:3980656/1‑100 (MQ=255)
                                                                        cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                          <  1:2954512/100‑1 (MQ=255)
                                                                        cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                          >  1:2985262/1‑100 (MQ=255)
                                                                          tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTg                        <  1:453930/100‑1 (MQ=255)
                                                                          tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACAACCAGCCCCTg                        >  2:829049/1‑100 (MQ=255)
                                                                           cAATATCGCCGGGTACGACTCAGACGTACCCGGCACTCCATCAATAGATTACAACGTAATTGTCCGCGCATTACCCTCAGGCTTCACAACCAGCCCCTGc                       <  2:2920966/100‑1 (MQ=255)
                                                                             ggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc                     >  1:1589104/6‑100 (MQ=37)
                                                                                atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                  >  2:3779626/3‑100 (MQ=38)
                                                                                 tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                 >  1:3055658/2‑100 (MQ=39)
                                                                                 tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                 >  2:3100155/2‑100 (MQ=38)
                                                                                 tgCCAGGTACGTCTGCGTCGTACCCGGCATTCCATCAATAGATTGCAACGTAATTGCCAGCGCATTACCTTGAGGCTTCACTCCCAGCACCCGCTCCCtt                 >  2:4158506/2‑100 (MQ=255)
                                                                                 cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                 <  1:3860173/100‑1 (MQ=255)
                                                                                  gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTAGAACGTAATTGTCAGCGCATTCCATTGAGGCTTCACCACGAGCCTCTGCTCACttt                >  2:2325457/1‑100 (MQ=38)
                                                                                  gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGGGGCTTCACCACCAGCCCCTGCTCACttt                >  1:246157/1‑100 (MQ=38)
                                                                                  gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                >  1:2417985/1‑100 (MQ=39)
                                                                                  gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGGAATGGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTTCTCACttt                >  1:3342718/1‑100 (MQ=38)
                                                                                   ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc               >  1:2684546/1‑100 (MQ=39)
                                                                                   ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc               >  1:3778157/1‑100 (MQ=39)
                                                                                   ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc               <  1:1679548/100‑1 (MQ=255)
                                                                                    cGGGTACGTCTGAGTCGTACCCGGCATTCCATCTATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa              >  2:2765158/1‑100 (MQ=255)
                                                                                         aCGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTTATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg         <  2:2681966/100‑1 (MQ=255)
                                                                                         aCGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg         <  1:1508577/100‑1 (MQ=255)
                                                                                          cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc        <  1:2538068/100‑1 (MQ=255)
                                                                                          cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc        <  2:541678/100‑1 (MQ=255)
                                                                                          cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc        >  1:2873923/1‑100 (MQ=255)
                                                                                            aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGa      <  1:2050851/100‑1 (MQ=255)
                                                                                             cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa     >  2:4110621/1‑100 (MQ=255)
                                                                                             cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa     >  1:2777255/1‑100 (MQ=255)
                                                                                              tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc    <  2:1240584/100‑1 (MQ=255)
                                                                                               cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc   >  2:3182525/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
AAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCG  >  NC_000913/3832096‑3832291

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: