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.

CGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAG  >  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)
          tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                          <  2:1774053/100‑1 (MQ=255)
          tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                          <  1:4295089/100‑1 (MQ=255)
          tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑a                                                                                          <  1:2563599/100‑1 (MQ=255)
           cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                         <  1:2058192/100‑1 (MQ=255)
            ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                        >  2:1815326/1‑100 (MQ=255)
            ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgtat                                                                                         <  2:2178009/100‑4 (MQ=255)
              ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAT‑CAGCCTCACGGGAGGCTGTATGTc                                                                                      <  1:3367177/100‑1 (MQ=38)
                gcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                    >  2:469994/1‑100 (MQ=255)
                 cgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTGTATGTcgt                                                                                    <  1:74599/100‑3 (MQ=38)
                  ggCAGGTGAAACGGTCTGTCCACTGCAATGGAGATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAc                                                                                  <  1:3767311/100‑1 (MQ=255)
                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAc                                                                                  <  1:2280162/100‑1 (MQ=255)
                   gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                 <  2:3244464/100‑1 (MQ=255)
                   gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                 >  1:2218207/1‑100 (MQ=255)
                    cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                 >  1:1030444/1‑99 (MQ=255)
                    cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAT‑CAGCCTCACGGGAGGCTGTATGTcgttaat                                                                                <  2:233100/100‑7 (MQ=25)
                     aaGTGAAACGGTCTGTACACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGACAAAccc                                                                                <  1:1746484/99‑1 (MQ=255)
                         gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                           >  2:2250879/1‑100 (MQ=255)
                         gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                           >  1:252615/1‑100 (MQ=255)
                         gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                           >  2:1064829/1‑100 (MQ=255)
                           aaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACt                                                                         >  2:426017/1‑100 (MQ=255)
                              ggTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCa                                                                      >  2:2576080/1‑100 (MQ=255)
                               gTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCAc                                                                     <  1:3363898/100‑1 (MQ=255)
                                 cTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCg                                                                   <  2:3430045/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCgg                                                                  <  2:4051041/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCgg                                                                  <  1:679247/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCgg                                                                  <  1:679087/100‑1 (MQ=255)
                                  tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCgg                                                                  >  2:553539/1‑100 (MQ=255)
                                   gTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGa                                                                 >  1:2802947/1‑100 (MQ=255)
                                    tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGAc                                                                <  2:2691682/100‑1 (MQ=255)
                                      cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACgg                                                              >  1:829044/1‑100 (MQ=255)
                                         tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           <  2:872099/100‑1 (MQ=255)
                                         tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           <  2:862687/100‑1 (MQ=255)
                                         tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           <  2:670394/100‑1 (MQ=255)
                                          gCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTccc                                                          >  2:2592378/1‑100 (MQ=255)
                                          gCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTccc                                                          >  1:857486/1‑100 (MQ=255)
                                          gCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                           >  1:3083273/1‑99 (MQ=255)
                                              tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCg                                                      >  1:1323405/1‑100 (MQ=255)
                                               gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                     >  2:1048499/1‑100 (MQ=255)
                                               gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                     >  2:3931366/1‑100 (MQ=255)
                                               gCATATTGGTGATGTGCATACAAGGCACCTTTATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGc                                                     >  2:2967286/1‑100 (MQ=255)
                                                cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                                    <  1:3545186/100‑1 (MQ=255)
                                                cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                                    <  2:3834122/100‑1 (MQ=255)
                                                cATATTGGTGATTTGCATACAAGGCACCGTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcc                                                    <  1:1548375/100‑1 (MQ=255)
                                                   atTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  2:2086122/1‑100 (MQ=255)
                                                    ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  1:286827/1‑99 (MQ=255)
                                                    ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                 >  2:3012459/1‑99 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              >  1:3015432/1‑100 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              >  1:949068/1‑100 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              <  2:1582188/100‑1 (MQ=255)
                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg                                              >  1:142880/1‑100 (MQ=255)
                                                       gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                             <  1:3171715/100‑1 (MQ=255)
                                                       gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                             >  1:609712/1‑100 (MQ=255)
                                                        tGATCTGCATACAAGGCACCTTAATTAACCACATACAGCCGCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTggg                                            <  2:933595/100‑1 (MQ=255)
                                                         gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                           <  1:22393/100‑1 (MQ=255)
                                                         gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                           <  2:4058298/100‑1 (MQ=255)
                                                          aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                          >  2:1856065/1‑100 (MQ=255)
                                                               cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                     <  2:27875/100‑1 (MQ=255)
                                                               cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                     >  2:710095/1‑100 (MQ=255)
                                                                aTACAAGTCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCCAAAGGGGCGAgggg                                    <  1:2364424/100‑2 (MQ=255)
                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggtt                                    >  2:258520/1‑98 (MQ=255)
                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt                                    >  2:3320589/1‑100 (MQ=255)
                                                                   cAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                 <  1:4033028/100‑1 (MQ=255)
                                                                    aaGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                <  1:1195018/100‑1 (MQ=255)
                                                                        cACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                            >  1:1877273/1‑100 (MQ=255)
                                                                        cACCTTAATTAACGACATACAGCCTCCCGTGAAGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga                            >  2:1896013/1‑100 (MQ=255)
                                                                         aCCTTAATTAACGATATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                           >  2:3440430/1‑100 (MQ=255)
                                                                         aCCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag                           >  2:2522687/1‑100 (MQ=255)
                                                                          ccTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                          <  1:2434347/100‑1 (MQ=255)
                                                                           cTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                         >  2:1001681/1‑100 (MQ=255)
                                                                                ttaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                    <  2:127113/100‑1 (MQ=255)
                                                                                ttaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                    >  1:3244267/1‑100 (MQ=255)
                                                                                 taaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                   <  2:4220736/100‑1 (MQ=255)
                                                                                 taaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                   <  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)
                                                                                  aaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc                  <  2:1206904/100‑1 (MQ=255)
                                                                                  aaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc                  >  2:238006/1‑100 (MQ=255)
                                                                                   aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc                 >  2:3704074/1‑100 (MQ=255)
                                                                                   aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc                 >  2:2506053/1‑100 (MQ=255)
                                                                                    agACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCCCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc                  >  1:1230970/2‑98 (MQ=255)
                                                                                      gacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc               >  2:1384184/6‑100 (MQ=25)
                                                                                      gacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc               >  2:3800807/6‑100 (MQ=25)
                                                                                      aCATACCCCCTCCCGTGAGGCTC‑ATGTCAAACCCAGAAACTGAACCGAGCGGGCCCCGCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGc               <  1:2083594/99‑1 (MQ=25)
                                                                                      aCATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGAAGGTCCCCTCGCCCCTTTGGGGAGAAGGTTCGGGAGCGGGGAACACCCGc               >  1:4390811/1‑99 (MQ=38)
                                                                                       cATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc             >  2:2689133/1‑100 (MQ=255)
                                                                                        catcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc             >  1:3638897/4‑100 (MQ=25)
                                                                                        catcAGCCTCACGGGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGGGGGGAAAACCCg                >  2:325112/4‑96 (MQ=37)
                                                                                         tACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc           >  1:4150890/1‑100 (MQ=255)
                                                                                          aCAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCt          >  1:1180134/1‑100 (MQ=255)
                                                                                             gCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt       >  1:1558603/1‑100 (MQ=255)
                                                                                             gCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCACCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt       >  2:418479/1‑100 (MQ=255)
                                                                                             gCCTCACGGGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt       >  1:1711088/1‑100 (MQ=255)
                                                                                              acTCCCGTGCGGCTG‑ATGTCAAACACAAAAACTGCACCGGACGGTCCCCTTGCCCCTTTGTGGAGATGGGGATGTATAGGAGAAAACCCGCACTCTAttt      >  1:361734/2‑100 (MQ=255)
                                                                                                  ccGTGGGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg  <  1:1856058/100‑1 (MQ=255)
                                                                                                  ccGTGCGGCTG‑ATGTCAAACCCAAAAACGGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa   <  1:3144389/99‑1 (MQ=255)
                                                                                                  ccGGGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg  <  2:2132568/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
CGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAG  >  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: ‑