Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,262,665 G→A S277S (AGC→AGT 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,262,6650GA98.6% 258.0 / ‑4.0 69S277S (AGC→AGTfruBfused fructose‑specific PTS enzymes: IIA component/HPr component
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/30);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCGCTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTGCCGTTTCG  >  NC_000913/2262519‑2262812
                                                                                                                                                  |                                                                                                                                                   
ccggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTgg                                                                                                                                                 <  2:937025/151‑1 (MQ=255)
 cggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTgg                                                                                                                                                 >  1:386827/1‑150 (MQ=255)
 cggtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTGGTCGGCGCATCATCACTGGTc                                                                                                                                               >  1:362585/1‑151 (MQ=255)
   gtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTc                                                                                                                                               >  2:667895/1‑150 (MQ=255)
   gtACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTc                                                                                                                                               >  1:447735/1‑150 (MQ=255)
     aCCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTcagc                                                                                                                                            >  2:362412/1‑151 (MQ=255)
         tCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAgcg                                                                                                                                        >  2:961385/1‑151 (MQ=255)
           aaGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGcc                                                                                                                                      >  2:18436/1‑151 (MQ=255)
           aaGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGcc                                                                                                                                      >  1:612486/1‑151 (MQ=255)
            aGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCa                                                                                                                                     <  1:271191/151‑1 (MQ=255)
                    cACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGtt                                                                                                                             >  1:122501/1‑151 (MQ=255)
                     aCGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTg                                                                                                                            >  1:96146/1‑151 (MQ=255)
                      cGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGc                                                                                                                           <  2:511497/151‑1 (MQ=255)
                      cGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGc                                                                                                                           >  1:810549/1‑151 (MQ=255)
                      cGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGc                                                                                                                           >  1:853050/1‑151 (MQ=255)
                        gTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCg                                                                                                                         >  2:871057/1‑151 (MQ=255)
                        gTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCg                                                                                                                         >  2:571374/1‑151 (MQ=255)
                         tAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGc                                                                                                                        >  1:934770/1‑151 (MQ=255)
                         tAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGc                                                                                                                        >  1:437400/1‑151 (MQ=255)
                          aaTATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCa                                                                                                                       <  2:712351/151‑1 (MQ=255)
                           atatCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCAt                                                                                                                      >  1:800529/1‑151 (MQ=255)
                            tatCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATc                                                                                                                     >  1:637310/1‑151 (MQ=255)
                               cACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGc                                                                                                                  >  2:284150/1‑151 (MQ=255)
                                aCTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCt                                                                                                                 <  1:659862/151‑1 (MQ=255)
                                      aaaTTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGc                                                                                                           >  1:926548/1‑151 (MQ=255)
                                           gTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGc                                                                                                       <  2:374812/150‑1 (MQ=255)
                                           gTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGc                                                                                                       <  2:688729/150‑1 (MQ=255)
                                           gTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCg                                                                                                      <  1:391739/151‑1 (MQ=255)
                                              tAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGt                                                                                                    <  1:2532/150‑1 (MQ=255)
                                                 tGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGc                                                                                                >  2:606308/1‑151 (MQ=255)
                                                  ggTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCt                                                                                               <  1:287469/151‑1 (MQ=255)
                                                     aTTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCttt                                                                                             <  1:264490/150‑1 (MQ=255)
                                                     aTTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTa                                                                                            <  2:268398/151‑1 (MQ=255)
                                                      ttGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTAt                                                                                           >  2:653063/1‑151 (MQ=255)
                                                       tGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTAt                                                                                           <  1:56614/150‑1 (MQ=255)
                                                        gACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTAtt                                                                                          <  2:224683/150‑1 (MQ=255)
                                                        gACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTg                                                                                         >  1:973852/1‑151 (MQ=255)
                                                            aGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCg                                                                                      >  2:770792/1‑150 (MQ=255)
                                                              cATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGc                                                                                   >  2:477526/1‑151 (MQ=255)
                                                                tGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCa                                                                                  >  2:937750/1‑150 (MQ=255)
                                                                 ggTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCaa                                                                                 <  2:262957/150‑1 (MQ=255)
                                                                 ggTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAAt                                                                                >  1:461512/1‑151 (MQ=255)
                                                                  gTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATa                                                                               <  2:495281/151‑1 (MQ=255)
                                                                    aCCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATaa                                                                              >  2:317164/1‑150 (MQ=255)
                                                                        ggACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCa                                                                          <  1:931102/150‑1 (MQ=255)
                                                                         gACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGc                                                                        <  2:509340/151‑1 (MQ=255)
                                                                          aCGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTTAGCTTTATTGTCGAGCAATAAATCAGCg                                                                       >  1:1044815/1‑151 (MQ=255)
                                                                          aCGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCg                                                                       <  1:673154/151‑1 (MQ=255)
                                                                            gAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCgag                                                                     <  1:821130/151‑1 (MQ=255)
                                                                                     ggCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAg                                                            <  1:1005534/151‑1 (MQ=255)
                                                                                       ccGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGaa                                                          <  2:657924/151‑1 (MQ=255)
                                                                                           gtTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACcg                                                       <  2:926703/150‑1 (MQ=255)
                                                                                           gtTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACcg                                                       <  1:961225/150‑1 (MQ=255)
                                                                                            ttCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACcgc                                                      <  1:453412/150‑1 (MQ=255)
                                                                                             tCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACcgcg                                                     >  1:627289/1‑150 (MQ=255)
                                                                                                    gcACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTg                                             <  1:569171/151‑1 (MQ=255)
                                                                                                           aaaCTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCGCTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGTGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGtc                                       >  1:409027/1‑150 (MQ=255)
                                                                                                             aCTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGtcatt                                    >  2:201569/1‑151 (MQ=255)
                                                                                                              cTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGtcatt                                    >  2:248854/1‑150 (MQ=255)
                                                                                                                        aaCACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCcaca                          <  1:18434/150‑1 (MQ=255)
                                                                                                                        aaCACGTCGTCGGTCGGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCacac                         >  1:983873/1‑151 (MQ=255)
                                                                                                                                gtcggtcgGCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGc                 >  2:916640/1‑151 (MQ=255)
                                                                                                                                  cggtcggCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCAt               >  2:984509/1‑151 (MQ=255)
                                                                                                                                       ggCGCATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTg          <  2:92195/151‑1 (MQ=255)
                                                                                                                                        gcgcATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTGc         >  1:658300/1‑151 (MQ=255)
                                                                                                                                         cgcATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTGcc        >  1:195604/1‑151 (MQ=255)
                                                                                                                                          gcATCATCACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTGCCg       <  1:85907/151‑1 (MQ=255)
                                                                                                                                            atcatcACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTGCCGtt     >  2:795170/1‑151 (MQ=255)
                                                                                                                                               atcACTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATAGTCATTCATCGCCACACTCACCAGCATGGCTGCCGTTTCg  <  2:612589/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCGGTACCATCAAGGTTTGTCACGGTAATATCACTGTTAAATTGTTTAATGGTATTGACCAGCATGGTACCTGGACGAGCATGCAGGCCGTGTTCATTGCGCACCACAAACTCCGCGCTTAACACGTCGTCGGTCGGCGCATCATCGCTGGTCAGCAGCGCCAGCAACGTTGCCGCATCCGCTTTCAGCAAGCGGTCAGCTTTATTGTCGAGCAATAAATCAGCGAGACGCTTAAGAACCGCGATGGGCTGATCGTCATTCATCGCCACACTCACCAGCATGGCTGCCGTTTCG  >  NC_000913/2262519‑2262812

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