Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,316,957 G→C 5.9% intergenic (‑117/+56) ribC ← / ← treSA riboflavin synthase/trehalose synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,316,9570GC 5.9% 261.5 / 6.0 84intergenic (‑117/+56)ribC/treSAriboflavin synthase/trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base C (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

GGGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAATCGCCT  >  NC_002947/3316809‑3317106
                                                                                                                                                    |                                                                                                                                                     
gggCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCgggg                                                                                                                                                      <  2:1309846/150‑1 (MQ=255)
    ggACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAgg                                                                                                                                                  >  2:710954/1‑150 (MQ=255)
    ggACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAgg                                                                                                                                                  >  2:2394930/1‑150 (MQ=255)
        ggCTTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAGCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAggcggc                                                                                                                                              <  1:874117/150‑1 (MQ=255)
         gCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGCCGCCTGCCCCGCGCAggcgcca                                                                                                                                             <  1:791846/150‑4 (MQ=255)
          ccTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGt                                                                                                                                            >  2:2261433/1‑150 (MQ=255)
           cTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTg                                                                                                                                           >  2:766301/1‑150 (MQ=255)
           cTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTg                                                                                                                                           <  1:710954/150‑1 (MQ=255)
            tGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGc                                                                                                                                          >  1:1918863/1‑150 (MQ=255)
              gACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCtt                                                                                                                                        >  1:1143755/1‑150 (MQ=255)
                      cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTtgt                                                                                                                                      >  1:2345772/1‑144 (MQ=255)
                      cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTtgt                                                                                                                                      <  2:2345772/144‑1 (MQ=255)
                      cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGa                                                                                                                                >  1:1951317/1‑150 (MQ=255)
                       ggTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGACCGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCCCCGCGCAggcgcca                                                                                                                                             <  2:72941/136‑4 (MQ=255)
                       ggTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCCCCGCGCAggcgcca                                                                                                                                             >  1:72941/1‑133 (MQ=255)
                        gtgtACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGACGTAAACGCTATATTAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGaaa                                                                                                                              <  1:1482831/150‑1 (MQ=255)
                           tACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGGGGCGTGCTTGTGGGCGAAAGCg                                                                                                                           <  2:1911321/150‑1 (MQ=255)
                            aCATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCgg                                                                                                                          <  1:1807430/150‑1 (MQ=255)
                            aCATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCgg                                                                                                                          <  1:1604746/150‑1 (MQ=255)
                             cATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGt                                                                                                                         >  2:2290497/1‑150 (MQ=255)
                              aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCgggg                                                                                                                                                      <  1:2143054/120‑1 (MQ=255)
                              aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCgggg                                                                                                                                                      >  2:2143054/1‑120 (MQ=255)
                               tGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGttt                                                                                                                       >  2:173038/1‑150 (MQ=255)
                                ggAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCCGGCGGCGTGCTTGTGGGCGAAACCGGTTTg                                                                                                                      <  1:2290497/150‑1 (MQ=255)
                                        tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGc                                                                                                                     <  2:2109213/143‑1 (MQ=255)
                                        tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGc                                                                                                                     >  1:2109213/1‑143 (MQ=255)
                                        tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCa                                                                                                              <  1:173038/150‑1 (MQ=255)
                                        tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCa                                                                                                              >  2:529009/1‑150 (MQ=255)
                                          cTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATg                                                                                                            <  1:1888633/150‑1 (MQ=255)
                                           tGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGATGTAAATGCAAGACCAGGCGTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGAAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATgg                                                                                                           <  2:2365628/150‑1 (MQ=255)
                                              gCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGATTGCGACGGCATGGCCa                                                                                                        >  1:930889/1‑150 (MQ=255)
                                                   gaaaagaaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCa                                                                                                   <  2:2164630/150‑1 (MQ=255)
                                                      aagaaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAggcggc                                                                                                                                              <  1:2307861/104‑1 (MQ=255)
                                                      aagaaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAggcggc                                                                                                                                              >  2:2307861/1‑104 (MQ=255)
                                                        gaaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTc                                                                                              >  1:1580694/1‑150 (MQ=255)
                                                              ttCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTgg                                                                                                                                    >  2:1799495/1‑106 (MQ=255)
                                                              ttCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTgg                                                                                                                                    <  1:1799495/106‑1 (MQ=255)
                                                                ccGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATg                                                                                       <  2:971401/149‑1 (MQ=255)
                                                                ccGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATg                                                                                       >  1:971401/1‑149 (MQ=255)
                                                                  gAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGccc                                                                                    >  2:912330/1‑150 (MQ=255)
                                                                  gAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGccc                                                                                    >  1:1389435/1‑150 (MQ=255)
                                                                      gCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCgctg                                                                                <  2:1951317/150‑1 (MQ=255)
                                                                           gCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGtt                                                                           <  1:2354059/150‑1 (MQ=255)
                                                                            cccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTg                                                                          >  2:1197788/1‑150 (MQ=255)
                                                                             ccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGTTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         <  1:1154295/150‑1 (MQ=255)
                                                                             ccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAATATGCCCGCTGCTGTTGa                                                                         >  2:1265638/1‑150 (MQ=255)
                                                                             ccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         >  2:1054167/1‑150 (MQ=255)
                                                                             ccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTCGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGGGATTAGAAGCAGCGACGGCATGACCGCGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         >  2:1243394/1‑150 (MQ=255)
                                                                                 aGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGAt                                                                     <  1:2261433/150‑1 (MQ=255)
                                                                                    gcAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCg                                                                  >  2:2219322/1‑150 (MQ=255)
                                                                                     cAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGc                                                                 <  1:1197788/150‑1 (MQ=255)
                                                                                     cAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGc                                                                 <  2:1918863/150‑1 (MQ=255)
                                                                                      aGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCa                                                                <  1:1831398/150‑1 (MQ=255)
                                                                                        tAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAgg                                                              >  1:740229/1‑150 (MQ=255)
                                                                                          aaTGTAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGca                                                            >  2:1759225/1‑150 (MQ=255)
                                                                                            tGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGcaca                                                          >  1:455306/1‑150 (MQ=255)
                                                                                            tGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTCGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGcaca                                                          >  2:424457/1‑150 (MQ=255)
                                                                                                     aGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGt                                                 >  2:51929/1‑150 (MQ=255)
                                                                                                         ttCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATg                                                                                                            <  1:623950/87‑1 (MQ=255)
                                                                                                         ttCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATg                                                                                                            >  2:623950/1‑87 (MQ=255)
                                                                                                             tggctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         >  2:2288924/1‑118 (MQ=255)
                                                                                                             tggctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         <  1:2288924/118‑1 (MQ=255)
                                                                                                              ggctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGccc                                                                                    <  1:1058064/106‑1 (MQ=255)
                                                                                                              ggctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGccc                                                                                    >  2:1058064/1‑106 (MQ=255)
                                                                                                               gctggcAAGCCGGCTTGCGCGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAg                                       >  2:882379/1‑150 (MQ=255)
                                                                                                               actcgcAAGCCGGCTTGCCAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGAt                                                                     <  2:749786/116‑1 (MQ=255)
                                                                                                               actcgcAAGCCGGCTTGCCAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGAt                                                                     >  1:749786/5‑120 (MQ=255)
                                                                                                                ctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAgg                                      <  1:529009/150‑1 (MQ=255)
                                                                                                                   gcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTc                                   >  2:373965/1‑150 (MQ=255)
                                                                                                                     aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                                 <  2:2112548/150‑1 (MQ=255)
                                                                                                                     aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                                 >  2:846711/1‑150 (MQ=255)
                                                                                                                     aaGCCGGCTTGCGAGTCTGGCGCCGGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTGGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCGCGTACGGGTCCAGGTTCAc                                 >  2:1074863/1‑150 (MQ=255)
                                                                                                                     aaGCCGGCTTGCGAGTCTGGCGCCCGCGCGGGGGAGGCGGCGTGCTGGTGGGCGAAACCGGTTTGCGACGGCACCGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACTGCGCCTCGTACGGGTCCAGGTTCAc                                 >  2:240491/1‑150 (MQ=255)
                                                                                                                      aGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGAACACCGCCTCGTACGGGTCCAGGTTCAcc                                <  1:1243394/150‑1 (MQ=255)
                                                                                                                       gCCGGCTTGCGAGGCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCa                               >  2:2311987/1‑150 (MQ=255)
                                                                                                                        ccGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCAt                              >  2:830565/1‑150 (MQ=255)
                                                                                                                         cGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAATATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATc                             <  1:1265638/150‑1 (MQ=255)
                                                                                                                          ggCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGCTTCAAACATGCCCGCTGCTGTTGACGATCCGCACGCACAGCGCCTCGTTCTGGTCCAGGTTCACCATCa                            <  2:2104418/150‑1 (MQ=255)
                                                                                                                           gctcgCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAg                           <  1:424457/146‑1 (MQ=255)
                                                                                                                             ttGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGgc                         <  2:1159560/150‑1 (MQ=255)
                                                                                                                             ttGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGgc                         <  1:373965/150‑1 (MQ=255)
                                                                                                                                  aGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         >  2:1822377/1‑97 (MQ=255)
                                                                                                                                  aGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                         <  1:1822377/97‑1 (MQ=255)
                                                                                                                                      tGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAg                                                         >  1:13730/1‑109 (MQ=255)
                                                                                                                                      tGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAg                                                         <  2:13730/109‑1 (MQ=255)
                                                                                                                                       gccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAg                           <  2:1219937/136‑1 (MQ=255)
                                                                                                                                       gccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTc               >  2:1819697/3‑150 (MQ=255)
                                                                                                                                                    ggCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCAAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAATCGCCt  >  2:1846719/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAATCGCCT  >  NC_002947/3316809‑3317106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑