Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,316,953 G→C 5.7% intergenic (‑113/+60) 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,9530GC 5.7% 276.5 / 6.4 87intergenic (‑113/+60)ribC/treSAriboflavin synthase/trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (43/39);  new base C (2/3);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GGGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTC  >  NC_002947/3316809‑3317093
                                                                                                                                                |                                                                                                                                            
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)
           cTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTgcgc                                                                                                                                             >  2:2261776/1‑135 (MQ=255)
           cTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTgcgc                                                                                                                                             <  1:2261776/135‑1 (MQ=255)
           cTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTg                                                                                                                              <  1:710954/150‑1 (MQ=255)
           cTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTg                                                                                                                              >  2:766301/1‑150 (MQ=255)
            tGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGc                                                                                                                             >  1:1918863/1‑150 (MQ=255)
              gACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCtt                                                                                                                           >  1:1143755/1‑150 (MQ=255)
                      cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTtgt                                                                                                                         <  2:2345772/144‑1 (MQ=255)
                      cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTtgt                                                                                                                         >  1:2345772/1‑144 (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:1604746/150‑1 (MQ=255)
                            aCATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCgg                                                                                                             <  1:1807430/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                                                                                                        >  1:2109213/1‑143 (MQ=255)
                                        tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGc                                                                                                        <  2:2109213/143‑1 (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                                                        >  1:749786/5‑120 (MQ=255)
                                                                                                               actcgcAAGCCGGCTTGCCAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGAt                                                        <  2:749786/116‑1 (MQ=255)
                                                                                                                ctggcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAgg                         <  1:529009/150‑1 (MQ=255)
                                                                                                                   gcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTc                      >  2:373965/1‑150 (MQ=255)
                                                                                                                     aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                    >  2:846711/1‑150 (MQ=255)
                                                                                                                     aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                    <  2:2112548/150‑1 (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                                                            <  1:1822377/97‑1 (MQ=255)
                                                                                                                                  aGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                            >  2:1822377/1‑97 (MQ=255)
                                                                                                                                      tGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAg                                            <  2:13730/109‑1 (MQ=255)
                                                                                                                                      tGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAg                                            >  1:13730/1‑109 (MQ=255)
                                                                                                                                       gccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAg              <  2:1219937/136‑1 (MQ=255)
                                                                                                                                       gccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTc  >  2:1819697/3‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GGGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTC  >  NC_002947/3316809‑3317093

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: