Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,241,409 T→C Q87R (CAG→CGG)  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,241,4090TC98.7% 285.1 / ‑4.1 78Q87R (CAG→CGG) galSgalactose‑ and fucose‑inducible galactose regulon transcriptional isorepressor; mgl operon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (45/32);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

ATCCATAAATTGCGCCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCTGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCTCACTCACGGCTT  >  NC_000913/2241261‑2241556
                                                                                                                                                    |                                                                                                                                                   
aTCCATAAATTTCGCCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACTCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTCTTGCCGATTAGCCCGTATTTCTGAtgccgc                                                                                                                                                   <  1:756195/151‑1 (MQ=255)
   cATAAATTGCGCCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGa                                                                                                                                                <  2:157611/151‑1 (MQ=255)
      aaaTTGCGCCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCg                                                                                                                                             >  1:876064/1‑151 (MQ=255)
        aTTGCGCCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGAc                                                                                                                                           >  2:1021068/1‑151 (MQ=255)
           gcgcCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGAcc                                                                                                                                          >  1:743693/1‑149 (MQ=255)
           gcgcCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGAcc                                                                                                                                          >  2:1108426/1‑149 (MQ=255)
            cgcCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGAcc                                                                                                                                          >  1:994629/1‑148 (MQ=255)
             gcCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGa                                                                                                                                       <  1:55818/150‑1 (MQ=255)
             gcCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGAt                                                                                                                                      <  1:1283577/151‑1 (MQ=255)
              ccAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGAcc                                                                                                                                          >  2:499354/1‑146 (MQ=255)
              ccAGTTCATCGACACTCAATGCTTTTGAGCGAACCATCAACGCATTACTACGCTGGCGAATTAACACCTTAATGGCGGGACGCGCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTACTTCTGTTTCCGCTGAGCGAcc                                                                                                                                          >  1:627783/1‑146 (MQ=255)
                   tCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCAcc                                                                                                                                 >  2:351085/1‑150 (MQ=255)
                     aTCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCt                                                                                                                              <  1:134325/151‑1 (MQ=255)
                     aTCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAAAACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCt                                                                                                                              >  1:706827/1‑151 (MQ=255)
                          cACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTAc                                                                                                                         >  1:373748/1‑151 (MQ=255)
                          cACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTAc                                                                                                                         <  2:403410/151‑1 (MQ=255)
                            cTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTAcc                                                                                                                        >  1:1174941/1‑150 (MQ=255)
                              cAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCa                                                                                                                       <  2:1197739/149‑1 (MQ=255)
                              cAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTAAATGGCGTAAAGCTCTTTTTCCGCTTCATGAGAGCACTTGCCGATTAGCACGTATTTCTGATGCCGCAGAGCGACCAGATCCACCGCTTTTACCAgc                                                                                                                     >  1:763086/1‑151 (MQ=255)
                                aTGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAgcgc                                                                                                                   <  1:917619/151‑1 (MQ=255)
                                 tGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCTGCTGAGCGACCAGATCCACCGCTTTTACCAgcgc                                                                                                                   >  1:825313/1‑150 (MQ=255)
                                  gCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAgcgcgc                                                                                                                 <  1:192300/151‑1 (MQ=255)
                                  gCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAgcgcg                                                                                                                  <  1:639346/150‑1 (MQ=255)
                                     tttGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGcc                                                                                                                >  1:1082177/1‑149 (MQ=255)
                                      ttGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGa                                                                                                              >  2:1282050/1‑150 (MQ=255)
                                      ttGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGa                                                                                                              >  1:1014503/1‑150 (MQ=255)
                                       tGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGaa                                                                                                             <  2:1121943/150‑1 (MQ=255)
                                           tGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACg                                                                                                        >  1:159608/1‑151 (MQ=255)
                                              aCAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCAt                                                                                                     <  2:418142/151‑1 (MQ=255)
                                                    aaCGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGaaa                                                                                               <  1:104330/151‑1 (MQ=255)
                                                      cGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAg                                                                                                  <  1:35162/146‑1 (MQ=255)
                                                       gCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGt                                                                                            >  2:761403/1‑151 (MQ=255)
                                                          ttACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTcc                                                                                          >  2:957247/1‑150 (MQ=255)
                                                          ttACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTcc                                                                                          >  2:662099/1‑150 (MQ=255)
                                                          ttACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTcc                                                                                          >  1:713948/1‑150 (MQ=255)
                                                           tACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCAt                                                                                        >  2:964686/1‑151 (MQ=255)
                                                              aaCGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATcac                                                                                     >  1:765553/1‑151 (MQ=255)
                                                               aCGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATcacc                                                                                    <  2:30380/151‑1 (MQ=255)
                                                               aCGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATcacc                                                                                    >  2:373655/1‑151 (MQ=255)
                                                                   tGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATcaccacc                                                                                 >  2:1272025/1‑150 (MQ=255)
                                                                      cGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGc                                                                             >  1:1249546/1‑151 (MQ=255)
                                                                           tAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCCCCACGCCAAt                                                                         <  2:1103373/150‑1 (MQ=255)
                                                                            aaCACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGTCCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCATAAACGTCCATCACCACCACGCCAATg                                                                        <  1:1018715/150‑1 (MQ=255)
                                                                                ccTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTCGCACGTATTTCTGATGCCGATGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGgttcc                                                                   >  1:914414/1‑148 (MQ=255)
                                                                                 cTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCg                                                                  >  1:780974/1‑151 (MQ=255)
                                                                                   cAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCt                                                                >  2:645751/1‑151 (MQ=255)
                                                                                       ggCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAAcc                                                            >  2:321036/1‑151 (MQ=255)
                                                                                        gCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCTACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCt                                                           <  1:82386/151‑1 (MQ=255)
                                                                                          gTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTg                                                          >  2:1078327/1‑150 (MQ=255)
                                                                                             aCGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGtt                                                      <  1:761367/151‑1 (MQ=255)
                                                                                              cGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGtt                                                      <  2:986644/150‑1 (MQ=255)
                                                                                                 tctTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGcc                                                   >  2:1248616/1‑150 (MQ=255)
                                                                                                 tctTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGcc                                                   >  2:620965/1‑150 (MQ=255)
                                                                                                 tctTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGcc                                                   >  2:897950/1‑150 (MQ=255)
                                                                                                    ttttCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAgc                                                <  1:269719/150‑1 (MQ=255)
                                                                                                        ccGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCt                                            <  1:1273795/150‑1 (MQ=255)
                                                                                                        ccGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTg                                           <  1:1021030/151‑1 (MQ=255)
                                                                                                          gCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTggg                                         >  2:854517/1‑151 (MQ=255)
                                                                                                          gCTTCATGATAGCTATAGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTg                                           >  1:623578/1‑149 (MQ=255)
                                                                                                              cATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCAt                                      <  1:760925/150‑1 (MQ=255)
                                                                                                               aTGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTg                                    >  2:7958/1‑151 (MQ=255)
                                                                                                               aTGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTg                                    >  2:1249970/1‑151 (MQ=255)
                                                                                                                tGATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTgg                                   >  2:986878/1‑151 (MQ=255)
                                                                                                                 gATAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTgg                                   <  1:286725/150‑1 (MQ=255)
                                                                                                                  aTAGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCg                                 >  1:1206335/1‑151 (MQ=255)
                                                                                                                    aGCTATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGtt                               >  1:100388/1‑151 (MQ=255)
                                                                                                                       tATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGtt                               >  1:518224/1‑148 (MQ=255)
                                                                                                                       tATTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTgg                            >  1:252258/1‑151 (MQ=255)
                                                                                                                        aTTGCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTgg                            <  2:143997/150‑1 (MQ=255)
                                                                                                                           gCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTTTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGa                        <  1:431726/151‑1 (MQ=255)
                                                                                                                           gCCGATTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGa                        <  1:259286/151‑1 (MQ=255)
                                                                                                                               aTTAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAAt                    >  1:885433/1‑151 (MQ=255)
                                                                                                                                 tAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATc                   <  1:5998/150‑1 (MQ=255)
                                                                                                                                 tAGCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATc                   <  2:373775/150‑1 (MQ=255)
                                                                                                                                   gCACGTATTTCTGATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAg                <  1:847005/151‑1 (MQ=255)
                                                                                                                                               gATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCtcactcac       <  2:438156/148‑1 (MQ=255)
                                                                                                                                               gATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCTCACTCACgg     <  1:1081731/150‑1 (MQ=255)
                                                                                                                                               gATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCTCACTCACGGc    >  2:1175447/1‑151 (MQ=255)
                                                                                                                                               gATGCCGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCTCACTCACGGc    >  2:658723/1‑151 (MQ=255)
                                                                                                                                                 tgccgcTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCTCACTCACGGCtt  >  2:528644/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCCATAAATTGCGCCAGTTCATCGTCACTCAATGCTTTTGAGTGAACAATCAACGCATTACAACGCTGGCGAATTAACACCTCAATGGCGTGACGCTCTTTTTCCGCTTCATGATAGCTATTGCCGATTAGCACGTATTTCTGATGCTGCTGAGCGACCAGATCCACCGCTTTTACCAGCGCGCCGAAAAACGCATCAGAAACGTCCATCACCACCACGCCAATGGTGTCGCTAACCTGAGTTGCCAGCGCCTGGGCATTGGCGTTTGGCCGATAATCCAGCTCACTCACGGCTT  >  NC_000913/2241261‑2241556

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