Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,052,897 C→G 6.5% G208R (GGC→CGC)  ubiH ← 2‑octaprenyl‑6‑methoxyphenol hydroxylase, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,052,8970CG 6.5% 320.2 / 12.2 93G208R (GGC→CGC) ubiH2‑octaprenyl‑6‑methoxyphenol hydroxylase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base G (3/3);  total (47/46)
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.39e-01

ATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGCTGCTA  >  NC_000913/3052749‑3053043
                                                                                                                                                    |                                                                                                                                                  
aTTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGcc                                                                                                                                                    >  4:415501/1‑149 (MQ=255)
 ttttCCCAAGTCGCCAGCCAAAGGCCGACTGGCGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACAGCGGCAGCGGGCCa                                                                                                                                                   >  4:774635/1‑149 (MQ=255)
  tttCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCa                                                                                                                                                   >  2:155502/1‑148 (MQ=255)
    tCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATGGCCCGCTGGCGATg                                                                                                                                                 <  2:143518/148‑1 (MQ=255)
     cccAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGt                                                                                                                                                <  1:185970/148‑1 (MQ=255)
      ccAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTg                                                                                                                                              >  4:557339/1‑149 (MQ=255)
        aaGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCg                                                                                                                                            <  4:535848/149‑1 (MQ=255)
         aGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGt                                                                                                                                           <  3:832274/149‑1 (MQ=255)
            cGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTaa                                                                                                                                         <  2:415223/148‑1 (MQ=255)
             gccagccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAg                                                                                                                                       >  3:84486/1‑149 (MQ=255)
                 gccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGtt                                                                                                                                   <  4:696588/149‑1 (MQ=255)
                 gccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGt                                                                                                                                    <  2:800703/148‑1 (MQ=255)
                  ccaAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGtt                                                                                                                                   <  2:550828/148‑1 (MQ=255)
                    aaaGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCa                                                                                                                                 >  2:889929/1‑148 (MQ=255)
                        gCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgc                                                                                                                            <  4:84486/149‑1 (MQ=255)
                        gCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgc                                                                                                                            >  3:204358/1‑149 (MQ=255)
                          cGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgcgc                                                                                                                          <  3:867269/149‑1 (MQ=255)
                            aCTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAgcgcgc                                                                                                                        >  3:754778/1‑149 (MQ=255)
                              tagAGTTCACGGCAAAACTTCTCGTCACTCCACGTCAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGccc                                                                                                                      <  3:739799/147‑1 (MQ=255)
                                gAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGCCCATGTTGCGAAAAGCGTTCAAAAGGGCGCCCtc                                                                                                                    >  4:323101/1‑148 (MQ=255)
                                 aGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCtt                                                                                                                    <  2:425570/148‑1 (MQ=255)
                                   ttCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCAt                                                                                                                 >  3:30019/1‑149 (MQ=255)
                                    tCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCAt                                                                                                                 >  1:531408/1‑148 (MQ=255)
                                      aCGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATgcg                                                                                                              >  3:744409/1‑149 (MQ=255)
                                         gCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCaa                                                                                                           >  3:443408/1‑149 (MQ=255)
                                         gCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCa                                                                                                            >  2:301923/1‑148 (MQ=255)
                                           aaaaCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCGACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACg                                                                                                         >  4:219360/1‑149 (MQ=255)
                                             aaCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGa                                                                                                       >  3:430919/1‑149 (MQ=255)
                                              aCTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGa                                                                                                       >  2:11426/1‑148 (MQ=255)
                                                ttCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGt                                                                                                    >  3:111794/1‑149 (MQ=255)
                                                 tctcGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTa                                                                                                   <  4:111794/149‑1 (MQ=255)
                                                 tctcGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTa                                                                                                   <  4:267508/149‑1 (MQ=255)
                                                  ctcGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAg                                                                                                  <  4:430919/149‑1 (MQ=255)
                                                      tCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCaa                                                                                               >  2:399054/1‑148 (MQ=255)
                                                      tCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCaa                                                                                               >  2:571631/1‑148 (MQ=255)
                                                       cACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGaa                                                                                                      >  1:598063/1‑140 (MQ=255)
                                                       cACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGaa                                                                                                      <  2:598063/140‑1 (MQ=255)
                                                            cACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTgg                                                                                         >  2:791334/1‑148 (MQ=255)
                                                              cGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCaa                                                                                      <  3:415501/149‑1 (MQ=255)
                                                              cGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCa                                                                                       >  2:319186/1‑148 (MQ=255)
                                                               gACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCaa                                                                                      <  1:791334/148‑1 (MQ=255)
                                                               gACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCaa                                                                                      >  1:714586/1‑148 (MQ=255)
                                                               gACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCaa                                                                                      >  1:24214/1‑148 (MQ=255)
                                                                acaacaCCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAAt                                                                                     >  2:437779/1‑148 (MQ=255)
                                                                acaacaCCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATc                                                                                    <  3:819293/149‑1 (MQ=255)
                                                                   acacCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCAc                                                                                  >  1:212324/1‑148 (MQ=255)
                                                                    cacCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACg                                                                                 <  2:680325/148‑1 (MQ=255)
                                                                       ctctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGtt                                                                                           >  3:771611/1‑135 (MQ=255)
                                                                       ctctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCCGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAAAGtt                                                                                           <  4:771611/135‑1 (MQ=255)
                                                                        tctTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATAGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAg                                                                            >  4:756919/1‑149 (MQ=255)
                                                                              cgcCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAgttgtt                                                                       >  2:614404/1‑148 (MQ=255)
                                                                              cgcCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAgttgtt                                                                       <  1:345831/148‑1 (MQ=255)
                                                                              cgcCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAgttgtt                                                                       <  1:708907/148‑1 (MQ=255)
                                                                               gcCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTc                                                                      >  1:709874/1‑148 (MQ=255)
                                                                                  gTTCCAGTGGATGACACCAAGACAGCGAACAGCGTCCGGCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTaa                                                                  <  3:323101/149‑1 (MQ=255)
                                                                                           gATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGa                                                                                                       <  4:888254/103‑1 (MQ=255)
                                                                                           gATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGa                                                                                                       >  3:888254/1‑103 (MQ=255)
                                                                                            aTGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCctgc                                                         <  1:190723/148‑1 (MQ=255)
                                                                                            aTGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCctgc                                                         <  2:82105/148‑1 (MQ=255)
                                                                                               acacCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCATGGCCAGTTGTTCGTAAGGCTCCTGCTGcc                                                     <  3:774635/149‑1 (MQ=255)
                                                                                                 accagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAg                                                   >  3:655367/1‑149 (MQ=255)
                                                                                                   cagaccagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTc                                                 <  4:443408/149‑1 (MQ=255)
                                                                                                      accagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCaa                                               <  2:861005/148‑1 (MQ=255)
                                                                                                      accagCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCaa                                               <  2:593899/148‑1 (MQ=255)
                                                                                                         agCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAcgcc                                           <  4:744409/149‑1 (MQ=255)
                                                                                                          gCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAcgccg                                          >  3:396202/1‑149 (MQ=255)
                                                                                                           cGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATc                                                                                    <  2:784833/106‑1 (MQ=255)
                                                                                                           cGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATc                                                                                    >  1:784833/1‑106 (MQ=255)
                                                                                                           cGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAAcgccg                                          >  2:253972/1‑148 (MQ=255)
                                                                                                               cAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCAcg                                      <  1:127545/148‑1 (MQ=255)
                                                                                                                  cgtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCgg                                   <  1:301923/148‑1 (MQ=255)
                                                                                                                   gtccgtcAGACATCGGCGACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTg                                 <  3:219360/149‑1 (MQ=255)
                                                                                                                   gtccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTg                                 >  4:730958/1‑149 (MQ=255)
                                                                                                                    tccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCa                                                                                                <  1:172237/85‑1 (MQ=255)
                                                                                                                    tccgtcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCa                                                                                                >  2:172237/1‑85 (MQ=255)
                                                                                                                        tcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTaa                            >  3:193947/1‑149 (MQ=255)
                                                                                                                        tcAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTa                             >  2:367262/1‑148 (MQ=255)
                                                                                                                         cAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAACAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTaa                            >  1:394124/1‑148 (MQ=255)
                                                                                                                          aGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACg                          >  3:203994/1‑149 (MQ=255)
                                                                                                                             tttacgCAACATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCt                        >  2:279867/6‑148 (MQ=255)
                                                                                                                              aTCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGa                      >  3:186341/1‑149 (MQ=255)
                                                                                                                               tCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGaa                     <  4:669860/149‑1 (MQ=255)
                                                                                                                               tCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGa                      <  1:614404/148‑1 (MQ=255)
                                                                                                                                   cAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATgg                  <  1:220034/148‑1 (MQ=255)
                                                                                                                                    aaCATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGt                                  <  4:181458/131‑1 (MQ=255)
                                                                                                                                    aaCATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGt                                  >  3:181458/1‑131 (MQ=255)
                                                                                                                                     aCATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCgg                                   >  1:44120/1‑129 (MQ=255)
                                                                                                                                     aCATGGCCCGCTGGCGATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCgg                                   <  2:44120/129‑1 (MQ=255)
                                                                                                                                             aGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgc       <  3:22185/149‑1 (MQ=255)
                                                                                                                                             aGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAg        <  1:293674/148‑1 (MQ=255)
                                                                                                                                              gCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgc       <  1:11426/148‑1 (MQ=255)
                                                                                                                                              gCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGGTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAgct      <  3:557339/149‑1 (MQ=255)
                                                                                                                                                  gCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGCTGCTa  >  4:76490/1‑149 (MQ=255)
                                                                                                                                                  gCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGCTGCTa  >  4:479033/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATTTTCCCAAGTCGCCAGCCAAAGGCCGACTGGAGTTCACGGCAAAACTTCTCGTCACTCCACGACAACACCTCTTCGCGCCGTTCCAGTGGATGACACCAGACCAGCGAACAGCGTCCGTCAGACATCGGCAACATCGCCAGCGGGCCATGTTGCGTAAAGCGTTCAAAAGCGCGCCCTTCATGCGCAACGGAAGTAGCAACGTTGGCAATCACGGCCAGTTGTTCGTAAGGCTCCTGCTGCCAGTCAACGCCGCACGCGGTGGCTAACGCTGAATGGGTGCCATCAGCTGCTA  >  NC_000913/3052749‑3053043

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑