Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,336,325 C→T 5.2% R107H (CGT→CAT)  yfaA ← DUF2138 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,336,3250CT 5.2% 245.5 / 2.2 76R107H (CGT→CAT) yfaADUF2138 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base T (2/2);  total (36/40)
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.25e-01

CAGTTTGCTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCC  >  NC_000913/2336177‑2336471
                                                                                                                                                    |                                                                                                                                                  
cagTTTGCTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCCCAGTGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGCTCGTGTTCATAGACAATac                                                                                                                                                    >  4:195696/1‑149 (MQ=255)
  gTTTGCTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATacga                                                                                                                                                  >  6:46758/1‑149 (MQ=255)
   tttGCTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATacgac                                                                                                                                                 >  7:85173/1‑149 (MQ=255)
     tGCTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATacgacg                                                                                                                                                >  1:23406/1‑148 (MQ=255)
       cTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGaa                                                                                                                                              <  1:233039/148‑1 (MQ=255)
       cTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAAt                                                                                                                                             >  8:74404/1‑149 (MQ=255)
       cTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAAt                                                                                                                                             >  5:211189/1‑149 (MQ=255)
        ttAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAAtg                                                                                                                                            >  7:9764/1‑149 (MQ=255)
           aaCCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAAtgctg                                                                                                                                         >  3:38256/1‑149 (MQ=255)
           aaCCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAAtgctg                                                                                                                                         >  1:37147/1‑149 (MQ=255)
             ccACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGcc                                                                                                                                       >  1:126572/1‑149 (MQ=255)
              cACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCt                                                                                                                                      >  5:20387/1‑149 (MQ=255)
               aCTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCtt                                                                                                                                     <  4:38256/149‑1 (MQ=255)
                   cgcTGGATCAGCACCATATAATTTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCa                                                                                                                                   <  8:145679/147‑1 (MQ=255)
                    gcTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATg                                                                                                                                >  2:209864/1‑149 (MQ=255)
                        gATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGccc                                                                                                                             <  2:37147/148‑1 (MQ=255)
                        gATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCa                                                                                                                            <  3:231149/149‑1 (MQ=255)
                            aGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGAc                                                                                                                         <  5:46758/148‑1 (MQ=255)
                            aGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGAc                                                                                                                         <  8:9764/148‑1 (MQ=255)
                                       aTGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGACCAGACGATCGGCATGAt                                                                                                             >  3:137020/1‑149 (MQ=255)
                                          tgAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTc                                                                                                           <  2:23406/148‑1 (MQ=255)
                                          tgAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTc                                                                                                           <  7:249548/148‑1 (MQ=255)
                                          tgAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCt                                                                                                          >  4:77533/1‑149 (MQ=255)
                                            aaaGATGGCCTTGTTTAAAGTGCCACAGCGCTGCCTGCGCGGGATGATCTAAGAGTGACGAAAAGAGCATATCTTTCAGCGAGAGATCGTGTTCATAGACAATACGACGAATGCTGCCATCAATGACCAGACGATCGCCATGATACTGa                                                                                                        >  4:153395/1‑149 (MQ=255)
                                               gATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGataa                                                                                                     >  8:96379/1‑149 (MQ=255)
                                                 tGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGataata                                                                                                   <  7:74404/149‑1 (MQ=255)
                                                    ccTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGaa                                                                                                <  5:230010/149‑1 (MQ=255)
                                                     cTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGaa                                                                                                >  2:51065/1‑148 (MQ=255)
                                                     cTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGaa                                                                                                >  1:200033/1‑148 (MQ=255)
                                                           tATCGTGCCACAGCGCTGCCTGCGCGGGCTCATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACTaa                                                                                          >  2:29805/1‑148 (MQ=255)
                                                            aTCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAAt                                                                                        >  7:158962/1‑149 (MQ=255)
                                                                tGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGa                                                                                                                          >  8:4268/1‑111 (MQ=255)
                                                                tGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGa                                                                                                                          <  7:4268/111‑1 (MQ=255)
                                                                 gCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCACAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCg                                                                                   <  8:85173/149‑1 (MQ=255)
                                                                  ccACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCTGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGc                                                                                  <  6:8876/149‑1 (MQ=255)
                                                                   cacaGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATATGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCt                                                                                 >  1:2591/1‑149 (MQ=255)
                                                                    acaGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATATGACGAATGCTGCCTTCAATGCCCAGACGATCg                                                                                                                    <  7:47774/113‑1 (MQ=255)
                                                                    acaGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATATGACGAATGCTGCCTTCAATGCCCAGACGATCg                                                                                                                    >  8:47774/1‑113 (MQ=255)
                                                                    acaGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCt                                                                                 >  6:141278/1‑148 (MQ=255)
                                                                     caGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGa                                                                                                        >  5:183318/1‑124 (MQ=255)
                                                                     caGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGa                                                                                                        <  6:183318/124‑1 (MQ=255)
                                                                     caGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCa                                                                               <  1:209864/149‑1 (MQ=255)
                                                                      aGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCt                                                                                 >  8:199101/1‑146 (MQ=255)
                                                                             cctgcGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCg                                                                       >  2:125926/1‑149 (MQ=255)
                                                                             cctgcGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCg                                                                       >  1:166198/1‑149 (MQ=255)
                                                                             cctgcGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAAAGAAATCTTCGCTCAGAACATc                                                                        <  1:20626/148‑1 (MQ=255)
                                                                                 cgcgGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTg                                                                     <  7:147928/147‑1 (MQ=255)
                                                                                  gcgGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTg                                                                     >  7:145311/1‑146 (MQ=255)
                                                                                  gcgGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGgag                                                                  <  6:211189/149‑1 (MQ=255)
                                                                                      gCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAAc                                                              <  8:158962/149‑1 (MQ=255)
                                                                                          aTCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGa                                                           <  1:218347/148‑1 (MQ=255)
                                                                                               agagTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGt                                                      <  3:77533/148‑1 (MQ=255)
                                                                                                gagTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTga                                                    <  6:20387/149‑1 (MQ=255)
                                                                                                      cGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCtt                                                                                     >  4:27494/1‑110 (MQ=255)
                                                                                                            gagCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTgg                                        <  2:166198/149‑1 (MQ=255)
                                                                                                              gCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGa                                      <  2:200033/149‑1 (MQ=255)
                                                                                                              gCTTATCTATCAGCGTGAGATCGTGTACATAGACAATACGACGTATGCTGCCTTCAATGGCCAGACGATCGGCATGATTCTGATAATAGATAACGAAATCTTCGCTCAGAAAATCGTTGATAAACGGAATGGTGAGGAGATCTTTggg                                       <  3:153395/148‑1 (MQ=255)
                                                                                                               cTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTggg                                       >  7:216195/1‑147 (MQ=255)
                                                                                                                     tttCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATg                                                        >  7:14412/1‑124 (MQ=255)
                                                                                                                     tttCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATg                                                        <  8:14412/124‑1 (MQ=255)
                                                                                                                     tttCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGc                               >  7:157905/1‑149 (MQ=255)
                                                                                                                     tttCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGc                               >  3:153204/1‑149 (MQ=255)
                                                                                                                      ttCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCt                              >  7:205052/1‑149 (MQ=255)
                                                                                                                         aGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGTAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCag                           <  2:126572/149‑1 (MQ=255)
                                                                                                                          gCGTGAGATCGTGTTCATAGACAATATGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCaga                          <  2:2591/149‑1 (MQ=255)
                                                                                                                          gCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCaga                          <  5:141278/149‑1 (MQ=255)
                                                                                                                          gCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCaga                          <  8:216195/149‑1 (MQ=255)
                                                                                                                            gTGAGATCGTGTTCATAGACAATACGACTAATGCAGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCTTGGAGAAACGGAATGGAGAGGACATCTTTg                                         >  8:5545/1‑132 (MQ=255)
                                                                                                                                  tCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCt                  <  5:7945/149‑1 (MQ=255)
                                                                                                                                  tCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCt                  <  8:157905/149‑1 (MQ=255)
                                                                                                                                    gtgtTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGt                <  8:213035/149‑1 (MQ=255)
                                                                                                                                     tgtTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTc               <  8:145311/149‑1 (MQ=255)
                                                                                                                                       ttCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCga             <  1:51065/149‑1 (MQ=255)
                                                                                                                                          aTAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCgaga           >  5:231335/1‑148 (MQ=255)
                                                                                                                                          aTAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCgaga           >  1:144648/1‑148 (MQ=255)
                                                                                                                                             gACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAaga       <  3:94771/149‑1 (MQ=255)
                                                                                                                                              aCATTACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGaaa         >  7:189940/1‑146 (MQ=255)
                                                                                                                                              aCATTACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAaga       <  8:189940/148‑1 (MQ=255)
                                                                                                                                                 aTACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTc   <  7:96379/149‑1 (MQ=255)
                                                                                                                                                  tACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGTTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAACGAGGTcc  >  1:66166/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CAGTTTGCTTAAACCACTGCGCTGGATCAGCACCATATAATGTGAAAGATGGCCTTGTTTATCGTGCCACAGCGCTGCCTGCGCGGGCTGATCTAAGAGTGACGAAAAGAGCTTATCTTTCAGCGTGAGATCGTGTTCATAGACAATACGACGAATGCTGCCTTCAATGCCCAGACGATCGGCATGATTCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCC  >  NC_000913/2336177‑2336471

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: