Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,191 C→G 16.6% intergenic (+25/+28) 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,1910CG16.6% 579.3 / 105.1 225intergenic (+25/+28)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (84/102);  new base G (23/14);  total (107/118)
Fisher's exact test for biased strand distribution p-value = 7.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

Only 100 of 231 total aligned reads displayed.

CAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACC  >  NC_000913/3832092‑3832290
                                                                                                   |                                                                                                   
cagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTc                                                                                                     <  1:2935902/100‑1 (MQ=255)
cagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTc                                                                                                     <  1:832427/100‑1 (MQ=255)
    aaGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGAc                                                                                                 <  2:2572389/100‑1 (MQ=255)
     aGTGACTCAAACCATGAATGCTGGCATCGTTTCCCGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                <  1:3688935/100‑1 (MQ=255)
     aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                <  1:3920648/100‑1 (MQ=255)
        gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGtgc                                                                                             >  1:588003/1‑98 (MQ=255)
          cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                           >  1:1091233/1‑100 (MQ=255)
          cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                           >  1:682157/1‑100 (MQ=255)
          cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                           >  1:2134286/1‑100 (MQ=255)
          cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                           >  2:4021648/1‑100 (MQ=255)
           tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                          <  1:2409767/100‑1 (MQ=255)
           tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCg                                                                                          >  1:2288554/1‑100 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  2:1174373/100‑1 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  2:3582060/100‑1 (MQ=39)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         >  1:1397435/1‑100 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACAGCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         <  2:751263/100‑1 (MQ=255)
             aaaCCATGAATGCTGGCATCGTTTCCAGCGGTTAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                        <  2:1384825/100‑1 (MQ=38)
             aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                        <  1:2975549/100‑1 (MQ=39)
              aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcg                                                                                       <  1:154126/100‑2 (MQ=38)
              aaCCATGAATGCTGGCATCGTTTCCAGCCGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcg                                                                                       <  2:824720/100‑2 (MQ=38)
               aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                      <  1:2618610/100‑1 (MQ=255)
                ccgtGAATGCTGGCATCGTTTCCAGCGCTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGCGTACGACTCAGACGTACCCgg                                                                                         <  2:1920852/93‑1 (MQ=255)
                    gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                 <  2:279885/100‑1 (MQ=255)
                     aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  1:151649/100‑1 (MQ=255)
                     aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  1:1033591/100‑1 (MQ=255)
                       tGCTGGCATCGTTTCCAGCGGTGAATTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCACCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  2:1076253/98‑1 (MQ=255)
                       tGCTGGCATCGTTTCCAGCGGTGAAGTTTCATCGCCGGAGCCGGTGCGGTAAGCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                 <  2:2831812/97‑1 (MQ=255)
                       tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                              <  2:2303445/100‑1 (MQ=255)
                       tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                              >  1:137988/1‑100 (MQ=255)
                       tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCAGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                 <  2:2076553/97‑1 (MQ=255)
                        gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                             >  1:42200/1‑100 (MQ=255)
                            gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                         >  2:1349061/1‑100 (MQ=255)
                                cGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACa                                                                     >  1:2049428/1‑100 (MQ=255)
                                 gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                    >  2:3432388/1‑100 (MQ=255)
                                 gTTTCCAGCGGTGAAGTTACATCGACGGAGCCCGGGCGGTAAACGTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                    <  2:3894458/100‑1 (MQ=255)
                                  tttCCAGCGGTGAAGTTCCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                    <  2:1868636/99‑1 (MQ=255)
                                  tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                   <  1:3656682/100‑1 (MQ=255)
                                   ttCCGGCGGTGCTATTACCTCGACGGGGCCGCGACTGTCAAGATCAATGTCGCCGGGTACGACGCAGCCGTCCCCGGCATGCCATCAATAGATTACAACg                                                                  <  2:4090751/100‑1 (MQ=255)
                                   ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                                  <  2:4034426/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)
                                         cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                            <  1:1926515/100‑1 (MQ=255)
                                           gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                          <  2:559754/100‑1 (MQ=255)
                                                gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcgc                                                     >  2:1851460/1‑100 (MQ=255)
                                                 ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCa                                                    <  2:1946577/100‑1 (MQ=255)
                                                  tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                   <  2:2928480/100‑1 (MQ=255)
                                                   aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                                  >  1:3689557/1‑100 (MQ=255)
                                                    cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                                  >  2:1490036/1‑99 (MQ=255)
                                                    cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                 >  2:3687404/1‑100 (MQ=255)
                                                    cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                 >  2:1198538/1‑100 (MQ=255)
                                                      tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                               >  2:1990100/1‑100 (MQ=255)
                                                       cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                              >  1:3034960/1‑100 (MQ=255)
                                                       cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                              >  1:3158364/1‑100 (MQ=255)
                                                        gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                             >  1:1285076/1‑100 (MQ=255)
                                                        gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                             >  2:3846938/1‑100 (MQ=255)
                                                          cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGCCGTACCCGGCATTCCCTCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                            >  1:2697554/1‑99 (MQ=255)
                                                          cGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGa                                           <  2:681880/100‑1 (MQ=255)
                                                          cGGAGCCAGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATATATTACAACGTAATTGTCAGCGCATTCCCTTGa                                           <  2:10229/100‑1 (MQ=255)
                                                            gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                         <  1:2875900/100‑1 (MQ=255)
                                                            gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                         <  1:1182487/100‑1 (MQ=255)
                                                            gAGCCGGTGCGGTAAACATCAAGCTCGCCGGGTACGACTCAGGCGTACCCGGCATTCCATCATTAGCTTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                         <  2:1246797/100‑1 (MQ=255)
                                                             aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTATTTGTCAGCGCATTCCCTTGAgg                                         <  2:2434569/99‑1 (MQ=255)
                                                               ccGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                      <  2:3820266/100‑1 (MQ=255)
                                                                cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTc                                     <  1:2090289/100‑1 (MQ=255)
                                                                 ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                      <  2:4140678/98‑1 (MQ=255)
                                                                 ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                    <  2:3000055/100‑1 (MQ=255)
                                                                 ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                    >  1:1580564/1‑100 (MQ=255)
                                                                 ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                    <  1:3296331/100‑1 (MQ=255)
                                                                  gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                   >  1:3968302/1‑100 (MQ=255)
                                                                  gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                    <  2:42200/99‑1 (MQ=255)
                                                                   tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                                  >  1:1346761/1‑100 (MQ=255)
                                                                    gCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                                 <  1:3609401/100‑1 (MQ=255)
                                                                     cGGTAAACATCAATCTCGCAGGGGACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAATGTCAGCGCATTCCCTTGAGGATTcaccac                                >  2:2191181/1‑100 (MQ=255)
                                                                      ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                               <  2:1087629/100‑1 (MQ=255)
                                                                          aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                           >  2:3980656/1‑100 (MQ=255)
                                                                          aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGATTCACCACCAGc                            >  2:1623044/1‑99 (MQ=255)
                                                                            cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                         >  2:3188435/1‑100 (MQ=255)
                                                                              tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTg                       <  1:453930/100‑1 (MQ=255)
                                                                               cAATATCGCCGGGTACGACTCAGACGTACCCGGCACTCCATCAATAGATTACAACGTAATTGTCCGCGCATTACCCTCAGGCTTCACAACCAGCCCCTGc                      <  2:2920966/100‑1 (MQ=255)
                                                                                    atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                 >  2:3779626/3‑100 (MQ=38)
                                                                                     tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                >  2:3100155/2‑100 (MQ=38)
                                                                                     tgCCAGGTACGTCTGCGTCGTACCCGGCATTCCATCAATAGATTGCAACGTAATTGCCAGCGCATTACCTTGAGGCTTCACTCCCAGCACCCGCTCCCtt                >  2:4158506/2‑100 (MQ=255)
                                                                                     cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                <  1:3860173/100‑1 (MQ=255)
                                                                                     cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                <  2:441839/100‑1 (MQ=255)
                                                                                      gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGGGGCTTCACCACCAGCCCCTGCTCACttt               >  1:246157/1‑100 (MQ=38)
                                                                                      gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt               >  2:1755441/1‑100 (MQ=39)
                                                                                      gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt               >  1:2417985/1‑100 (MQ=39)
                                                                                        cGGGTACGTCTGAGTCGTACCCGGCATTCCATCTATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa             >  2:2765158/1‑100 (MQ=255)
                                                                                         gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg            >  1:3284648/1‑100 (MQ=255)
                                                                                         gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg            <  1:3586084/100‑1 (MQ=255)
                                                                                           gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc          >  1:1750165/1‑100 (MQ=255)
                                                                                             aCGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTTATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg        <  2:2681966/100‑1 (MQ=255)
                                                                                             aCGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg        <  2:829176/100‑1 (MQ=255)
                                                                                              cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc       >  1:1661144/1‑100 (MQ=255)
                                                                                              cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc       <  1:2538068/100‑1 (MQ=255)
                                                                                              cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc       >  1:2873923/1‑100 (MQ=255)
                                                                                                aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg      >  1:3049603/1‑99 (MQ=255)
                                                                                                 cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa    >  2:4110621/1‑100 (MQ=255)
                                                                                                 cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa    >  1:2777255/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACC  >  NC_000913/3832092‑3832290

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: