Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,240,754 (C)5→6 coding (915/1041 nt) galS ← galactose‑ and fucose‑inducible galactose regulon transcriptional isorepressor, mgl operon transcriptional repressor, autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,240,7491.C96.7% 208.2 / ‑0.6 61coding (920/1041 nt)galSgalactose‑ and fucose‑inducible galactose regulon transcriptional isorepressor, mgl operon transcriptional repressor, autorepressor
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (27/31);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGA  >  NC_000913/2240614‑2240880
                                                                                                                                         |                                                                                                                                     
caGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGccc                                                                                                                                    <  1:1276237/139‑1 (MQ=255)
     tGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctg                                                                                                                               <  2:955120/139‑4 (MQ=255)
      gAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              <  2:700588/139‑5 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCa                                                                                                                        >  2:1323892/1‑139 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCa                                                                                                                        >  2:800759/1‑139 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCa                                                                                                                        >  1:290032/1‑139 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              >  2:529750/1‑128 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              >  2:546941/1‑128 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              <  1:546941/132‑5 (MQ=255)
             gTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              <  1:529750/132‑5 (MQ=255)
                 cATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTc                                                                                                                   <  1:1323892/139‑1 (MQ=255)
                     gCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              >  2:1340823/1‑120 (MQ=255)
                     gCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              <  1:1340823/124‑5 (MQ=255)
                      ccTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTgg                                                                                                              >  2:1012796/1‑139 (MQ=255)
                              gttgAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              >  2:182695/1‑111 (MQ=255)
                              gttgAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGcccccctgc                                                                                                                              <  1:182695/115‑5 (MQ=255)
                               ttgAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCAt                                                                                                >  1:1452587/1‑139 (MQ=255)
                                      aGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCa                                                                                                 >  2:347233/1‑135 (MQ=255)
                                      aGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCa                                                                                                 <  1:347233/135‑1 (MQ=255)
                                      aGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCCAATTTCGCCATTGa                                                                                             <  2:290032/139‑1 (MQ=255)
                                         gATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGc                                                                                          >  1:141412/1‑139 (MQ=255)
                                           tCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                        <  1:800759/139‑1 (MQ=255)
                                            cgccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           >  2:1109095/1‑135 (MQ=255)
                                            cgccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           <  1:1109095/135‑1 (MQ=255)
                                             gccgctgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGTCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             <  2:916135/132‑1 (MQ=255)
                                             gccgccgcATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGa                                                                                             >  1:916135/1‑132 (MQ=255)
                                                     aTTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACg                                                                              <  2:141412/139‑1 (MQ=255)
                                                       tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCATGGCCAGTTCGTTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  1:976648/139‑1 (MQ=255)
                                                       tCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:1131943/139‑1 (MQ=255)
                                                         tGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACg                                                                          <  1:142525/139‑1 (MQ=255)
                                                                ttGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTa                                                                   >  2:1225472/1‑139 (MQ=255)
                                                                        gAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCggg                                                           <  2:989425/139‑1 (MQ=255)
                                                                         agtGGCGACGCACTAACGTCGGCATAAAACAGTGGCAGGCACGAGGATTAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTACGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGt                                                          <  1:890499/137‑1 (MQ=255)
                                                                          aTGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTc                                                         <  1:1388253/139‑1 (MQ=255)
                                                                                 cGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           <  2:1365943/98‑1 (MQ=255)
                                                                                 cGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAg                                                                                           >  1:1365943/1‑98 (MQ=255)
                                                                                 cGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTaa                                                  >  1:1264937/1‑139 (MQ=255)
                                                                                   cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCg                                                                                                                  >  2:603156/1‑74 (MQ=255)
                                                                                   cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCg                                                                                                                  <  1:603156/74‑1 (MQ=255)
                                                                                   cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACg                                                >  1:282651/1‑139 (MQ=255)
                                                                                   cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACg                                                >  2:849494/1‑139 (MQ=255)
                                                                                   cACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGACAACGCACGGTCGTTAATTGCGGGTCGGTGTAACg                                                >  1:740252/1‑139 (MQ=255)
                                                                                     cTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACggg                                              <  2:1478075/139‑1 (MQ=255)
                                                                                      tAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGc                                             <  1:1225472/139‑1 (MQ=255)
                                                                                      tAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGACCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCCATGGGATAACGCACGGTCGTTAATTGCGGGGCGGTGTAACggct                                             >  2:890499/1‑137 (MQ=255)
                                                                                       aaCGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCg                                            <  1:1012796/139‑1 (MQ=255)
                                                                                         cGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                          <  1:299098/139‑1 (MQ=255)
                                                                                         cGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                          <  1:127898/139‑1 (MQ=255)
                                                                                             ggCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGa                                      <  2:1264937/139‑1 (MQ=255)
                                                                                              gCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTaa                                                  >  1:1209204/1‑126 (MQ=255)
                                                                                              gCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTaa                                                  <  2:1209204/126‑1 (MQ=255)
                                                                                                aTAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAata                                   <  1:1068133/139‑1 (MQ=255)
                                                                                                         tggctggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACg                                                                          >  1:659988/1‑91 (MQ=255)
                                                                                                         tggctggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACg                                                                          <  2:659988/91‑1 (MQ=255)
                                                                                                          ggctggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:790985/88‑1 (MQ=255)
                                                                                                          ggctggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            >  1:790985/1‑88 (MQ=255)
                                                                                                             tggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCga                      <  1:849494/139‑1 (MQ=255)
                                                                                                              ggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCa                                                                                                                        <  1:447924/42‑1 (MQ=255)
                                                                                                              ggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCa                                                                                                                        >  2:447924/1‑42 (MQ=255)
                                                                                                              ggcACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCgat                     >  1:1334800/1‑139 (MQ=255)
                                                                                                                cACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            <  2:983514/82‑1 (MQ=255)
                                                                                                                cACGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCa                                                                            >  1:983514/1‑82 (MQ=255)
                                                                                                                 aCGAGGATCAATATTGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCgatgat                  >  2:1056898/1‑139 (MQ=255)
                                                                                                                          aaatttccttgctgcCCCCTGCAAAGCCCGTTTCGTTGGCAATTTCGCCAATTGAAGCAATGGGTAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGt        <  1:88590/130‑1 (MQ=255)
                                                                                                                              ttGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAAc     >  1:1387242/1‑139 (MQ=255)
                                                                                                                               tGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGc                                                                                                    <  2:301964/44‑1 (MQ=255)
                                                                                                                               tGCCTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGc                                                                                                    >  1:301964/1‑44 (MQ=255)
                                                                                                                                 ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                          <  2:1116765/99‑1 (MQ=255)
                                                                                                                                 ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGAt                                          >  1:1116765/1‑99 (MQ=255)
                                                                                                                                 ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCg                       <  1:962338/118‑1 (MQ=255)
                                                                                                                                 ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCg                       >  2:962338/1‑118 (MQ=255)
                                                                                                                                 ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGa  >  2:1118867/1‑139 (MQ=255)
                                                                                                                                 ccTGCTGCCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGa  >  2:128832/1‑139 (MQ=255)
                                                                                                                                   tgctgcCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                        <  1:652422/51‑1 (MQ=255)
                                                                                                                                   tgctgcCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCaa                                                                                        >  2:652422/1‑51 (MQ=255)
                                                                                                                                         |                                                                                                                                     
CAGATTGAAAGCGGTTACATCGCCTGATTTGTTGAGTTAGTGATCGCCGCCGCATTCTGGCGTGTTGCGACAGAATGGCGACGCACTAACGTCGGCATAAAACAGTGGCTGGCACGAGGATCAATATTGCCTGCTGCCCCCTGCAAGGCCAGTTCGGTGGCTAATTTCGCCATTGAAGCAATGGGATAACGCACGGTCGTTAATTGCGGGTCGGTGTAACGGGCGATGGGAATATCATCGAAACCGATGATTGAGAGATGTAACGGA  >  NC_000913/2240614‑2240880

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: