Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,939,333 G→C 15.1% intergenic (+23/‑35) fucU → / → fucR L‑fucose mutarotase/l‑fucose operon activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,939,3330GC15.1% 465.5 / 67.4 165intergenic (+23/‑35)fucU/fucRL‑fucose mutarotase/l‑fucose operon activator
Reads supporting (aligned to +/- strand):  ref base G (65/75);  new base C (2/23);  total (67/98)
Fisher's exact test for biased strand distribution p-value = 2.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 166 total aligned reads displayed.

AGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGA  >  NC_000913/2939235‑2939432
                                                                                                  |                                                                                                   
 gCCTTTGCGATCGTTATCACAGGGGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                   <  1:2623835/100‑1 (MQ=255)
 gCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                   >  1:2469606/1‑100 (MQ=255)
  ccTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAcg                                                                                                  >  2:286299/1‑100 (MQ=255)
       gCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggc                                                                                             >  1:1934582/1‑100 (MQ=255)
        cGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcg                                                                                            >  2:2837585/1‑100 (MQ=255)
        cGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcg                                                                                            >  1:847012/1‑100 (MQ=255)
        cGATCGTTATCACAGGAGAACGAGCGAAGTACGGGAATATTCTATTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcg                                                                                            <  1:1698538/100‑1 (MQ=255)
         gATCGTTATCACAGGCGAACGCGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcgg                                                                                           <  2:1482820/100‑1 (MQ=255)
         gATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgt                                                                                           <  2:71839/100‑2 (MQ=255)
         gATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcgg                                                                                           <  2:2669284/100‑1 (MQ=255)
         gATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggcg                                                                                            <  2:2379910/99‑1 (MQ=255)
          gctgatatggctcaacagcatgtctactatcgcttgcaggcgtgccgctttcatatttttcaacccatcattcgataacCGCCCGCATGACGCGGGCGGt                                                                                          >  1:1532813‑M1/80‑100 (MQ=255)
          aTCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGt                                                                                          <  2:1448782/100‑1 (MQ=255)
          aTCGTTATCACAGGCGAACGAGCGAAGTAAGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGt                                                                                          <  2:978076/100‑1 (MQ=255)
           tCGTTATCACAGGCGAACGAGCGGAGTACGGGAATATTCTGTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGt                                                                                          <  1:1092372/99‑1 (MQ=255)
            cGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtacc                                                                                        <  2:1357341/100‑5 (MQ=255)
             gTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtaccg                                                                                       <  2:2365547/100‑6 (MQ=255)
             gTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtacc                                                                                        <  1:1323253/99‑5 (MQ=255)
             gTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTAt                                                                                       >  2:1872751/1‑100 (MQ=255)
             gTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTAAAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtacc                                                                                        <  1:2248186/99‑5 (MQ=255)
              ttATCACAGGCGCACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                      >  2:1906316/1‑100 (MQ=255)
              ttATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAACAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtac                                                                                         <  2:2161917/97‑4 (MQ=255)
              ttATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATc                                                                                      <  2:2020694/100‑1 (MQ=255)
                aTCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGa                                                                                    >  1:2572538/1‑100 (MQ=255)
                 tCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGaa                                                                                   <  1:728755/100‑1 (MQ=255)
                  cacaGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAat                                                                                  >  2:91754/1‑100 (MQ=255)
                  cacaGGCGAACGAGCGAAGTACGGGAATAATATTATAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtaccggtatg                                                                                  <  1:2791247/100‑11 (MQ=25)
                      ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                              >  1:502429/1‑100 (MQ=255)
                      ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                              <  1:2011565/100‑1 (MQ=255)
                      ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                              >  2:2011572/1‑100 (MQ=255)
                         gAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                            <  1:2703157/99‑1 (MQ=255)
                         gAACGAGCGAAGTACGGGAATATTAAAATAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATgggg                                                                           <  1:175208/100‑1 (MQ=255)
                         gAAAGAGCGAAGTACGGGAAGAAAAAAAAAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATgggg                                                                           <  1:174880/100‑1 (MQ=255)
                              aGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATATCATACCGGTACGCCCGCATGACGCGGGAGGTTATCGAATGATGGGGTGaa                                                                       >  1:1613057/1‑99 (MQ=255)
                               gCGAAGTACGGGAAAATTATTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGaaaa                                                                     <  1:221375/100‑1 (MQ=255)
                                cGAAGTACGGGAATATTCTTTGAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATCACACGGGCGGTTAGCGACTGATGGGGTgca                                                                       >  1:2656527/1‑95 (MQ=255)
                                      aCGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaa                                                              >  1:2815376/1‑100 (MQ=255)
                                      aCGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaa                                                              <  1:2086381/100‑1 (MQ=255)
                                       cGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                             <  1:527834/100‑1 (MQ=255)
                                       cGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                             >  2:2672613/1‑100 (MQ=255)
                                         ggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTTATACCGGTACGACCGCATGACGCGGGCGGTTATCGAATGATGGGGCGAAAAATATGAAAGc                                                           >  1:1814828/1‑100 (MQ=255)
                                         ggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                           <  1:111980/100‑1 (MQ=255)
                                          gAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCg                                                          >  2:2552952/1‑100 (MQ=255)
                                           aaTATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCgg                                                         <  1:1812464/100‑1 (MQ=255)
                                            atatTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGc                                                        <  2:2611803/100‑1 (MQ=255)
                                              atTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCgg                                                         >  1:2309553/1‑97 (MQ=255)
                                              atTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCAc                                                      <  1:2162522/100‑1 (MQ=255)
                                               ttCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACg                                                     <  2:166923/100‑1 (MQ=255)
                                                tCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGc                                                    <  2:2158738/100‑1 (MQ=255)
                                                  ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCa                                                  <  1:91753/100‑1 (MQ=255)
                                                  ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAAAATGAAAGCGGCACGCCa                                                  >  2:578525/1‑100 (MQ=255)
                                                  ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACACGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCa                                                  >  1:421743/1‑100 (MQ=255)
                                                    ttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGc                                                <  1:2001860/100‑1 (MQ=255)
                                                     tAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCa                                               >  1:2595196/1‑100 (MQ=255)
                                                       aaaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGCGAAAAATATGAAAGCGGCACGCCAGCAAg                                             >  1:1588639/1‑100 (MQ=255)
                                                        aaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGc                                            >  2:1139184/1‑100 (MQ=255)
                                                        aaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGc                                            <  1:1107584/100‑1 (MQ=255)
                                                        aaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGc                                            <  1:2117113/100‑1 (MQ=255)
                                                           aaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGAt                                         <  2:2794002/100‑1 (MQ=255)
                                                              gggTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGt                                      >  1:1281503/1‑100 (MQ=255)
                                                              gggTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGt                                      <  1:1026542/100‑1 (MQ=255)
                                                                gTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCg                                    >  1:1174216/1‑100 (MQ=255)
                                                                  aaCACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAGGCGGCACGCCAGCAAGCGATAGTCGAc                                  <  2:2479672/100‑1 (MQ=255)
                                                                  aaCACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAc                                  <  2:534742/100‑1 (MQ=255)
                                                                    cacCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCgtcaa                                >  1:2289518/1‑96 (MQ=255)
                                                                    cacCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                >  2:2064015/1‑100 (MQ=255)
                                                                      ccGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgc                              >  2:1992616/1‑100 (MQ=255)
                                                                      ccGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgc                              <  2:1296916/100‑1 (MQ=255)
                                                                        gTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctg                            <  2:1027716/100‑1 (MQ=255)
                                                                        gTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCACCAAGCGATAGTCGACctgct                             >  1:973727/1‑99 (MQ=255)
                                                                         tAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgc                           <  1:1639013/100‑1 (MQ=255)
                                                                          aaTCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgct                          <  2:1125012/100‑1 (MQ=255)
                                                                             ctcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGaa                       >  1:411201/1‑100 (MQ=255)
                                                                              tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAc                      <  2:693575/100‑1 (MQ=255)
                                                                               attcgataacCGCCCGCGTCATGCGGGCGGTTATCGAAGGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAcc                     <  1:1787785‑M1/90‑1 (MQ=255)
                                                                                aTACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCa                    <  2:2571572/100‑1 (MQ=255)
                                                                                 tACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCAt                   >  1:1342633/1‑100 (MQ=255)
                                                                                   nataacCGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGGACCATAc                 <  1:1343098‑M1/94‑1 (MQ=255)
                                                                                   gataacCGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAc                 <  1:1543187‑M1/94‑1 (MQ=255)
                                                                                   ccGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAc                 <  2:2323108/100‑1 (MQ=255)
                                                                                   ccGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAc                 >  2:1495525/1‑100 (MQ=255)
                                                                                    cGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAcc                <  1:2672602/100‑1 (MQ=255)
                                                                                    ataacCGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAcc                <  1:2068747‑M1/95‑1 (MQ=255)
                                                                                      aacCGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAg              <  2:376315‑M1/97‑1 (MQ=255)
                                                                                      aacCGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAg              <  2:1200237‑M1/97‑1 (MQ=255)
                                                                                       acCGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGc             <  1:2847417‑M1/98‑1 (MQ=255)
                                                                                        aCGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGcc            <  1:2433968/100‑1 (MQ=255)
                                                                                         CGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCt           <  1:60277‑M1/100‑1 (MQ=255)
                                                                                         CGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCt           <  2:398192‑M1/100‑1 (MQ=255)
                                                                                          GCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGACTg          <  1:1162687‑M1/100‑1 (MQ=255)
                                                                                          GCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTg          >  2:686818‑M1/1‑100 (MQ=255)
                                                                                           CCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGa         >  1:286297‑M1/1‑100 (MQ=255)
                                                                                           CCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAATCGACCTGCTGCTGAACCATACCAGCCTGa         <  1:336273‑M1/100‑1 (MQ=255)
                                                                                            CCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGAc        <  1:1048414‑M1/100‑1 (MQ=255)
                                                                                             CGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGAcc       >  1:2445598‑M1/1‑100 (MQ=255)
                                                                                              GCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCa      <  1:2336198‑M1/100‑1 (MQ=255)
                                                                                                 tGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACgg   <  1:322290/100‑1 (MQ=255)
                                                                                                 tGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACgg   >  1:2482317/1‑100 (MQ=255)
                                                                                                 tGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACgg   >  2:865033/1‑100 (MQ=255)
                                                                                                  gACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGa  >  1:278032/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGA  >  NC_000913/2939235‑2939432

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: