Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,915 C→T 5.1% T130M (ACG→ATG)  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,9150CT 5.1% 219.8 / 2.4 80T130M (ACG→ATG) frsAfermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base T (2/2);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

ACCGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGT  >  NC_000913/256767‑257060
                                                                                                                                                    |                                                                                                                                                 
aCCGTCGATGATTTATTCGATACGGTGATTGGCTTCCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTAc                                                                                                                                                   <  2:103993/149‑1 (MQ=255)
 ccGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACg                                                                                                                                                  >  1:641916/1‑149 (MQ=255)
 ccGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACg                                                                                                                                                  >  2:6464/1‑149 (MQ=255)
   gacgacGATTTATACGATACGGTGATTGGCTCCCGTGGCGGCAACTGGATTTATGCGTGGGCCACCCAGGCGATGGGGTGGCAACAAACAGCCGGTGCGGAAGACGAGCCGCAACTCAGTGGGCGTCACTGGCTGCATGCGGCTACGtt                                                                                                                                                <  2:745795/149‑1 (MQ=255)
      gatgatTTATACGATACGGTGATTGGCTCCCGTGGCGGCAACTGGAGTTATGAGGGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGGCGTCACTGGCTGCATGCGGCTACGTTGTa                                                                                                                                             <  1:454859/149‑1 (MQ=255)
        tgatTTATACGATACGGTGATTGGCTACCTTGGCGGCAACTGGATTTATGAGTGGGCCAACCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTAAGTGGTCGTCACTGGATGCATGCGGCTACGTTGTaca                                                                                                                                           >  2:772445/1‑149 (MQ=255)
         gatTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTacaa                                                                                                                                          <  1:766350/149‑1 (MQ=255)
         gatTTATACGATACGGGGATTGGCTACCGGGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTacaa                                                                                                                                          <  1:527627/149‑1 (MQ=255)
          atTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGTATTTATGAGTGGGCCACCCAGTCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGAGCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTacatc                                                                                                                                         <  1:345986/149‑3 (MQ=255)
             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:406476/1‑149 (MQ=255)
                                cTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGa                                                                                                                   >  1:529288/1‑149 (MQ=255)
                                          ggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATACGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcatg                                                                                                                              >  1:1020333/1‑125 (MQ=255)
                                          ggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATACGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcatg                                                                                                                              <  2:1020333/128‑4 (MQ=255)
                                             aaCTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCt                                                                                                      <  1:422448/149‑1 (MQ=255)
                                             aaCTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCt                                                                                                      >  2:270283/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                 
ACCGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGT  >  NC_000913/256767‑257060

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: