Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,263,236 C→G 6.4% W87S (TGG→TCG) ‡ fruB ← fused fructose‑specific PTS enzymes: IIA component/HPr component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,263,2360CG 6.4% 282.9 / 8.1 92W87S (TGG→TCG) ‡fruBfused fructose‑specific PTS enzymes: IIA component/HPr component
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base G (2/4);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCT  >  NC_000913/2263087‑2263385
                                                                                                                                                     |                                                                                                                                                     
cgcgAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGGAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCAcccc                                                                                                                                                       <  1:829791/150‑1 (MQ=255)
        tcttctGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCGGGCAGCAATACCGATTGCCACGTACGCTACCAGACCGTCACCCCAGGTGACg                                                                                                                                               >  1:251392/1‑150 (MQ=255)
            ctgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtca                                                                                                                                           <  2:1178770/150‑5 (MQ=255)
            ctgctgTTGTTGCTGACGTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCtg                                                                                                                                           >  1:1071704/1‑149 (MQ=255)
             tgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtcag                                                                                                                                          <  1:1375554/150‑6 (MQ=255)
             tgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtcag                                                                                                                                          <  1:1488018/150‑6 (MQ=255)
             tgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTc                                                                                                                                          >  1:1431446/1‑150 (MQ=255)
             tgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTc                                                                                                                                          >  2:780627/1‑150 (MQ=255)
             tgctgTTGTTGCTGACTTCAGTTGGTCAGCAACGGAATCATCGCTCAGTACGGGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtcag                                                                                                                                          <  2:414119/150‑6 (MQ=255)
                   tgttgCTGGCTTCTGTTGTTCAGCAACGGATTCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCGCAGATGTTCATCCGAGCTGGCAGCAATAACGATTGACACGTACGCTAACTGACCGTCACACCAGGTGACGCCTTCGGGGaa                                                                                                                                    >  2:921405/1‑150 (MQ=255)
                   tgttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGCCCGTCCCCCCAGGTGACGCCTTCCGGGaa                                                                                                                                    >  2:408167/1‑150 (MQ=255)
                    gttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGGCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAAc                                                                                                                                   <  2:1529583/150‑1 (MQ=255)
                     ttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACt                                                                                                                                  <  1:472004/150‑1 (MQ=255)
                      tgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTg                                                                                                                                 <  1:322295/150‑1 (MQ=255)
                      ggCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGGCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATGGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTg                                                                                                                                 <  2:1538605/149‑1 (MQ=255)
                        cTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGCCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTCACGCCTTCCGGGAACTGaa                                                                                                                               >  2:318541/1‑150 (MQ=255)
                          gACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAAt                                                                                                                             <  1:399906/150‑1 (MQ=255)
                          gACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAAt                                                                                                                             <  1:406200/150‑1 (MQ=255)
                                     ttCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGCCGCCTTCCGGGAACAGAAATACCTGAACGcc                                                                                                                  >  2:784838/1‑150 (MQ=255)
                                       cAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCgg                                                                                                                >  2:8370/1‑150 (MQ=255)
                                       cAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCgg                                                                                                                >  1:502181/1‑150 (MQ=255)
                                        aGCAACGGAATCATCGCGCAGTACGTGGGTCAGGTGGCGTAGCTGGCCAAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTTACGCCTTCCGGGAACTGAAATACCTGAACGCCGGt                                                                                                               <  1:1556990/150‑1 (MQ=255)
                                          cAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGttt                                                                                                             <  1:460445/150‑1 (MQ=255)
                                           aaCGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGtttt                                                                                                            >  2:609748/1‑150 (MQ=255)
                                            aCGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTc                                                                                                           <  2:566718/150‑1 (MQ=255)
                                            aCGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTc                                                                                                           >  2:348627/1‑150 (MQ=255)
                                             cGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                          <  1:1631804/150‑1 (MQ=255)
                                             cGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                          >  1:491440/1‑150 (MQ=255)
                                             cGGAATCATCGCTCAGTACGTGGGGCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                          <  1:374685/150‑1 (MQ=255)
                                                  tcatcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGGCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTGTTCAGCAcc                                                                                                     <  1:97917/150‑1 (MQ=255)
                                                   catcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGACCGCCGGTTTTCAGccccg                                                                                                    >  1:133741/1‑149 (MQ=255)
                                                   catcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCt                                                                                                    <  1:479729/150‑1 (MQ=255)
                                                       gCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTCAACGCAGGATATCAGCACCTGATc                                                                                                >  1:355685/1‑150 (MQ=255)
                                                        cTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACGCCAGTTGACGCCTTCCGGGAACTGATATACCTGAACCCCGCTTTTCAGCACCTGTTcg                                                                                               >  1:77040/1‑150 (MQ=255)
                                                           aGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                            <  2:1263610/150‑1 (MQ=255)
                                                                    gTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCCCCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCggga                                                                                   >  1:791815/1‑148 (MQ=255)
                                                                        gCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCCCCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGc                                                                               >  2:1469520/1‑150 (MQ=255)
                                                                        gCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGc                                                                               <  2:879876/150‑1 (MQ=255)
                                                                        gCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGc                                                                               >  2:1385612/1‑150 (MQ=255)
                                                                           ggCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGCTCGCGGGTGTCGGTAGTGCCgt                                                                            >  1:894208/1‑150 (MQ=255)
                                                                            gCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGAACGCTACCTGACCGTCACCCCAGTTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCgtg                                                                           <  2:1336658/150‑1 (MQ=255)
                                                                              gTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGtgtg                                                                         <  1:1352256/150‑1 (MQ=255)
                                                                               tAGCAGGCCCAGATTTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTgg                                                                        <  1:919195/150‑1 (MQ=255)
                                                                                   aGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGAGACGCCTTCCGGGAACTGAAAAACCTGAACGCCGGGTTTCAGCACCTCATCGCGGGTGTCGGGAGTGCCGTCTGGAATa                                                                    >  1:977509/1‑150 (MQ=255)
                                                                                    ggCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAg                                                                   <  1:8370/150‑1 (MQ=255)
                                                                                    ggCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAg                                                                   >  1:83931/1‑150 (MQ=255)
                                                                                        cAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAAt                                                               >  2:1402797/1‑150 (MQ=255)
                                                                                           aTGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATAcc                                                            <  2:668951/150‑1 (MQ=255)
                                                                                           aTGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCGACCTGACCGTCACCCCAGGGGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATAcc                                                            <  2:113528/150‑1 (MQ=255)
                                                                                            tGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCa                                                           >  1:166436/1‑150 (MQ=255)
                                                                                            tGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCa                                                           >  1:1232703/1‑150 (MQ=255)
                                                                                             gTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCGTTCCGGGATCTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCAt                                                          <  1:1204907/150‑1 (MQ=255)
                                                                                             gTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGATTAGCAATACCAtt                                                         <  1:346709/150‑1 (MQ=255)
                                                                                                    cGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAg                                                   <  1:1494101/150‑1 (MQ=255)
                                                                                                        cTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCCCCCCAGGTGCCGCCTTCCGGGAACTGAAATACCTGCACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCCTTGCCGAGGAAc                                               >  1:1725774/1‑150 (MQ=255)
                                                                                                        cTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAAc                                               >  1:1611141/1‑150 (MQ=255)
                                                                                                          ggCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCAGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCAGTGCCGAGGAACgg                                             >  2:1096293/1‑149 (MQ=255)
                                                                                                            cagcaATACCGATTGGCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCTGGAACTGAAAGCGCTGACCGCCGGTGTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGAGGAATAGCAATACCAATGCCGAGGAACGTTg                                           >  1:296114/1‑150 (MQ=255)
                                                                                                            cagcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTg                                           >  2:1117039/1‑150 (MQ=255)
                                                                                                            cagcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGACGAGGAACGTTg                                           >  2:241082/1‑150 (MQ=255)
                                                                                                               caATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGATCTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGGGGCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAgg                                        <  2:268552/150‑1 (MQ=255)
                                                                                                               caATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAgg                                        <  1:619012/150‑1 (MQ=255)
                                                                                                                aaTACCGATTGCCACGTACGCTACCTGACCGTCAGCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGt                                       <  2:149971/150‑1 (MQ=255)
                                                                                                                aaTACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCCTTGCCGAGGAACGTTGAGGt                                       >  1:1496939/1‑150 (MQ=255)
                                                                                                                   aCCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTg                                    <  2:1556360/150‑1 (MQ=255)
                                                                                                                    ccGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGc                                   <  2:1700383/150‑1 (MQ=255)
                                                                                                                      gATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGCACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTg                                 >  1:342963/1‑150 (MQ=255)
                                                                                                                      gATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTg                                 >  1:619489/1‑150 (MQ=255)
                                                                                                                        ttGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGtt                               <  2:83931/150‑1 (MQ=255)
                                                                                                                        ttGCCACGTACGCTACCTGACCGGCAGCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGATTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGtt                               <  1:318541/150‑1 (MQ=255)
                                                                                                                          gCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTcg                             >  2:149640/1‑150 (MQ=255)
                                                                                                                            cACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCAGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAAAAGCAATACCAATGCCGAGGAACGTTGAGGTTTGCTGTTcgcg                           >  2:722911/1‑150 (MQ=255)
                                                                                                                              cGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCGGTGTGGATTAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTcgcgcg                         >  1:1315096/1‑150 (MQ=255)
                                                                                                                               gTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCgcgcgc                        >  2:173437/1‑150 (MQ=255)
                                                                                                                                tACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGTGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGcc                       <  1:57739/150‑1 (MQ=255)
                                                                                                                                   gCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATATCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGc                    <  1:921405/150‑1 (MQ=255)
                                                                                                                                     tACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCAt                  >  2:611940/1‑150 (MQ=255)
                                                                                                                                      aCCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATg                 <  1:424036/150‑1 (MQ=255)
                                                                                                                                       ccTGACCGTCACCCCAGTTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGCTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGc                <  1:84335/150‑1 (MQ=255)
                                                                                                                                       ccTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                          <  1:460383/60‑1 (MQ=255)
                                                                                                                                       ccTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                          >  2:460383/1‑60 (MQ=255)
                                                                                                                                        cTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGcc               <  1:1402797/150‑1 (MQ=255)
                                                                                                                                        cTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGcc               >  2:900900/1‑150 (MQ=255)
                                                                                                                                        cTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGcc               <  1:1666084/150‑1 (MQ=255)
                                                                                                                                           ggcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCAtt            >  1:228445/3‑150 (MQ=255)
                                                                                                                                            ccGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTg           <  2:1413285/150‑1 (MQ=255)
                                                                                                                                            ccGGCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTg           <  2:1188324/150‑1 (MQ=255)
                                                                                                                                             cGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGATGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGa          >  1:833358/1‑150 (MQ=255)
                                                                                                                                             cGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGa          >  2:1587547/1‑150 (MQ=255)
                                                                                                                                               tCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACg        <  1:490314/150‑1 (MQ=255)
                                                                                                                                               gcACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGCACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACg        <  1:822345/149‑1 (MQ=255)
                                                                                                                                                cACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGATCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGt       <  2:1181258/150‑1 (MQ=255)
                                                                                                                                                cACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGt       >  2:997861/1‑150 (MQ=255)
                                                                                                                                                  ccccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAg     >  1:525461/1‑150 (MQ=255)
                                                                                                                                                  ccccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAg     <  1:587007/150‑1 (MQ=255)
                                                                                                                                                  ccccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAg     <  2:1232703/150‑1 (MQ=255)
                                                                                                                                                   cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGc    >  2:1519251/1‑150 (MQ=255)
                                                                                                                                                     cAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCt  >  1:72848/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCCT  >  NC_000913/2263087‑2263385

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: