Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,157,633 A→T 5.7% intergenic (+59/‑236) ycfH → / → ptsG putative metal‑dependent hydrolase YcfH/glucose‑specific PTS enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,157,6330AT 5.7% 328.3 / 7.4 105intergenic (+59/‑236)ycfH/ptsGputative metal‑dependent hydrolase YcfH/glucose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base A (43/56);  new base T (4/2);  total (47/58)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 108 total aligned reads displayed.

TGTTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGCGCT  >  NC_000913/1157484‑1157781
                                                                                                                                                     |                                                                                                                                                    
tgtTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGaaaa                                                                                                                                                      >  2:794061/1‑150 (MQ=255)
 gtTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAAt                                                                                                                                                     >  1:143996/1‑150 (MQ=255)
      cGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGc                                                                                                                                                >  2:1246370/1‑150 (MQ=255)
       gTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCg                                                                                                                                               >  1:272009/1‑150 (MQ=255)
         tGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGt                                                                                                                                             <  1:404452/150‑1 (MQ=255)
         tGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACAAGTAAGGTTCACCGCCGAAAATTGGGCGGt                                                                                                                                             >  2:253065/1‑150 (MQ=255)
          gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                                                                            >  2:54788/1‑150 (MQ=255)
          gaagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                                                                            <  1:55725/150‑1 (MQ=255)
           aagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGa                                                                                                                                           >  2:912617/1‑150 (MQ=255)
           aagaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGa                                                                                                                                           <  2:441818/150‑1 (MQ=255)
            agaaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGTTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGaa                                                                                                                                          >  2:164624/1‑150 (MQ=255)
              aaCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATa                                                                                                                                        <  1:636470/150‑1 (MQ=255)
               aCTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaactt                                                                                                                                       <  1:1137867/150‑4 (MQ=255)
                 tGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAAcc                                                                                                                                     <  1:1261438/150‑1 (MQ=255)
                 tGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAAcc                                                                                                                                     >  2:1079880/1‑150 (MQ=255)
                  ggCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCa                                                                                                                                    <  1:1397159/150‑1 (MQ=255)
                  ggCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAAAGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCa                                                                                                                                    >  1:1456570/1‑150 (MQ=255)
                     gcAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGt                                                                                                                                 >  2:275058/1‑150 (MQ=255)
                       aGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGttt                                                                                                                               <  1:168149/150‑1 (MQ=255)
                        ggTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTg                                                                                                                              <  1:1087246/150‑1 (MQ=255)
                        ggTAACCACCGATAACTTCGCCCGGCTGTATCACATCGACGCTTCCCGCCGTCAATCCATCCGTTGAAGGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTg                                                                                                                              <  2:1500599/150‑1 (MQ=255)
                           aaCCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGaaa                                                                                                                           <  2:1005252/150‑1 (MQ=255)
                            accaccGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAat                                                                                                                          >  2:1526945/1‑150 (MQ=255)
                                ccGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                                                                      <  2:1185121/150‑1 (MQ=255)
                                  gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                    <  1:683352/150‑1 (MQ=255)
                                  gATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                    <  2:458625/150‑1 (MQ=255)
                                  cctcccgtCGCCCGGCTGTTTCCCCTCGACGCCTCCCCCCCTCAATCCCTCCGTAGACTGAGGTTTTCTCCAGCTCGTAATTCATGGCTAAAACGAGTAAACTTCACCGCCGAACATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                                                                    <  2:304368/143‑1 (MQ=255)
                                    tAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAc                                                                                                                  <  2:1367885/150‑1 (MQ=255)
                                      aCTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACat                                                                                                                >  2:102048/1‑150 (MQ=255)
                                         tCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATg                                                                                                             <  1:84429/150‑1 (MQ=255)
                                              cGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTg                                                                                                        >  1:1157908/1‑150 (MQ=255)
                                              cGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTg                                                                                                        >  1:1157907/1‑150 (MQ=255)
                                                 cTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCa                                                                                                     <  1:938243/150‑1 (MQ=255)
                                                  tGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCaa                                                                                                    >  2:47116/1‑150 (MQ=255)
                                                   gTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAACCAGTAAAATTAAACGCCCAAAATTTGGCCGTTAATAACCAAGTTTTAAATATTGTTACATATGTTTTTTAaaaa                                                                                                  >  2:576445/1‑150 (MQ=255)
                                                    tttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCaaaa                                                                                                  <  1:1142622/150‑1 (MQ=255)
                                                     ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAt                                                                                                 <  1:1282876/150‑1 (MQ=255)
                                                     ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAt                                                                                                 >  2:1376386/1‑150 (MQ=255)
                                                     ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAt                                                                                                 <  2:162602/150‑1 (MQ=255)
                                                      tCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtg                                                                                                <  1:1079880/150‑1 (MQ=255)
                                                      tCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtg                                                                                                <  1:1103022/150‑1 (MQ=255)
                                                        acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                                                              <  1:135150/150‑1 (MQ=255)
                                                         caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                                                                      <  2:1029329/125‑1 (MQ=255)
                                                         caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                                                                      >  1:1029329/1‑125 (MQ=255)
                                                         caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGc                                                                                             >  1:1538046/1‑150 (MQ=255)
                                                             gACGCTTCCCGCCTTCAATCCCTCCGTTGAATGAGTTTTGTTAAAGCTCGTAATTAATGGCTAAACCGCGTAAAGTTCACCGCCGAAAAGTGGGCGGGGAGTTCCCACGTTTGAAATATTGTGCCATATGTTTTGTCAAAATGTGCACCt                                                                                         >  2:1230408/1‑150 (MQ=255)
                                                               cGCTTCCCGCCTTCAATCCTTCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTtc                                                                                       <  2:1239289/150‑1 (MQ=255)
                                                               cGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTtc                                                                                       <  1:1237720/150‑1 (MQ=255)
                                                                 cTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTtctc                                                                                     <  1:54788/150‑1 (MQ=255)
                                                                   tCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAACATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCa                                                                                   <  1:1217999/150‑1 (MQ=255)
                                                                    cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCaaa                                                                                  >  2:821689/1‑150 (MQ=255)
                                                                    cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCaa                                                                                  >  2:1217005/1‑150 (MQ=255)
                                                                         cTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATc                                                                             >  2:1274929/1‑150 (MQ=255)
                                                                            cAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATAATGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGa                                                                          <  1:951218/150‑1 (MQ=255)
                                                                               tccatccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGt                                                                       <  1:990236/150‑1 (MQ=255)
                                                                               tccatccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGt                                                                       <  1:198432/150‑1 (MQ=255)
                                                                               tccatccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGt                                                                       >  1:299322/1‑150 (MQ=255)
                                                                                 catccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTg                                                                     <  1:1217005/150‑1 (MQ=255)
                                                                                 catccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTg                                                                     <  1:1011694/150‑1 (MQ=255)
                                                                                   tccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGGCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGaa                                                                   <  1:1246370/150‑1 (MQ=255)
                                                                                   tccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGaa                                                                   <  1:395168/150‑1 (MQ=255)
                                                                                           aTGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAAACACGTTTGAAATATTGTGACATAAGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATa                                                           >  1:636398/1‑150 (MQ=255)
                                                                                               gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCg                                                       >  1:1525487/1‑150 (MQ=255)
                                                                                                 ttttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTc                                                     <  2:1538046/150‑1 (MQ=255)
                                                                                                 ttttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATACGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTc                                                     <  2:865152/150‑1 (MQ=255)
                                                                                                   ttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaa                                                   <  1:1001527/150‑1 (MQ=255)
                                                                                                       aaaGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaaacaa                                               >  1:1300117/1‑150 (MQ=255)
                                                                                                          gCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAAtt                                            >  1:14710/1‑150 (MQ=255)
                                                                                                          gCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAAtt                                            >  2:1349206/1‑150 (MQ=255)
                                                                                                             cGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGc                                         >  2:331685/1‑150 (MQ=255)
                                                                                                              gTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCa                                        >  1:602421/1‑150 (MQ=255)
                                                                                                                      tGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAatta                                <  2:272009/150‑1 (MQ=255)
                                                                                                                        gCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAattatt                              >  1:357560/1‑150 (MQ=255)
                                                                                                                        gCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAGGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCAGTCAAACAAATTGGCACTGAattatt                              >  1:618491/1‑150 (MQ=255)
                                                                                                                          tAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTAtttt                            <  2:1157907/150‑1 (MQ=255)
                                                                                                                           aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                           >  1:70226/1‑150 (MQ=255)
                                                                                                                           aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                           >  2:672422/1‑150 (MQ=255)
                                                                                                                           aaaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTa                           >  2:1416077/1‑150 (MQ=255)
                                                                                                                            aaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTAc                          <  1:164624/150‑1 (MQ=255)
                                                                                                                            aaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTAc                          <  1:576445/150‑1 (MQ=255)
                                                                                                                              aCGAGTAAAGTTGACCTCCGAAAATTTGGCTGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTGAAAATGTGCAACTTCTCCAATGATCTGAAGTTTAAACGTGATAGCCGTCAAACTAATTGGCACTGAATTATTTTActc                        <  2:699093/150‑1 (MQ=255)
                                                                                                                              aCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTActc                        <  1:912617/150‑1 (MQ=255)
                                                                                                                                gAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCtg                      <  1:1349206/150‑1 (MQ=255)
                                                                                                                                 aGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCtgt                     <  2:196367/150‑1 (MQ=255)
                                                                                                                                  gTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCtgtg                    <  2:143996/150‑1 (MQ=255)
                                                                                                                                     ttattCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaa                 >  1:1164304/4‑150 (MQ=255)
                                                                                                                                     aaGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaa                 >  2:359405/1‑150 (MQ=255)
                                                                                                                                     aaGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaa                 >  1:1324563/1‑150 (MQ=255)
                                                                                                                                        ttCACCGCCCAATTTTCTGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGtaataa              >  2:1090422/1‑150 (MQ=255)
                                                                                                                                         tCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaa             >  1:546227/1‑150 (MQ=255)
                                                                                                                                         tCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaa             >  1:605735/1‑150 (MQ=255)
                                                                                                                                         tCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaa             >  2:522163/1‑150 (MQ=255)
                                                                                                                                          cACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaat            <  1:275058/150‑1 (MQ=255)
                                                                                                                                          cACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaat            <  1:253065/150‑1 (MQ=255)
                                                                                                                                           aCCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTaataaata           <  1:1376386/150‑1 (MQ=255)
                                                                                                                                               ccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAgg       <  2:602421/150‑1 (MQ=255)
                                                                                                                                                cgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAggg      <  1:973320/150‑1 (MQ=255)
                                                                                                                                                cgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAggg      <  1:982051/150‑1 (MQ=255)
                                                                                                                                                   aaaTTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGgcgc   <  1:47116/150‑1 (MQ=255)
                                                                                                                                                    aaTTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGCGCt  >  2:132146/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGTTGCCGTTGAAGAACTGGCGCAGGTAACCACCGATAACTTCGCCCGTCTGTTTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACAAATTGGCACTGAATTATTTTACTCTGTGTAATAAATAAAGGGCGCT  >  NC_000913/1157484‑1157781

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