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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,939,3330GC 6.8% 173.0 / 2.8 65intergenic (+23/‑35)fucU/fucRL‑fucose mutarotase/l‑fucose operon activator
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base C (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGT  >  NC_000913/2939186‑2939475
                                                                                                                                                   |                                                                                                                                              
cTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGa                                                                                                                                               <  1:814492/149‑1 (MQ=255)
cTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGa                                                                                                                                               >  1:239488/1‑149 (MQ=255)
     atcatcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgg                                                                                                                                          >  1:720678/1‑149 (MQ=255)
       catcCTCATCAATCGTTTTGCGTTTTATGAACTGTCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACgcgggc                                                                                                                                        <  1:395976/149‑1 (MQ=255)
        atcCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACTCCCGCATGACgcgtgcg                                                                                                                                       >  1:1010689/1‑149 (MQ=255)
           cgaaTCAATCGTGTTGCGTGTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTATCACCGTAATCTCATACCGGTACGCCCGCGTCATgcgggcgtac                                                                                                                                    <  2:792885/146‑4 (MQ=255)
               tCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCg                                                                                                                                <  2:720678/149‑1 (MQ=255)
                 aaTCGTTTAGCGTTTTATGAACGTGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACTAGCGAAGTACGTGAATATTCTTTTAAAAAAATGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGTGCGGTTATCGaa                                                                                                                              <  2:927623/149‑1 (MQ=255)
                  aTCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACCAGCTAAGTACCGGAATATTCTTTTAAAAAAAGGGCTAACACCGTAATCTCAGACCGGTACGCCCGGATGACGCGGGCGGTTATCGAat                                                                                                                             >  2:987678/1‑149 (MQ=255)
                        ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATggg                                                                                                                       <  2:486834/149‑1 (MQ=255)
                           cGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTg                                                                                                                    >  2:308869/1‑149 (MQ=255)
                            gTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGa                                                                                                                   >  1:430382/1‑149 (MQ=255)
                            gTTTTATGAACGGGAGCAAAAAGCCTTTGCGATCGTTATCACAGGCTATCGAGCGAAGTACGTGAATATTCTTTTAAAATAAGGGGTAACACCGTAATCTCATACCGGTACGCCAGCATGACGCGGGCGGTTATCGAATGATGGGGGGa                                                                                                                   <  2:574473/149‑1 (MQ=255)
                             ttttATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGTGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGaa                                                                                                                  <  1:429186/149‑1 (MQ=255)
                                 aTGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAat                                                                                                              <  1:319398/149‑1 (MQ=255)
                                   gAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAatat                                                                                                            <  1:308869/149‑1 (MQ=255)
                                        ggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGACCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAg                                                                                                       <  2:799834/149‑1 (MQ=255)
                                          cgcAAAAAGCCTTTGAGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCg                                                                                                     <  2:828271/149‑1 (MQ=255)
                                             aaaaaGCCTTTGCGATCGTTATCACAGGCGCACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATAAAGGTGTGAAAAATATGAAAGCGGcc                                                                                                  >  1:198117/1‑148 (MQ=255)
                                             aaaaaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGcc                                                                                                  >  1:239373/1‑148 (MQ=255)
                                             aaaaaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCAATTAAAAAAAGCGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGCTTATCGAATGATGGGGTAAAAAATATGAAAGCGGCa                                                                                                  >  2:444065/1‑149 (MQ=255)
                                               aaaGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACg                                                                                                >  2:416576/1‑149 (MQ=255)
                                                   ccTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACTCCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAg                                                                                            <  1:813522/149‑1 (MQ=255)
                                                     tttGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCa                                                                                          >  2:363739/1‑149 (MQ=255)
                                                          gATCGTTATCACAGGCGAACGAGCTAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGa                                                                                     <  1:187324/149‑1 (MQ=255)
                                                          gATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGGGAAAAATATGAAAGCGGCACGCCAGCAAGCGa                                                                                     <  1:205566/149‑1 (MQ=255)
                                                                 ttCACAGGCGTAGGAGCTAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAATAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGa                                                                              <  2:29915/148‑1 (MQ=255)
                                                                  tCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGAc                                                                             <  2:621985/149‑1 (MQ=255)
                                                                     caGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCAAACCCGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGAGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGACGCCctg                                                                          >  1:781156/1‑149 (MQ=255)
                                                                      aGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgc                                                                         >  1:767100/1‑149 (MQ=255)
                                                                       ggCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgct                                                                        <  2:28705/149‑1 (MQ=255)
                                                                         cGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgc                                                                      <  1:1026634/149‑1 (MQ=255)
                                                                           aaCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgctg                                                                    >  1:176734/1‑149 (MQ=255)
                                                                               aGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTACTCTCATACCGGGACGCCCGCATGTCGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAAcc                                                                <  2:475436/149‑1 (MQ=255)
                                                                                gCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCa                                                               >  1:500694/1‑149 (MQ=255)
                                                                                  gAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTAAGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATATTCGACCTGCTGCTGAACCATa                                                             <  2:768456/149‑1 (MQ=255)
                                                                                       aCGGGAATATTCTTTTAAAAACTGGGGTAACACCGTAATCTCCTCCCGGTACGCCCGCATCACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGGCGACCTGCTGCTGAACCATACCAGc                                                        <  2:433545/149‑1 (MQ=255)
                                                                                       aCGGGAATATGCTTGAAAAACAAGGGGTAACACCGTTATCTCAGACCGGTACGCCCGCCTGACGCGGGCGGTTATCGAATGAGGGGGGGAAAAATATGAAAGCGGCACGACAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGc                                                        <  1:577964/149‑1 (MQ=255)
                                                                                          ggAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAACGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGACCCATGCCAGCCtc                                                     >  1:574325/1‑148 (MQ=255)
                                                                                            aaTATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGAc                                                   <  1:359658/149‑1 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGCAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCAGAACCATCCCAGCCTGCCCACGGaa                                            >  1:524203/1‑149 (MQ=255)
                                                                                                   ttttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGGGATCGAATGATGGGGGGAAAATTATGAAGGCGGCACGCCAGCCAGCGATAGTCGACCTGCTGCTGAACCCTCCCAGCCTGACCACGGat                                            >  2:944912/1‑148 (MQ=255)
                                                                                                    tttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATTAAAGCGGCACGCAAGAAAGCGATATTCGACCTGCTCCTGAACCATACCAGCCTGCCCACGCAAg                                           >  2:69450/1‑149 (MQ=255)
                                                                                                    tttAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAg                                           >  2:72328/1‑149 (MQ=255)
                                                                                                        aaaaaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctct                                       >  2:391815/1‑149 (MQ=255)
                                                                                                          aaaaGGGGTAACACCGTAATCTCATTCCGGTACGCCCGCATGACGCGGGCGGTTATCGCAGGATCGGGTGAAAAATATGATAGCGGCACGCCAGCTAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGCCCACGGAAGctcgct                                     >  1:619105/1‑149 (MQ=255)
                                                                                                          aaaaGGGGGAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctct                                     <  2:112155/149‑1 (MQ=255)
                                                                                                           aaagctctAACACCGTATGCTCTCACCGGTACGCCCGCATGACGCTGGCGGTTATCGAAAGATGGGGTGAAAAATATGAAAGCGTCACGCCAGCACGCTATACGCGACCTGCTGCTGAACCATCCCAGCCTGACCACGGAAGctctctc                                    <  2:318020/142‑1 (MQ=255)
                                                                                                           aaaGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGctctctc                                    <  2:446780/149‑1 (MQ=255)
                                                                                                            aaGTGGTAAAACCGGAATCTCATACCTGTACGCCCGCATGACGCGGGCTGTAATCGCATGATTGGCTGAAAAATATGAACGCGGCACGCAAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCGGACCACGGAAGctctctct                                   <  1:655148/149‑1 (MQ=255)
                                                                                                               gggTAACACCGTAATCTCATACCGGTACGCCCGCAAGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGaa                                <  1:340426/149‑1 (MQ=255)
                                                                                                                      acCGTAATCTCATACAGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTaa                         <  2:893454/149‑1 (MQ=255)
                                                                                                                            aTCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCa                   <  1:290247/149‑1 (MQ=255)
                                                                                                                              ctcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCGCTCTCTGAACAGCTAAAGGTCAGt                 <  2:210193/149‑1 (MQ=255)
                                                                                                                               tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTa                >  2:714397/1‑149 (MQ=255)
                                                                                                                               tcATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTa                >  2:983013/1‑149 (MQ=255)
                                                                                                                                 aTACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTaaa              >  1:57214/1‑149 (MQ=255)
                                                                                                                                  tACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAg             >  1:518099/1‑149 (MQ=255)
                                                                                                                                        accGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGAGAGTCGACCTGCTGCTGAACCATACCAGCCTGACCCCGGAAGCTCTCTCTGAACAGCTAAAGGGCAGGAAAGACAccc       >  1:185359/3‑148 (MQ=255)
                                                                                                                                         ccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt      >  1:427919/2‑149 (MQ=255)
                                                                                                                                         ccGCCCGCGTCATGCGGGCGGTTATCGAATGATGGGGAGAAAAATATGAAAGCGGCACGCCAGCAAGAGATACTCGACCTGCTGCTGAACAATACCAGCCTGACCACGGAAGCTCACTCTGAACAGCTAAAGGTCAGTAAAGAAAACAt      >  1:299135/2‑149 (MQ=255)
                                                                                                                                         aCGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAATATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAt      <  1:444065/149‑1 (MQ=255)
                                                                                                                                           gCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtc    >  1:499064/1‑149 (MQ=255)
                                                                                                                                             ccGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAACCCATtcgt  >  1:879322/1‑149 (MQ=255)
                                                                                                                                             ccGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCAGACCAGCCTGTCCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgt  >  1:216243/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGT  >  NC_000913/2939186‑2939475

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