Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,905,639 (A)9→8 intergenic (‑42/+92) bglF ← / ← bglG fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component/transcriptional antiterminator of the bgl operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,905,6310A.98.6% 334.3 / ‑1.5 71intergenic (‑34/+100)bglF/bglGfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component/transcriptional antiterminator of the bgl operon
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (32/38);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACAT  >  NC_000913/3905482‑3905775
                                                                                                                                                     |                                                                                                                                                
tGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCa                                                                                                                                                  <  1:925404/150‑1 (MQ=255)
tGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCa                                                                                                                                                  >  2:923444/1‑150 (MQ=255)
  cgcTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaa                                                                                                                                               <  1:56553/151‑1 (MQ=255)
     tttgctttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaa                                                                                                                                             >  1:60736/1‑150 (MQ=255)
       tgctttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaac                                                                                                                                          >  1:702082/1‑150 (MQ=255)
        gctttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaac                                                                                                                                          <  1:934041/150‑2 (MQ=255)
        gctttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaac                                                                                                                                          <  2:275838/150‑2 (MQ=255)
          tttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaacctg                                                                                                                                       >  2:1081756/1‑147 (MQ=255)
           ttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaacctg                                                                                                                                       <  2:900864/150‑5 (MQ=255)
           ttCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCaaaaaaaacctg                                                                                                                                       >  2:622497/1‑146 (MQ=255)
            tCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGa                                                                                                                                     >  2:769766/1‑150 (MQ=255)
             cATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGaaa                                                                                                                                   <  2:1197522/151‑1 (MQ=255)
                  tttaatttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAAc                                                                                                                              >  2:861132/1‑151 (MQ=255)
                     aatttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGc                                                                                                                           >  2:656191/1‑151 (MQ=255)
                      atttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGcc                                                                                                                          >  1:512432/1‑151 (MQ=255)
                         taaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTg                                                                                                                        <  2:1008165/150‑1 (MQ=255)
                         taaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTg                                                                                                                        <  1:420152/150‑1 (MQ=255)
                          aaaaCGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGAc                                                                                                                      <  1:14726/151‑1 (MQ=255)
                            aaCGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGt                                                                                                                    <  1:356836/151‑1 (MQ=255)
                              cGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTg                                                                                                                   >  1:621013/1‑150 (MQ=255)
                               gTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGaa                                                                                                                 <  2:462186/151‑1 (MQ=255)
                                 aagcgcgTTGCGCAATGCATCAGACTCTCAAGGTTATCTGCGCCCGCGACGCCTGCGGCTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAATACCTGAAAAAAACGGCCTGGCGTGAATc                                                                                                               <  2:1184571/148‑1 (MQ=255)
                                 aaTCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGTCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATc                                                                                                               <  1:958548/151‑1 (MQ=255)
                                       gTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCa                                                                                                          >  1:564561/1‑150 (MQ=255)
                                         tGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt                                                                                                       >  1:834406/1‑151 (MQ=255)
                                          gcgcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt                                                                                                       <  1:253777/150‑1 (MQ=255)
                                           cgcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtttt                                                                                                     <  2:463944/151‑1 (MQ=255)
                                            gcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttt                                                                                                    <  1:298160/151‑1 (MQ=255)
                                            gcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttt                                                                                                    <  1:162290/151‑1 (MQ=255)
                                             cAATGCATTAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtttttt                                                                                                   >  2:69814/1‑151 (MQ=255)
                                              aaTGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttttt                                                                                                  <  2:970698/151‑1 (MQ=255)
                                               aTGCATCCGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTc                                                                                                 >  1:562606/1‑151 (MQ=255)
                                                  cATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAgg                                                                                              >  2:590886/1‑151 (MQ=255)
                                                  cATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAgg                                                                                              <  1:1111785/151‑1 (MQ=255)
                                                     cAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGttt                                                                                           <  1:874566/151‑1 (MQ=255)
                                                      aGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGtttt                                                                                          <  2:309853/151‑1 (MQ=255)
                                                      aGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATGTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGttt                                                                                           <  2:505082/151‑1 (MQ=255)
                                                        aCTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGc                                                                                        <  2:206043/151‑1 (MQ=255)
                                                           cacaATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGccc                                                                                      <  2:105581/150‑1 (MQ=255)
                                                           cacaATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCg                                                                                     >  1:900267/1‑151 (MQ=255)
                                                               aTGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCtt                                                                                  >  1:454123/1‑150 (MQ=255)
                                                                 gTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGt                                                                               >  2:848506/1‑151 (MQ=255)
                                                                 gTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGt                                                                               >  1:882925/1‑151 (MQ=255)
                                                                     tCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCgg                                                                           >  2:769182/1‑151 (MQ=255)
                                                                        gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCg                                                                            <  1:353460/147‑1 (MQ=255)
                                                                        gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTa                                                                         <  1:1147819/150‑1 (MQ=255)
                                                                         cgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAAc                                                                       <  1:380613/151‑1 (MQ=255)
                                                                              ccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATcc                                                                  >  2:358034/1‑151 (MQ=255)
                                                                                    ccTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                                            <  1:412375/151‑1 (MQ=255)
                                                                                     cTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTc                                                           >  1:474734/1‑151 (MQ=255)
                                                                                        cGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGt                                                        >  2:1100864/1‑151 (MQ=255)
                                                                                         gACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGta                                                       <  1:1226229/151‑1 (MQ=255)
                                                                                         gACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGt                                                        >  1:684783/1‑150 (MQ=255)
                                                                                            tATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaat                                                     >  2:704651/1‑150 (MQ=255)
                                                                                             aTTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaat                                                     <  1:194029/149‑1 (MQ=255)
                                                                                                tttCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAatat                                                >  2:184509/1‑151 (MQ=255)
                                                                                                     ggCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTc                                             <  1:395907/149‑1 (MQ=255)
                                                                                                      gCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAgt                                          <  2:974610/151‑1 (MQ=255)
                                                                                                       cTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAgtg                                         >  1:247983/1‑151 (MQ=255)
                                                                                                          aCTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTc                                      <  2:620317/151‑1 (MQ=255)
                                                                                                              cGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTg                                  <  1:833733/151‑1 (MQ=255)
                                                                                                               gTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTgc                                 >  1:1079648/1‑151 (MQ=255)
                                                                                                                   tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAc                             >  1:1131568/1‑151 (MQ=255)
                                                                                                                   tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAc                             >  1:173333/1‑151 (MQ=255)
                                                                                                                    aaCTTGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg                            >  1:1181904/1‑151 (MQ=255)
                                                                                                                        tGCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg                          <  2:358907/149‑1 (MQ=255)
                                                                                                                         gCCCTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGct                        <  1:590994/150‑1 (MQ=255)
                                                                                                                          cccTCTACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGctct                      <  1:190748/151‑1 (MQ=255)
                                                                                                                               tACCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATAtt                 >  2:821299/1‑151 (MQ=255)
                                                                                                                                aCCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCtata                   <  1:144384/148‑1 (MQ=255)
                                                                                                                                aCCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTatat                  <  2:888953/149‑1 (MQ=255)
                                                                                                                                aCCGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATAtt                 <  1:710775/150‑1 (MQ=255)
                                                                                                                                  cGCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTAt              <  2:74086/151‑1 (MQ=255)
                                                                                                                                   gCTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATg             >  2:637251/1‑151 (MQ=255)
                                                                                                                                    cTTTGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATgg            >  1:432905/1‑151 (MQ=255)
                                                                                                                                       tGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCt          <  1:69820/150‑1 (MQ=255)
                                                                                                                                       tGCGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCt          <  1:564372/150‑1 (MQ=255)
                                                                                                                                         cGGCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTaa        >  1:1192664/1‑150 (MQ=255)
                                                                                                                                          ggCAAAACTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTaaa       >  2:926238/1‑150 (MQ=255)
                                                                                                                                             aaaaCTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAAc    >  2:837216/1‑150 (MQ=255)
                                                                                                                                              aaaCTCC‑AAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACAt  <  2:979008/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACAT  >  NC_000913/3905482‑3905775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑