Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,926 A→G 5.5% I134V (ATT→GTT)  frsA → fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913256,9260AG 5.5% 213.3 / 2.9 73I134V (ATT→GTT) frsAfermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base G (2/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

TATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATCACCGG  >  NC_000913/256780‑257074
                                                                                                                                                  |                                                                                                                                                    
tataCGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACAtt                                                                                                                                                    >  1:369460/1‑149 (MQ=255)
 atacgatacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTg                                                                                                                                                   <  1:24734/149‑1 (MQ=255)
   acgataccGTGATTGGAGAGCGGGGCGGCAAAGGGGTTCATGAGTGGGCATCGCAGGCGAGTGTGTGGCAACAATTAGCGTGTGCGGAAGACGATCCGCTACTCATTGGTCGTCACTGGCTGCATGCGTCTACGTTGTACAACATTgaa                                                                                                                                                 <  2:918337/149‑3 (MQ=255)
     gatacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTtccgc                                                                                                                                               >  1:655037/1‑149 (MQ=255)
       aacgGTGATTGGCTACCGTGGCGGCAACTGGAGTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcat                                                                                                                                             <  1:637493/148‑3 (MQ=255)
          ggTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATc                                                                                                                                          <  1:752720/149‑1 (MQ=255)
          ggTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGATGTACAACATTGCCGCCTAta                                                                                                                                          >  2:424222/1‑148 (MQ=255)
           gttaTTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCAACCAGGCGATGGTGTGGCAACAAAAAGCATGTTAGGAAGACGATCCGCAACTCAGTGGTCTGCACTTACTTCATGCGGCTACGTTGTACAACATTGCCGCCTATca                                                                                                                                         >  2:764061/4‑148 (MQ=255)
                   cTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGa                                                                                                                                 >  1:529288/1‑149 (MQ=255)
                   cTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGa                                                                                                                                 >  1:406476/1‑149 (MQ=255)
                             ggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATACGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcatg                                                                                                                                            >  1:1020333/1‑125 (MQ=255)
                             ggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATACGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcatg                                                                                                                                            <  2:1020333/128‑4 (MQ=255)
                                aaCTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCt                                                                                                                    >  2:270283/1‑149 (MQ=255)
                                aaCTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCt                                                                                                                    <  1:422448/149‑1 (MQ=255)
                                 aCTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTg                                                                                                                   <  1:33216/149‑1 (MQ=255)
                                   tGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAACGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCGGGc                                                                                                                 >  2:320280/1‑149 (MQ=255)
                                   tGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGc                                                                                                                 >  2:83286/1‑149 (MQ=255)
                                        ttATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGc                                                                                                            >  2:946248/1‑149 (MQ=255)
                                           tGAGTGGGCCACACAGGCGAGGGTGTGTCAACAAATCGCCTGTGCGTAAGACGATCCGCAACGCATTGGTCGTCACTGGCTGCATGCGGCTACGTTGGACAACATTGCCGCCTATCCTCATCTGAAAGTAGATGACCTGGCCGAGCAAg                                                                                                         <  1:919874/149‑1 (MQ=255)
                                             aGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAgcg                                                                                                       >  1:269356/1‑149 (MQ=255)
                                                   ccACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCt                                                                                                 >  1:150510/1‑149 (MQ=255)
                                                    cACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCtt                                                                                                >  2:843929/1‑149 (MQ=255)
                                                     aCCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCttt                                                                                               <  1:843929/149‑1 (MQ=255)
                                                      cccAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTg                                                                                              <  1:300165/149‑1 (MQ=255)
                                                         aGGCGATGGTGTGGCAACAAAACGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTcc                                                                                           >  1:966222/1‑148 (MQ=255)
                                                         aGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCa                                                                                           <  2:147041/149‑1 (MQ=255)
                                                         aGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGATTTGTCa                                                                                           >  2:113051/1‑149 (MQ=255)
                                                          ggCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCaa                                                                                          <  1:592007/149‑1 (MQ=255)
                                                           gCGATGGTGAGGCAACAAAAAGCCTGTGCGGAAGACGATCCCCACCTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCGGCATATACTCATCTGAAAGGAGATGAGCTGGCCGAGCAGTAGCAGGCTTTGTCaaa                                                                                         >  1:848394/1‑149 (MQ=255)
                                                                  ggtgGCAACAAAAAGCCGGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGcc                                                                                  <  2:949179/148‑1 (MQ=255)
                                                                   gtgGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCCGGCTTTGTCAAACCGCGCCt                                                                                 >  1:581032/1‑149 (MQ=255)
                                                                    tgGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGATCGTCACTGGCTGCATGCGGCTACGTGGTACAACATTGCCGCCTATCCTCATCTGATAGGCGAAGCCCGGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCtt                                                                                >  2:748141/1‑148 (MQ=255)
                                                                      gCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATg                                                                              >  2:877100/1‑149 (MQ=255)
                                                                         acaaAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAag                                                                           >  2:378811/1‑149 (MQ=255)
                                                                         acaaAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAag                                                                           <  1:115850/149‑1 (MQ=255)
                                                                           aaaaaGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGCCCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAagag                                                                         >  1:706979/1‑149 (MQ=255)
                                                                            aaaaGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGCTGACCTGGCCGAGCAAGCGCAGGCTTTGTCAACCCGCGCCTATGAAGAgg                                                                        >  1:671681/1‑149 (MQ=255)
                                                                            aaaaGCCGGTGCGGAAGACGATCCGAAACTCAGTGGTCTTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAgg                                                                        <  1:597850/149‑1 (MQ=255)
                                                                             aaaGCCGGTGCGGCAGAGGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTAGGTTGTACAACATTGCCGCATATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGc                                                                       <  1:439445/149‑1 (MQ=255)
                                                                              aaGCATGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGATGTACAACAATGCCGCATATCCTCATCTGAAAGTAGATCACCTGGCCGATCAAGAGAAGGCTTTTTCAACACGAGCATATGAAGAGGcc                                                                      >  1:254656/1‑149 (MQ=255)
                                                                                 ccTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCt                                                                   <  1:445958/149‑1 (MQ=255)
                                                                                     tgCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTCCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGc                                                               <  2:39503/149‑1 (MQ=255)
                                                                                         gAAGACGATCCGCAACTCAGTGGTCGTCACTGGCGGCATGCGGCTACGTTGTACAACATTGGCGCCTTTCCTCATCTGAAAGGAGATGCCCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCt                                                           >  2:190989/1‑149 (MQ=255)
                                                                                          aaGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTa                                                          >  1:724810/1‑149 (MQ=255)
                                                                                           aGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTAc                                                         <  2:838852/149‑1 (MQ=255)
                                                                                              cGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCgg                                                      >  1:92821/1‑149 (MQ=255)
                                                                                               gATCCGCAACTCAGGTGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCggg                                                     <  1:703517/149‑1 (MQ=255)
                                                                                               gATCCGCAACTAAGTTGTCGTCACTGGCTGCATGCGGATACGTTGTACAACATTGCAGCCAATCCTCATCTGAAAGTAGATGAACTGGCCGATCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCggg                                                     >  2:1014321/1‑149 (MQ=255)
                                                                                                  ccGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCAc                                                  <  1:877100/149‑1 (MQ=255)
                                                                                                   cGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGCCCTGGCCGAGCAAGCGCAGGCTTTGTCACACCGCGCCTATGAAGAGGCCGCTCAGCGTCTCCCGGGCACg                                                 >  1:805683/1‑149 (MQ=255)
                                                                                                   cGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAACCCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACg                                                 >  2:782253/1‑149 (MQ=255)
                                                                                                   cGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACg                                                 >  2:804313/1‑149 (MQ=255)
                                                                                                    gCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGa                                                >  1:302779/1‑149 (MQ=255)
                                                                                                      aaCTCAGTGGTCTTAACTTGATGCATTCGGCTACGTTTTACAACATTGCCGCATATCCTAATATGAAATTAGATTACCTGGACGAGCAAGCGCAGGATTTGTCAAACCGCGCCTCTGAAGAGGCCGCTCAGCGTCTACCGGGCTCGATg                                              <  1:764061/149‑1 (MQ=255)
                                                                                                      aaCTCAGAGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCAGCTGAAAGGAGATGACCTGGCCTAGCAGGCGCAGGCTTTGTCAAAGCGCGCCTATGAAGAGGCCGCTCAGCGTCTTCCGGGCACGATg                                              >  2:603939/1‑149 (MQ=255)
                                                                                                       aCTCAGTGGGCGTCACTGCCTGCATGCGGCTACGTGGTACACCATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGTCCGCTCAGCGTCTACCGGGCACGATGc                                             <  2:748266/149‑1 (MQ=255)
                                                                                                        atCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTTCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGATCAAGCGCAGGATTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCg                                            <  2:559575/148‑1 (MQ=255)
                                                                                                            gTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAg                                        <  1:376707/149‑1 (MQ=255)
                                                                                                             tGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGa                                       >  2:170748/1‑149 (MQ=255)
                                                                                                                 cggcACTGGCGGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGa                                   <  2:1015719/146‑1 (MQ=255)
                                                                                                                  gtcACTGGCTGCATGCTGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGATCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAg                                  <  2:274081/149‑1 (MQ=255)
                                                                                                                  gtcACTGGCTGCATGCGGTTACGTTGTACAACATTGCCGCCTATCCTCAACTGAAAGGAGATGTCCTGGCAGAGCAAGAGCAGGCTTTGTCTAACCGCGACTATGATGATGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAg                                  <  2:278124/149‑1 (MQ=255)
                                                                                                                    cACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGtt                                >  1:535852/1‑149 (MQ=255)
                                                                                                                         gCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCCGCGTCTACCGGGCACGCGGCGGCAGATGGAGTTTCCCg                           >  1:74225/1‑149 (MQ=255)
                                                                                                                            gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAc                        >  2:429679/1‑149 (MQ=255)
                                                                                                                            ataggCGGCAATGTTGTACAACGTAGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAc                        >  1:346366/5‑149 (MQ=255)
                                                                                                                              aTGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAccc                      >  2:1020241/1‑149 (MQ=255)
                                                                                                                              aTGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGACAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGAAGCGGCAGATGGAGTTTACCGTAccc                      >  2:3604/1‑149 (MQ=255)
                                                                                                                                gCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTGTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCcgg                    <  1:314936/149‑1 (MQ=255)
                                                                                                                                 cGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGGCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCcggc                   <  1:6464/149‑1 (MQ=255)
                                                                                                                                  ggCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCcggcg                  <  1:386088/149‑1 (MQ=255)
                                                                                                                                    cTACGTTGTACAACATTCCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGGCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGt                <  1:192395/149‑1 (MQ=255)
                                                                                                                                        gTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTgcgc            >  2:285635/1‑149 (MQ=255)
                                                                                                                                             acaacaTTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATc       >  2:544870/1‑149 (MQ=255)
                                                                                                                                                acaTTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATCCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATCAcc    >  1:945381/1‑149 (MQ=255)
                                                                                                                                                 caTTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATCACCg   <  2:529288/149‑1 (MQ=255)
                                                                                                                                                 caTTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATCACCg   <  2:406476/149‑1 (MQ=255)
                                                                                                                                                  aTTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATCACCgg  >  2:16189/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGCGCCCATCACCGG  >  NC_000913/256780‑257074

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: