Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,033,723 G→C 15.6% intergenic (‑239/‑328) yedR ← / → yedS inner membrane protein/pseudogene, outer membrane protein homology; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,033,7230GC15.6% 364.0 / 58.0 137intergenic (‑239/‑328)yedR/yedSinner membrane protein/pseudogene, outer membrane protein homology; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (43/72);  new base C (9/12);  total (53/84)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

Only 100 of 143 total aligned reads displayed.

ATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCA  >  NC_000913/2033625‑2033821
                                                                                                  |                                                                                                  
  tAGCATTTATTACAATTATACCTTTTTTTTATTCTGTGTCCTGATGCTTCAGTTATTAGCCTTTGATCGTCTTGTTTATAGTTGTTGGGCCGGCatgatg                                                                                                 <  2:1649921/100‑1 (MQ=255)
   aGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGc                                                                                                <  1:2578980/100‑1 (MQ=255)
   aGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGc                                                                                                <  1:2136487/100‑1 (MQ=255)
    gCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTTTTAGCCTTTTCTCGTCTTTTTTATATTTTTTGTGCCGGCatgct                                                                                                  >  2:887891/1‑95 (MQ=255)
    gCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGcc                                                                                               <  2:1387765/100‑1 (MQ=255)
    gCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGcc                                                                                               <  1:2408388/100‑1 (MQ=255)
     cACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCg                                                                                              <  1:1392055/100‑1 (MQ=255)
      aCTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCgg                                                                                             <  2:314370/100‑1 (MQ=255)
      aCTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCgg                                                                                             >  1:1739056/1‑100 (MQ=255)
       cTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGc                                                                                            <  2:30891/100‑1 (MQ=255)
       cTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGc                                                                                            >  1:3447159/1‑100 (MQ=255)
        tGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCt                                                                                           >  1:1038359/1‑100 (MQ=255)
        tGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGcc                                                                                           <  1:1429745/100‑2 (MQ=255)
         gATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCtt                                                                                          >  1:2409104/1‑100 (MQ=255)
         gATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGccc                                                                                          <  1:2898316/100‑3 (MQ=255)
          aTTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGccca                                                                                         <  2:3147601/100‑4 (MQ=255)
          aTTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGccca                                                                                         <  2:3383469/100‑4 (MQ=255)
           ttACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGcccaa                                                                                        <  1:1813310/100‑5 (MQ=255)
            tACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCttttt                                                                                       >  2:2892545/1‑100 (MQ=255)
            tACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCttttt                                                                                       >  1:2356249/1‑100 (MQ=255)
            tACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGccca                                                                                         <  1:1275510/98‑4 (MQ=255)
             aCAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCtttttt                                                                                      <  2:422504/100‑1 (MQ=255)
             aCAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCtttttt                                                                                      <  2:2667719/100‑1 (MQ=255)
             aCAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCtttttt                                                                                      <  1:426898/100‑1 (MQ=255)
               atttATTCCTTTTTATTATTCTGTGTTTTGATTCTTTCGTTTTTAGCCTTTTATCTTCTTGTTTTTATTGTTTGGGCCGGCATGTTGCCGGCtttttttt                                                                                    <  2:22430/98‑1 (MQ=255)
               aaTTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCtttttttt                                                                                    <  1:3101176/100‑1 (MQ=255)
                aTTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTa                                                                                   <  1:1165128/100‑1 (MQ=255)
                aTTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATCATGCCGGcccaaaaaat                                                                                   <  1:1419716/100‑10 (MQ=39)
                 ttATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTAt                                                                                  >  1:568870/1‑100 (MQ=255)
                    taCCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGcc                                                                               >  1:2682311/1‑100 (MQ=255)
                    taCCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATAcc                                                                               >  1:3448357/1‑100 (MQ=255)
                     aCCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCt                                                                              >  2:1083683/1‑100 (MQ=255)
                      ccTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATTCCtt                                                                             >  2:1049258/1‑100 (MQ=255)
                      ccTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATTCCtt                                                                             >  2:1685926/1‑100 (MQ=255)
                      ccTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCtt                                                                             >  2:1104869/1‑100 (MQ=255)
                       cTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTc                                                                            <  1:3293654/100‑1 (MQ=255)
                           ttattaTTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTttatta                                                                        >  1:1803230/1‑100 (MQ=255)
                           ttattaTTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCAtt                                                                         <  1:2010365/99‑1 (MQ=255)
                           ttattaTTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTa                                                                        <  2:3448332/100‑1 (MQ=255)
                            tattatTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTaa                                                                       <  2:2768964/100‑1 (MQ=255)
                              ttattCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAAtg                                                                     <  2:354025/100‑1 (MQ=255)
                                attCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAAtgtg                                                                   <  1:3526737/100‑1 (MQ=255)
                                attCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAAtgtg                                                                   >  1:2482543/1‑100 (MQ=255)
                                  tCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGATTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTgcg                                                                 <  2:2550679/100‑1 (MQ=255)
                                    tgtgTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGcc                                                               <  1:2150940/100‑1 (MQ=255)
                                    tgtgTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGcc                                                               <  1:757279/100‑1 (MQ=255)
                                       gtCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGa                                                            >  1:1715311/1‑100 (MQ=255)
                                        tCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGAt                                                           <  2:2293508/100‑1 (MQ=255)
                                              tgCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACAcc                                                     <  2:1038351/100‑1 (MQ=255)
                                                cTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACAcc                                                     >  2:1182639/1‑98 (MQ=255)
                                                 ttCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGc                                                  <  2:3179945/100‑1 (MQ=255)
                                                    cGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGt                                               >  1:3317617/1‑100 (MQ=255)
                                                    cGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGtt                                              >  1:790401/1‑100 (MQ=255)
                                                        attaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTTTTTTTTTATTAATTTGCGCCTGATCACACCAGTCTTTTgg                                           >  2:682744/1‑100 (MQ=255)
                                                        attaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTgg                                           <  1:2736594/100‑1 (MQ=255)
                                                        attaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTgg                                           >  1:3529294/1‑100 (MQ=255)
                                                         ttaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGgc                                          <  1:2716136/100‑1 (MQ=255)
                                                         ttaGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGgc                                          <  2:2831866/100‑1 (MQ=255)
                                                          taGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGgcg                                         <  1:2676976/100‑1 (MQ=255)
                                                            gCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGca                                       <  2:245967/100‑1 (MQ=255)
                                                            gCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGca                                       <  1:1121059/100‑1 (MQ=255)
                                                            gCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGca                                       >  1:3639013/1‑100 (MQ=255)
                                                             ccTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGcaa                                      >  1:995122/1‑100 (MQ=255)
                                                               ttttATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGcaac                                     <  1:2963976/99‑1 (MQ=255)
                                                                tttATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGcaacaa                                   <  2:3578137/100‑1 (MQ=255)
                                                                     cGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCAtt                              <  1:2812293/100‑1 (MQ=255)
                                                                      gTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTg                             <  1:719191/100‑1 (MQ=255)
                                                                         ttGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATa                          <  2:163394/100‑1 (MQ=255)
                                                                             ttATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATAccccc                     <  1:2855281/100‑1 (MQ=255)
                                                                             tcatatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATAcccc                      >  1:2466800/3‑100 (MQ=255)
                                                                              tataTTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATAccccc                     <  1:1017103/100‑1 (MQ=255)
                                                                               atatTTTTTGGGCGGGCATGATGCCGGCTTTTTTTTATGTCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATAccccc                     <  2:102201/99‑1 (MQ=255)
                                                                               atatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTTCGCCTGATCACACCAGCCGTTTCGCGCAACAATCATTTATAccccc                     >  2:1099597/1‑99 (MQ=255)
                                                                               atatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCt                    >  1:1130490/1‑100 (MQ=255)
                                                                                tatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATAccccc                     <  2:1315822/98‑1 (MQ=255)
                                                                                tatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTa                   <  1:680569/100‑1 (MQ=255)
                                                                                tatTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTa                   >  1:3175865/1‑100 (MQ=255)
                                                                                 taaaaaaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTAt                  >  1:1470694/10‑100 (MQ=39)
                                                                                 atTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTAt                  >  2:917594/1‑100 (MQ=255)
                                                                                  ttttttttGCCGGCATGATGCCGGCTTTTTTTTATGCATTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTAtt                 <  2:3453021/100‑2 (MQ=255)
                                                                                  ttttttGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATg                 <  1:2084756/100‑1 (MQ=255)
                                                                                   aaaaaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGt                >  2:3353513/8‑100 (MQ=255)
                                                                                    aaaaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATCTGCGCCTGATCACACCACCCGTTTGGCGCACCAATCATTGCTCACCCCt                    >  2:1467151/7‑95 (MQ=255)
                                                                                     aaaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGttt              >  1:1474964/6‑100 (MQ=255)
                                                                                     aaaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGttt              >  2:1753828/6‑100 (MQ=255)
                                                                                      aaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTc             <  1:1753836/96‑1 (MQ=255)
                                                                                      aaaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTACGCCTGATCACACCAGCCGTTTGGCGCAACAATCAATGATACCCACTATGTTTc             >  2:3151297/5‑100 (MQ=255)
                                                                                       tGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTcc            <  1:840733/100‑1 (MQ=255)
                                                                                       tGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTcc            <  1:2520333/100‑1 (MQ=255)
                                                                                       tGGGACGGCATGATGCCGGCTTTTTTTTATGCCTTCATGAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTc             <  2:586172/99‑1 (MQ=255)
                                                                                       aaagCCGGCATCATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTcc            >  1:2593272/4‑100 (MQ=255)
                                                                                        gggCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCg           <  1:1104878/100‑1 (MQ=255)
                                                                                         ggCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGc          >  1:1943571/1‑100 (MQ=255)
                                                                                           ccGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTc        <  1:2339042/100‑1 (MQ=255)
                                                                                           ccGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTc        >  1:1536934/1‑100 (MQ=255)
                                                                                           ccGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTc        <  1:1849388/100‑1 (MQ=255)
                                                                                           ccGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTc        <  1:3643346/100‑1 (MQ=255)
                                                                                            cGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCa       >  1:577743/1‑100 (MQ=255)
                                                                                              gCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAAc     <  2:3641842/100‑1 (MQ=255)
                                                                                                 tgatgCCGGTTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCa  <  2:1529763/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
ATTAGCACTGATTACAATTATACCTTTTTATTATTCTGTGTCATGATGCTTCCGTTATTAGCCTTTTATCGTCTTGTTTATATTTTTTGGGCCGGCATGATGCCGGCTTTTTTTTATGCCTTCATTAATGTGCGCCTGATCACACCAGCCGTTTGGCGCAACAATCATTGATACCCCCTATGTTTCCGCTCAACTCA  >  NC_000913/2033625‑2033821

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: