Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,309,028 G→A 5.4% T105M (ACG→ATG)  alkB ← oxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,309,0280GA 5.4% 240.4 / 2.5 74T105M (ACG→ATG) alkBoxidative demethylase of N1‑methyladenine or N3‑methylcytosine DNA lesions
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AAAAATCGCGGGTAAGCCCAGAGAAACAGAAACAATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGA  >  NC_000913/2308880‑2309170
                                                                                                                                                    |                                                                                                                                              
aaaaaTCGCGGGTAAGCCCAGAGAAACAGAAACAATTGGCGCGCGCAGATCCGGGTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCccgccgccg                                                                                                                                                <  5:107795/149‑1 (MQ=255)
        cgGGTAAGCCCAGAGAAACAGAAACACTTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCa                                                                                                                                        <  1:418213/149‑1 (MQ=255)
           gTAAGCCCAGAGAAACAGAAACAATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGGTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGt                                                                                                                                     >  1:377506/1‑149 (MQ=255)
            tAAGCCCAGAGAAACAGAAACAATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGtt                                                                                                                                    >  4:298951/1‑149 (MQ=255)
                   agagAAACAGAAACAATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATa                                                                                                                              >  5:160689/1‑148 (MQ=255)
                      gaaacagaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAAtt                                                                                                                          <  5:304031/149‑1 (MQ=255)
                       aaacagaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTa                                                                                                                         >  7:111118/1‑149 (MQ=255)
                         acagaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTAt                                                                                                                        >  2:410702/1‑148 (MQ=255)
                          cagaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATg                                                                                                                       >  1:392468/1‑148 (MQ=255)
                           agaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGa                                                                                                                      >  3:336250/1‑148 (MQ=255)
                            gaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGaa                                                                                                                     >  4:179654/1‑148 (MQ=255)
                            gaaacaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGaa                                                                                                                     >  4:408585/1‑148 (MQ=255)
                               acaATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAATTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAAct                                                                                                                 >  7:436734/1‑149 (MQ=255)
                                   ttGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCtgt                                                                                                             <  6:160689/149‑1 (MQ=255)
                                         cgcgcAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATg                                                                                                       <  1:410702/149‑1 (MQ=255)
                                          gcgcAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATgg                                                                                                      >  7:149247/1‑149 (MQ=255)
                                            gcAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAGCTCTGTGGCATGGCg                                                                                                    >  6:304533/1‑149 (MQ=255)
                                              aGATCCGGTTCGTCTTTATCCTGATGCAGCGACAATTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCggg                                                                                                  <  8:436734/149‑1 (MQ=255)
                                                     gTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGt                                                                                           <  1:180796/149‑1 (MQ=255)
                                                     gTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGt                                                                                           <  7:327219/149‑1 (MQ=255)
                                                         gTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAt                                                                                       <  8:111118/149‑1 (MQ=255)
                                                           cTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGtt                                                                                                                                    <  5:103538/102‑1 (MQ=255)
                                                           cTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGtt                                                                                                                                    >  6:103538/1‑102 (MQ=255)
                                                           cTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAttt                                                                                     >  3:92645/1‑149 (MQ=255)
                                                              tATCCTGATGCAGCGCAATTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAACCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAtttgt                                                                                   <  4:298846/148‑1 (MQ=255)
                                                              tATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTg                                                                                                                                   <  6:241351/100‑1 (MQ=255)
                                                              tATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAttt                                                                                     <  2:392468/146‑1 (MQ=255)
                                                              tATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGAAAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTg                                                                                                                                   >  5:241351/1‑100 (MQ=255)
                                                               aTCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAAt                                                                                                                           <  7:123173/107‑1 (MQ=255)
                                                               aTCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAAt                                                                                                                           >  8:123173/1‑107 (MQ=255)
                                                                tCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAtttgtttg                                                                                <  4:336250/149‑1 (MQ=255)
                                                                tCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCATAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAtttgtttg                                                                                >  4:370622/1‑149 (MQ=255)
                                                                  cTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAtttgtttg                                                                                >  6:40112/1‑147 (MQ=255)
                                                                  cTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTAtttgtttg                                                                                <  5:40112/147‑1 (MQ=255)
                                                                  cTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCg                                                                              <  8:149247/149‑1 (MQ=255)
                                                                   tGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCg                                                                              >  7:280660/1‑148 (MQ=255)
                                                                    gATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGa                                                                            >  1:289412/1‑149 (MQ=255)
                                                                    gATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGa                                                                            >  5:98693/1‑149 (MQ=255)
                                                                     aTGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGAt                                                                           <  8:338250/149‑1 (MQ=255)
                                                                        cAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCaa                                                                        <  1:16348/149‑1 (MQ=255)
                                                                         aGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAAt                                                                       >  3:121930/1‑149 (MQ=255)
                                                                           cGACAGTTTCGCTCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATgg                                                                     >  7:158814/1‑149 (MQ=255)
                                                                            gACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATggg                                                                    <  4:354986/149‑1 (MQ=255)
                                                                             aCAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGACGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGc                                                                   <  4:323883/149‑1 (MQ=255)
                                                                                         cAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATaa                                                       >  7:371610/1‑149 (MQ=255)
                                                                                          aGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAAc                                                      <  8:158814/149‑1 (MQ=255)
                                                                                            gAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCt                                                    >  2:364650/1‑149 (MQ=255)
                                                                                             aGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCtt                                                   <  3:408585/149‑1 (MQ=255)
                                                                                               cGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGc                                                 <  4:92645/149‑1 (MQ=255)
                                                                                                 tAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGcc                                                >  4:166778/1‑148 (MQ=255)
                                                                                                   gCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGa                                              <  8:371610/148‑1 (MQ=255)
                                                                                                     ggTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATCGAGATAACCTTGCCGATgg                                           <  1:48822/149‑1 (MQ=255)
                                                                                                         gATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGggcg                                       >  2:60620/1‑149 (MQ=255)
                                                                                                         gATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCATAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATgg                                           <  2:362488/145‑1 (MQ=255)
                                                                                                         gATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCATAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATgg                                           >  1:362488/1‑145 (MQ=255)
                                                                                                             aGACAAGCATCTGGCTGGAAATCTGTATAGCCCGCCGCCGTAGCCGCACGTTGATATAAATTATGAAAACTCTGTGGCATGGCGGGGCACGGGTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCa                                   <  5:268979/149‑1 (MQ=255)
                                                                                                             aGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGgtcgtc                                     >  8:152382/1‑147 (MQ=255)
                                                                                                             aGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCATAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCa                                   <  3:370622/149‑1 (MQ=255)
                                                                                                              gACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAg                                  <  4:121930/149‑1 (MQ=255)
                                                                                                                    cATCTGTCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATTTGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGa                            <  2:20665/149‑1 (MQ=255)
                                                                                                                     aTCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAGCTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGAt                           <  5:304533/149‑1 (MQ=255)
                                                                                                                          gctggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCca                      >  5:401900/1‑149 (MQ=255)
                                                                                                                          gctggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAG‑CCAGATGCCca                      <  8:280660/148‑1 (MQ=255)
                                                                                                                           ctggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCtt                                                   <  6:386451/119‑1 (MQ=255)
                                                                                                                           ctggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCtt                                                   >  5:386451/1‑119 (MQ=255)
                                                                                                                            tggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGgtc                                        >  6:202328/1‑129 (MQ=255)
                                                                                                                            tggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGgtc                                        <  5:202328/129‑1 (MQ=255)
                                                                                                                            tggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCca                      >  6:311645/1‑147 (MQ=255)
                                                                                                                            tggAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCca                      <  5:311645/147‑1 (MQ=255)
                                                                                                                                aaTCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTg                >  1:107946/1‑149 (MQ=255)
                                                                                                                                  tCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGt              >  1:421222/1‑149 (MQ=255)
                                                                                                                                      gATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTc             <  3:179654/146‑1 (MQ=255)
                                                                                                                                        tAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGGGATAACCTTGCCGATGGGTCGTCCAGCCCAGAAGCCCACAGTTGGTCATCGc        <  6:98693/149‑1 (MQ=255)
                                                                                                                                        tAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGGCGTCCAGCCCAGATGCCCACAGTTGGTCATCGc        <  2:289412/149‑1 (MQ=255)
                                                                                                                                           cccGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCAc     >  4:322263/1‑149 (MQ=255)
                                                                                                                                            ccgccgccgTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCAcc    >  2:332067/1‑149 (MQ=255)
                                                                                                                                              gccgccaTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGa  >  6:63002/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AAAAATCGCGGGTAAGCCCAGAGAAACAGAAACAATTGGCGCGCGCAGATCCGGTTCGTCTTTATCCTGATGCAGCGACAGTTTCGCGCCAGGAGCGTAGCGGTTGATAAGACAAGCATCTGGCTGGAAATCTGGATAGCCCGCCGCCGTAGCCGCACGTTGACATAAATTATGAAAACTCTGTGGCATGGCGGGCCACGGTTTATTTGTTTGCGGATCAATGGGCGAATAGAGATAACCTTGCCGATGGGTCGTCCAGCCCAGATGCCCACAGTTGGTCATCGCCACCGA  >  NC_000913/2308880‑2309170

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: ‑