Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,188 A→T 16.1% intergenic (+22/+31) 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,1880AT16.1% 541.5 / 94.6 229intergenic (+22/+31)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (86/106);  new base T (23/14);  total (109/120)
Fisher's exact test for biased strand distribution p-value = 7.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-02

Only 100 of 232 total aligned reads displayed.

AAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGA  >  NC_000913/3832089‑3832287
                                                                                                   |                                                                                                   
  gcagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACt                                                                                                   >  2:85638/1‑100 (MQ=255)
  gcagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACt                                                                                                   >  1:183282/1‑100 (MQ=255)
   cagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTc                                                                                                  <  2:3387129/100‑1 (MQ=255)
   cagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTc                                                                                                  <  1:832427/100‑1 (MQ=255)
   cagcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTc                                                                                                  <  1:2935902/100‑1 (MQ=255)
     gcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAg                                                                                                <  1:2470116/100‑1 (MQ=255)
        aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                             <  2:3958372/100‑1 (MQ=255)
        aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                             <  1:3920648/100‑1 (MQ=255)
        aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                             <  1:2721566/100‑1 (MQ=255)
         gTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGt                                                                                            >  2:60260/1‑100 (MQ=255)
           gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGtgc                                                                                          >  1:588003/1‑98 (MQ=255)
           gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                          >  1:3956057/1‑100 (MQ=255)
            aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGGGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcc                                                                                         <  2:3814664/100‑1 (MQ=255)
             cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                        >  1:1091233/1‑100 (MQ=255)
             cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                        >  2:1980567/1‑100 (MQ=255)
              tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCg                                                                                       >  1:2288554/1‑100 (MQ=255)
               cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                      <  2:1192454/100‑1 (MQ=39)
               cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                      >  1:1397435/1‑100 (MQ=255)
                aaaCCATGAATGCTGGCATCGTTTCCAGCGGTTAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                     <  2:1384825/100‑1 (MQ=38)
                aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                     <  1:2975549/100‑1 (MQ=39)
                 aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcg                                                                                    <  1:154126/100‑2 (MQ=38)
                 aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                     <  2:413362/99‑1 (MQ=39)
                   ccgtGAATGCTGGCATCGTTTCCAGCGCTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGCGTACGACTCAGACGTACCCgg                                                                                      <  2:1920852/93‑1 (MQ=255)
                        aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                             <  1:1033591/100‑1 (MQ=255)
                        aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                             <  1:151649/100‑1 (MQ=255)
                        aaTGCTGGCATCGTCTCCAGCGGTGTAGATATTTCGATGGAGCCGGTGCGGCAAACATCAATCTCGCCGGGTACGACTCAGACGCTCCAGGTATTCCATc                                                                             <  2:3989460/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:1505030/100‑1 (MQ=255)
                          tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCAGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                              <  2:2076553/97‑1 (MQ=255)
                               gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                      >  2:2922934/1‑100 (MQ=255)
                               gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                      >  2:1349061/1‑100 (MQ=255)
                               gCATCGGTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGGGGTAAACATCAATCTCGCCTTGTACGACTCAGACGTACCCGGCATTCCATCAATa                                                                         <  2:2579037/97‑1 (MQ=255)
                                cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                     <  2:2288549/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)
                                              gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                       <  2:559754/100‑1 (MQ=255)
                                              gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                       <  2:4022513/100‑1 (MQ=255)
                                                gAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAg                                                     <  1:2097605/100‑1 (MQ=255)
                                                   gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcgc                                                  >  1:4144749/1‑100 (MQ=255)
                                                   gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcgc                                                  <  1:2411017/100‑1 (MQ=255)
                                                   gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcgc                                                  <  2:1004339/100‑1 (MQ=255)
                                                      aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                               >  1:3689557/1‑100 (MQ=255)
                                                       cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                              >  2:3687404/1‑100 (MQ=255)
                                                       cATCGACGCAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAGTAAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                              <  2:231111/100‑1 (MQ=255)
                                                         tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                            <  2:1171666/100‑1 (MQ=255)
                                                          cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                           >  1:3034960/1‑100 (MQ=255)
                                                          cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                           >  1:3158364/1‑100 (MQ=255)
                                                           gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                          >  2:3846938/1‑100 (MQ=255)
                                                           gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                          >  2:3586067/1‑100 (MQ=255)
                                                           gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                          >  1:1285076/1‑100 (MQ=255)
                                                           gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                          >  2:688963/1‑100 (MQ=255)
                                                              ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                       <  2:409006/100‑1 (MQ=255)
                                                              ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                       <  2:435690/100‑1 (MQ=255)
                                                               gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                      <  1:1182487/100‑1 (MQ=255)
                                                                aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTATTTGTCAGCGCATTCCCTTGAgg                                      <  2:2434569/99‑1 (MQ=255)
                                                                aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                     <  2:2712416/100‑1 (MQ=255)
                                                                  ccGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                   >  1:1974271/1‑100 (MQ=255)
                                                                    ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                 >  1:1580564/1‑100 (MQ=255)
                                                                    ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                 <  2:3000055/100‑1 (MQ=255)
                                                                     gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                <  1:3573908/100‑1 (MQ=255)
                                                                     gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                 <  2:42200/99‑1 (MQ=255)
                                                                      tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                               >  1:1346761/1‑100 (MQ=255)
                                                                       gCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                              >  2:1866688/1‑100 (MQ=255)
                                                                       gCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                              <  1:3609401/100‑1 (MQ=255)
                                                                         ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                            <  2:1087629/100‑1 (MQ=255)
                                                                          gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                           >  1:2879016/1‑100 (MQ=255)
                                                                          gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                           <  2:1644174/100‑1 (MQ=255)
                                                                           tAAACATCAATCTCGCCTGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAg                          <  2:3549748/100‑1 (MQ=255)
                                                                            aaaCATCAATCTCGCCGGGTACGACTCAGACGTACCCCGCATTCCATCAATAGCTCACAACGTAATTCTCAGCGCCTTCCCTTGCGGCTTCACCACCAGc                         >  2:1874818/1‑100 (MQ=255)
                                                                               cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTACTTGTCAGCGCATTCCCTTGAGGCATCACCACCAGcccc                      >  1:659274/1‑100 (MQ=255)
                                                                               cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                      >  2:3188435/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)
                                                                                   aaTCTCGCCGGGTACTACTCAGACTAACCCGGCATTCCATCAATACATTACAACGTACCTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                   <  2:3838489/99‑1 (MQ=255)
                                                                                        tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt             >  2:3100155/2‑100 (MQ=38)
                                                                                        tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt             >  1:3055658/2‑100 (MQ=39)
                                                                                        tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt             >  2:4105053/2‑100 (MQ=38)
                                                                                        tgCCAGGTACGTCTGCGTCGTACCCGGCATTCCATCAATAGATTGCAACGTAATTGCCAGCGCATTACCTTGAGGCTTCACTCCCAGCACCCGCTCCCtt             >  2:4158506/2‑100 (MQ=255)
                                                                                         gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGGGGCTTCACCACCAGCCCCTGCTCACttt            >  1:246157/1‑100 (MQ=38)
                                                                                         gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt            >  2:1755441/1‑100 (MQ=39)
                                                                                         gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGGAATGGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTTCTCACttt            >  1:3342718/1‑100 (MQ=38)
                                                                                         gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt            <  1:3876759/100‑1 (MQ=255)
                                                                                          ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc           >  1:2684546/1‑100 (MQ=39)
                                                                                           cGGGTACGTCTGAGTCGTACCCGGCATTCCATCTATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa          >  2:2765158/1‑100 (MQ=255)
                                                                                           cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa          >  2:126253/1‑100 (MQ=255)
                                                                                            gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg         >  1:3284648/1‑100 (MQ=255)
                                                                                              gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc       <  1:2304022/100‑1 (MQ=255)
                                                                                              gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc       >  1:1750165/1‑100 (MQ=255)
                                                                                                 cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    <  1:2538068/100‑1 (MQ=255)
                                                                                                 cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    >  1:1661144/1‑100 (MQ=255)
                                                                                                 cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    >  1:3793940/1‑100 (MQ=255)
                                                                                                 cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    <  1:3665237/100‑1 (MQ=255)
                                                                                                 cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc    >  1:2873923/1‑100 (MQ=255)
                                                                                                  gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg   >  1:3644640/1‑100 (MQ=255)
                                                                                                   aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGa  <  1:2050851/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGA  >  NC_000913/3832089‑3832287

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: