Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,194 A→T 17.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,1940AT17.8% 366.8 / 74.2 165intergenic (+28/+25)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (69/66);  new base T (16/13);  total (86/79)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 169 total aligned reads displayed.

CAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCG  >  NC_000913/3832095‑3832291
                                                                                                   |                                                                                                 
caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                   >  1:2678060/1‑100 (MQ=255)
caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                   <  2:1628882/100‑1 (MQ=255)
    tGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTa                                                                                               <  1:2534031/100‑1 (MQ=255)
       cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  1:3054367/1‑100 (MQ=255)
       cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  1:728589/1‑100 (MQ=255)
       cTCAAACCATGAATGATGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                            >  1:687890/1‑100 (MQ=255)
        tCAAACCATGAATGCTGGCATCGTTTGCAGCGGTGAAGTTACAGCGACGGAGCCGGTGCAGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                           <  2:1090815/100‑1 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                          <  2:2959563/100‑1 (MQ=39)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                          >  2:1204280/1‑100 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                          >  2:389330/1‑100 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                          >  1:1869367/1‑100 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                          >  1:3629093/1‑100 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCAGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                          <  2:80997/100‑1 (MQ=39)
          aaaCCATTAATGCTGGCATTGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGGGTCGTACCCGGc                                                                                         <  1:2390229/100‑1 (MQ=255)
          aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                         <  1:1101802/100‑1 (MQ=39)
          aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                         <  2:1702148/100‑1 (MQ=39)
          aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGc                                                                                         <  2:2327137/100‑1 (MQ=255)
           aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGc                                                                                         >  1:2516761/1‑99 (MQ=255)
            aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                       >  2:2577969/1‑100 (MQ=255)
             ccATGAATGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTTGCCGGGTACGTCTGAGTCGTACCCGGcgag                                                                                      <  1:1533721/100‑4 (MQ=37)
             ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAtt                                                                                      <  1:1460035/100‑1 (MQ=255)
              catgcatgCTGGCATAGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGCTAAACATCAATCTCGCCGGGTACGACTCAGACGTAGCCGGCATTc                                                                                     <  2:3407092/100‑1 (MQ=255)
              cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                     >  1:481422/1‑100 (MQ=255)
                tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                                   <  1:3219137/100‑1 (MQ=255)
                 gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAt                                                                                  <  2:728586/100‑1 (MQ=255)
                   aTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                <  1:2758374/100‑1 (MQ=255)
                    tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGTAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                               <  2:2457831/100‑1 (MQ=255)
                    tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCaa                                                                               <  2:3553931/100‑1 (MQ=255)
                     gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                              <  1:2921704/100‑1 (MQ=255)
                     gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                              <  1:1849941/100‑1 (MQ=255)
                      cTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGGAAACATCAATCTCGCCGGGTACGACCCAGACGTACCCCGAATTCCATCaa                                                                               >  2:1861809/1‑98 (MQ=255)
                        ggCATCGTTTCCAGCGGTGAAGTTACATCGACGTAGCCCGTGCGGTAAACATCAATCGCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                           <  2:2379567/100‑1 (MQ=255)
                          cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                         >  2:2489809/1‑100 (MQ=255)
                          cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                         <  2:1050592/100‑1 (MQ=255)
                           aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTa                                                                        <  2:2895981/100‑1 (MQ=255)
                            tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTAc                                                                       <  1:1790434/100‑1 (MQ=255)
                               tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                    <  2:261967/100‑1 (MQ=255)
                                ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCACTAGATTACAAc                                                                    >  2:1750725/1‑99 (MQ=255)
                                  ccAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                                 >  1:1290216/1‑100 (MQ=255)
                                   cAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAATCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                                >  2:1936794/1‑100 (MQ=255)
                                    aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTTAt                                                               >  2:3639304/1‑100 (MQ=255)
                                    aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                               >  1:3263771/1‑100 (MQ=255)
                                    aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                               <  1:3142542/100‑1 (MQ=255)
                                    aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCCCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                               <  1:79525/100‑1 (MQ=255)
                                      cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                             <  1:1223384/100‑1 (MQ=255)
                                      cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                             <  1:293214/100‑1 (MQ=255)
                                        gcgAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                           <  1:2357144/98‑1 (MQ=255)
                                        gTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTc                                                           <  2:3570425/100‑1 (MQ=255)
                                            aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                       <  1:3527186/100‑1 (MQ=255)
                                             gTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcgc                                                      >  2:266269/1‑100 (MQ=255)
                                               tACATCGACGGAGCCGGTGCGGTAAACATCAATTTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATGAATAGATTACAACGTAATTGTCAGCGCAt                                                    >  2:2195094/1‑100 (MQ=255)
                                               tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                    >  1:894964/1‑100 (MQ=255)
                                               tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                    <  2:556647/100‑1 (MQ=255)
                                                 cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCCTCAATAGATTACAACGTACTTGTCAGCGCATTc                                                  >  2:714977/1‑100 (MQ=255)
                                                 cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCCTTc                                                  >  2:3021991/1‑100 (MQ=255)
                                                  aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                 <  2:2724929/100‑1 (MQ=255)
                                                    cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                               >  1:2803938/1‑100 (MQ=255)
                                                     gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                              >  1:3254531/1‑100 (MQ=255)
                                                     gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                              >  1:2979169/1‑100 (MQ=255)
                                                        ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCCTGAg                                           <  1:1597124/100‑1 (MQ=255)
                                                         gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                          <  2:1058071/100‑1 (MQ=255)
                                                         gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                          <  1:2095910/100‑1 (MQ=255)
                                                          aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAGTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                         <  2:223459/100‑1 (MQ=255)
                                                          aGCCGGTGCCGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCCATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                         <  2:879210/100‑1 (MQ=255)
                                                                  cGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcaccac                                 >  1:2867549/1‑100 (MQ=255)
                                                                   ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCCTTCCCTTGAGGCTTcaccac                                 >  1:3275100/1‑99 (MQ=255)
                                                                     tAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAg                              <  2:2702664/100‑1 (MQ=255)
                                                                     tAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAg                              <  2:2444343/100‑1 (MQ=255)
                                                                       aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                            >  1:148363/1‑100 (MQ=255)
                                                                       aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                            <  2:586707/100‑1 (MQ=255)
                                                                            cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                       <  1:2233970/100‑1 (MQ=255)
                                                                            cAATCTCGCCGGCTACGACTCAGACGTACCCGGCATTCCATCCATAGATTACAACGTAAATGCCAGCACACTCCCTTGAGGCTTCACCACCAGCCCCtcc                       >  2:1779290/1‑98 (MQ=255)
                                                                             aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                      <  1:1754845/100‑1 (MQ=255)
                                                                              aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc                     <  1:1657993/100‑1 (MQ=255)
                                                                               tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                    >  1:1151322/1‑100 (MQ=255)
                                                                               gaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                    >  1:1249765/5‑100 (MQ=38)
                                                                                 tcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAc                   <  1:2255915/99‑1 (MQ=255)
                                                                                 tcGCCGGGTACAACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAGTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                  >  2:1910907/1‑100 (MQ=255)
                                                                                 atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                  >  2:3400784/3‑100 (MQ=38)
                                                                                 atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                  >  2:1831511/3‑100 (MQ=38)
                                                                                 atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                  >  1:1485867/3‑100 (MQ=38)
                                                                                 atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                  >  1:3424325/3‑100 (MQ=38)
                                                                                  tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                 >  1:1952493/2‑100 (MQ=38)
                                                                                  tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                 >  1:811777/2‑100 (MQ=39)
                                                                                  cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                 >  1:1618380/1‑100 (MQ=255)
                                                                                   gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                <  2:1485859/100‑1 (MQ=39)
                                                                                   accGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                >  2:1399773/2‑100 (MQ=39)
                                                                                    ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTCATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc               >  2:1434802/1‑100 (MQ=39)
                                                                                    ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCCGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                >  2:2550138/1‑99 (MQ=38)
                                                                                     cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa              >  2:2039723/1‑100 (MQ=255)
                                                                                     cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa              >  2:2506026/1‑100 (MQ=255)
                                                                                      gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg             >  1:2309106/1‑100 (MQ=255)
                                                                                        gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc           <  2:1158317/100‑1 (MQ=255)
                                                                                        gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc           <  2:3376262/100‑1 (MQ=255)
                                                                                        gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc           <  2:687887/100‑1 (MQ=255)
                                                                                          aggACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg         >  1:2470165/3‑100 (MQ=255)
                                                                                           cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc        <  1:2625832/100‑1 (MQ=255)
                                                                                            gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg       >  1:3561800/1‑100 (MQ=255)
                                                                                             aCTCAGACGTACCCGGCATTCCATCAATAGATTCCAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCCCACTTTCAGCCTGc        >  2:13250/1‑98 (MQ=255)
                                                                                                cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc   <  1:1916623/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCG  >  NC_000913/3832095‑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: