Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,410 C→A 16.0% intergenic (‑9/+105) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4100CA16.0% 503.1 / 74.5 181intergenic (‑9/+105)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base C (83/69);  new base A (20/9);  total (103/78)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 188 total aligned reads displayed.

CGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAG  >  NC_000913/4070312‑4070509
                                                                                                  |                                                                                                    
cgcgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTccc                                                                                                     >  2:3900383/1‑100 (MQ=255)
       tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGc                                                                                              >  1:1376128/1‑100 (MQ=255)
       tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGc                                                                                              <  2:4079810/100‑1 (MQ=255)
          tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                          <  2:1774053/100‑1 (MQ=255)
          tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                          <  1:4295089/100‑1 (MQ=255)
          tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGa                                                                                          <  1:2563599/100‑1 (MQ=255)
           cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                         <  1:2058192/100‑1 (MQ=255)
            ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATg                                                                                        >  2:1815326/1‑100 (MQ=255)
            ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtat                                                                                         <  2:2178009/100‑4 (MQ=255)
              ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATCAGCCTCACGGGAGGCTGTATGTc                                                                                      <  1:3367177/100‑1 (MQ=38)
                gcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCaa                                                                                    >  2:469994/1‑100 (MQ=255)
                 cgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTGTATGTcgt                                                                                    <  1:74599/100‑3 (MQ=38)
                  ggCAGGTGAAACGGTCTGTCCACTGCAATGGAGATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAc                                                                                  <  1:3767311/100‑1 (MQ=255)
                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAc                                                                                  <  1:2280162/100‑1 (MQ=255)
                   gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                 <  2:3244464/100‑1 (MQ=255)
                   gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                 >  1:2218207/1‑100 (MQ=255)
                    cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                 >  1:1030444/1‑99 (MQ=255)
                    cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATCAGCCTCACGGGAGGCTGTATGTcgttaat                                                                                <  2:233100/100‑7 (MQ=25)
                     aaGTGAAACGGTCTGTACACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGACAAAccc                                                                                <  1:1746484/99‑1 (MQ=255)
                         gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaaaa                                                                           >  2:2250879/1‑100 (MQ=255)
                         gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaaaa                                                                           >  1:252615/1‑100 (MQ=255)
                         gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaaaa                                                                           >  2:1064829/1‑100 (MQ=255)
                           aaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACt                                                                         >  2:426017/1‑100 (MQ=255)
                              ggTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCa                                                                      >  2:2576080/1‑100 (MQ=255)
                               gTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCAc                                                                     <  1:3363898/100‑1 (MQ=255)
                                 cTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCg                                                                   <  2:3430045/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCgg                                                                  <  2:4051041/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCgg                                                                  <  1:679247/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCgg                                                                  <  1:679087/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCgg                                                                  >  2:553539/1‑100 (MQ=255)
                                   gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGa                                                                 >  1:2802947/1‑100 (MQ=255)
                                    tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAc                                                                <  2:2691682/100‑1 (MQ=255)
                                      cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                                              >  1:829044/1‑100 (MQ=255)
                                         tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           <  2:872099/100‑1 (MQ=255)
                                         tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           <  2:862687/100‑1 (MQ=255)
                                         tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           <  2:670394/100‑1 (MQ=255)
                                          gCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTccc                                                          >  2:2592378/1‑100 (MQ=255)
                                          gCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTccc                                                          >  1:857486/1‑100 (MQ=255)
                                          gCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           >  1:3083273/1‑99 (MQ=255)
                                              tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCg                                                      >  1:1323405/1‑100 (MQ=255)
                                               gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                     >  2:1048499/1‑100 (MQ=255)
                                               gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                     >  2:3931366/1‑100 (MQ=255)
                                               gCATATTGGTGATGTGCATACAAGGCACCTTTATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                     >  2:2967286/1‑100 (MQ=255)
                                                cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                                    <  1:3545186/100‑1 (MQ=255)
                                                cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                                    <  2:3834122/100‑1 (MQ=255)
                                                cATATTGGTGATTTGCATACAAGGCACCGTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                                    <  1:1548375/100‑1 (MQ=255)
                                                   atTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  2:2086122/1‑100 (MQ=255)
                                                    ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  1:286827/1‑99 (MQ=255)
                                                    ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  2:3012459/1‑99 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              >  1:3015432/1‑100 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              >  1:949068/1‑100 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              <  2:1582188/100‑1 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              >  1:142880/1‑100 (MQ=255)
                                                       gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                             <  1:3171715/100‑1 (MQ=255)
                                                       gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                             >  1:609712/1‑100 (MQ=255)
                                                        tGATCTGCATACAAGGCACCTTAATTAACCACATACAGCCGCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTggg                                            <  2:933595/100‑1 (MQ=255)
                                                         gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                           <  1:22393/100‑1 (MQ=255)
                                                         gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                           <  2:4058298/100‑1 (MQ=255)
                                                          aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                          >  2:1856065/1‑100 (MQ=255)
                                                               cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                     <  2:27875/100‑1 (MQ=255)
                                                               cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                     >  2:710095/1‑100 (MQ=255)
                                                                aTACAAGTCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCCAAAGGGGCGAgggg                                    <  1:2364424/100‑2 (MQ=255)
                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggtt                                    >  2:258520/1‑98 (MQ=255)
                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                    >  2:3320589/1‑100 (MQ=255)
                                                                   cAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                 <  1:4033028/100‑1 (MQ=255)
                                                                    aaGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                <  1:1195018/100‑1 (MQ=255)
                                                                        cACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                            >  1:1877273/1‑100 (MQ=255)
                                                                        cACCTTAATTAACGACATACAGCCTCCCGTGAAGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                            >  2:1896013/1‑100 (MQ=255)
                                                                         aCCTTAATTAACGATATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                           >  2:3440430/1‑100 (MQ=255)
                                                                         aCCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                           >  2:2522687/1‑100 (MQ=255)
                                                                          ccTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                          <  1:2434347/100‑1 (MQ=255)
                                                                           cTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                         >  2:1001681/1‑100 (MQ=255)
                                                                                ttaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                    <  2:127113/100‑1 (MQ=255)
                                                                                ttaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                    >  1:3244267/1‑100 (MQ=255)
                                                                                 taaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                   <  2:4220736/100‑1 (MQ=255)
                                                                                 taaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                   <  1:1011625/100‑1 (MQ=255)
                                                                                 ggtttgacatcAGCCTCACGGGAGGCTGTATGTCAATCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                    >  2:3106904/11‑100 (MQ=21)
                                                                                  gtttgacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                   >  2:4094733/10‑100 (MQ=21)
                                                                                  gtttgacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                   >  2:4134152/10‑100 (MQ=25)
                                                                                  aaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc                  <  2:1206904/100‑1 (MQ=255)
                                                                                  aaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc                  >  2:238006/1‑100 (MQ=255)
                                                                                   aCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc                 >  2:3704074/1‑100 (MQ=255)
                                                                                   aCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc                 >  2:2506053/1‑100 (MQ=255)
                                                                                    agACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCCCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc                  >  1:1230970/2‑98 (MQ=255)
                                                                                      gacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc               >  2:1384184/6‑100 (MQ=25)
                                                                                      gacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc               >  2:3800807/6‑100 (MQ=25)
                                                                                      aCATACCCCCTCCCGTGAGGCTCATGTCAAACCCAGAAACTGAACCGAGCGGGCCCCGCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc               <  1:2083594/99‑1 (MQ=25)
                                                                                      aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAAGGTCCCCTCGCCCCTTTGGGGAGAAGGTTCGGGAGCGGGGAACACCCGc               >  1:4390811/1‑99 (MQ=38)
                                                                                       cATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc             >  2:2689133/1‑100 (MQ=255)
                                                                                        catcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc             >  1:3638897/4‑100 (MQ=25)
                                                                                        catcAGCCTCACGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGGGGGGAAAACCCg                >  2:325112/4‑96 (MQ=37)
                                                                                         tACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc           >  1:4150890/1‑100 (MQ=255)
                                                                                          aCAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCt          >  1:1180134/1‑100 (MQ=255)
                                                                                             gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt       >  1:1558603/1‑100 (MQ=255)
                                                                                             gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCACCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt       >  2:418479/1‑100 (MQ=255)
                                                                                             gCCTCACGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt       >  1:1711088/1‑100 (MQ=255)
                                                                                              acTCCCGTGCGGCTGATGTCAAACACAAAAACTGCACCGGACGGTCCCCTTGCCCCTTTGTGGAGATGGGGATGTATAGGAGAAAACCCGCACTCTAttt      >  1:361734/2‑100 (MQ=255)
                                                                                                  ccGTGGGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg  <  1:1856058/100‑1 (MQ=255)
                                                                                                  ccGTGCGGCTGATGTCAAACCCAAAAACGGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa   <  1:3144389/99‑1 (MQ=255)
                                                                                                  ccGGGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg  <  2:2132568/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
CGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAG  >  NC_000913/4070312‑4070509

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: