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)
            tCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGa                                                                                                                                     >  2:769766/1‑150 (MQ=255)
             cATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGaaa                                                                                                                                   <  2:1197522/151‑1 (MQ=255)
                  tttaatttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAAc                                                                                                                              >  2:861132/1‑151 (MQ=255)
                     aatttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGc                                                                                                                           >  2:656191/1‑151 (MQ=255)
                      atttaaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGcc                                                                                                                          >  1:512432/1‑151 (MQ=255)
                         taaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTg                                                                                                                        <  2:1008165/150‑1 (MQ=255)
                         taaAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTg                                                                                                                        <  1:420152/150‑1 (MQ=255)
                          aaaaCGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGAc                                                                                                                      <  1:14726/151‑1 (MQ=255)
                            aaCGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGt                                                                                                                    <  1:356836/151‑1 (MQ=255)
                              cGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTg                                                                                                                   >  1:621013/1‑150 (MQ=255)
                               gTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGaa                                                                                                                 <  2:462186/151‑1 (MQ=255)
                                 aagcgcgTTGCGCAATGCATCAGACTCTCAAGGTTATCTGCGCCCGCGACGCCTGCGGCTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAATACCTGAAAAAAACGGCCTGGCGTGAATc                                                                                                               <  2:1184571/148‑1 (MQ=255)
                                 aaTCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGTCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATc                                                                                                               <  1:958548/151‑1 (MQ=255)
                                       gTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCa                                                                                                          >  1:564561/1‑150 (MQ=255)
                                         tGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt                                                                                                       >  1:834406/1‑151 (MQ=255)
                                          gcgcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtt                                                                                                       <  1:253777/150‑1 (MQ=255)
                                           cgcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtttt                                                                                                     <  2:463944/151‑1 (MQ=255)
                                            gcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttt                                                                                                    <  1:298160/151‑1 (MQ=255)
                                            gcAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttt                                                                                                    <  1:162290/151‑1 (MQ=255)
                                             cAATGCATTAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAtttttt                                                                                                   >  2:69814/1‑151 (MQ=255)
                                              aaTGCATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAAttttttt                                                                                                  <  2:970698/151‑1 (MQ=255)
                                               aTGCATCCGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTc                                                                                                 >  1:562606/1‑151 (MQ=255)
                                                  cATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAgg                                                                                              >  2:590886/1‑151 (MQ=255)
                                                  cATCAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAgg                                                                                              <  1:1111785/151‑1 (MQ=255)
                                                     cAGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGttt                                                                                           <  1:874566/151‑1 (MQ=255)
                                                      aGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGtttt                                                                                          <  2:309853/151‑1 (MQ=255)
                                                      aGACTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATGTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGttt                                                                                           <  2:505082/151‑1 (MQ=255)
                                                        aCTCACAATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGc                                                                                        <  2:206043/151‑1 (MQ=255)
                                                           cacaATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGccc                                                                                      <  2:105581/150‑1 (MQ=255)
                                                           cacaATGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCg                                                                                     >  1:900267/1‑151 (MQ=255)
                                                               aTGTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCtt                                                                                  >  1:454123/1‑150 (MQ=255)
                                                                 gTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGt                                                                               >  2:848506/1‑151 (MQ=255)
                                                                 gTTATCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGt                                                                               >  1:882925/1‑151 (MQ=255)
                                                                     tCTGCGCCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCgg                                                                           >  2:769182/1‑151 (MQ=255)
                                                                        gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCg                                                                            <  1:353460/147‑1 (MQ=255)
                                                                        gcgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTa                                                                         <  1:1147819/150‑1 (MQ=255)
                                                                         cgcCCCCGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAAc                                                                       <  1:380613/151‑1 (MQ=255)
                                                                              ccGACTCCTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATcc                                                                  >  2:358034/1‑151 (MQ=255)
                                                                                    ccTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACt                                                            <  1:412375/151‑1 (MQ=255)
                                                                                     cTGCGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTc                                                           >  1:474734/1‑151 (MQ=255)
                                                                                        cGACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGt                                                        >  2:1100864/1‑151 (MQ=255)
                                                                                         gACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGta                                                       <  1:1226229/151‑1 (MQ=255)
                                                                                         gACTATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGt                                                        >  1:684783/1‑150 (MQ=255)
                                                                                            tATTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaat                                                     >  2:704651/1‑150 (MQ=255)
                                                                                             aTTTTTCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGtaat                                                     <  1:194029/149‑1 (MQ=255)
                                                                                                tttCTGGCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAatat                                                >  2:184509/1‑151 (MQ=255)
                                                                                                     ggCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTc                                             <  1:395907/149‑1 (MQ=255)
                                                                                                      gCTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAgt                                          <  2:974610/151‑1 (MQ=255)
                                                                                                       cTAACTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAgtg                                         >  1:247983/1‑151 (MQ=255)
                                                                                                          aCTCCGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTc                                      <  2:620317/151‑1 (MQ=255)
                                                                                                              cGTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTg                                  <  1:833733/151‑1 (MQ=255)
                                                                                                               gTCATAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTgc                                 >  1:1079648/1‑151 (MQ=255)
                                                                                                                   tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAc                             >  1:1131568/1‑151 (MQ=255)
                                                                                                                   tAACTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAc                             >  1:173333/1‑151 (MQ=255)
                                                                                                                    aaCTTGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcg                            >  1:1181904/1‑151 (MQ=255)
                                                                                                                        tGCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCAcgcg                          <  2:358907/149‑1 (MQ=255)
                                                                                                                         gCCCTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGct                        <  1:590994/150‑1 (MQ=255)
                                                                                                                          cccTCTACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGctct                      <  1:190748/151‑1 (MQ=255)
                                                                                                                               tACCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATAtt                 >  2:821299/1‑151 (MQ=255)
                                                                                                                                aCCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCtata                   <  1:144384/148‑1 (MQ=255)
                                                                                                                                aCCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTatat                  <  2:888953/149‑1 (MQ=255)
                                                                                                                                aCCGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATAtt                 <  1:710775/150‑1 (MQ=255)
                                                                                                                                  cGCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTAt              <  2:74086/151‑1 (MQ=255)
                                                                                                                                   gCTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATg             >  2:637251/1‑151 (MQ=255)
                                                                                                                                    cTTTGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATgg            >  1:432905/1‑151 (MQ=255)
                                                                                                                                       tGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCt          <  1:69820/150‑1 (MQ=255)
                                                                                                                                       tGCGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCt          <  1:564372/150‑1 (MQ=255)
                                                                                                                                         cGGCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTaa        >  1:1192664/1‑150 (MQ=255)
                                                                                                                                          ggCAAAACTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTaaa       >  2:926238/1‑150 (MQ=255)
                                                                                                                                             aaaaCTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAAc    >  2:837216/1‑150 (MQ=255)
                                                                                                                                              aaaCTCCAAAAAAAACCTGAAAAAAACGGCCTGACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACAt  <  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: