Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,324,132 T→G 9.4% intergenic (‑104/+88) treSA ← / ← PP_2919 trehalose synthase A/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,324,1320TG 9.4% 763.3 / 4.1 291intergenic (‑104/+88)treSA/PP_2919trehalose synthase A/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (111/130);  new base G (33/4);  total (149/142)
Fisher's exact test for biased strand distribution p-value = 5.91e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTGCGCACCG  >  NC_002947/3323983‑3324281
                                                                                                                                                     |                                                                                                                                                     
cgcgGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCt                                                                                                                                                       >  2:1405535/1‑150 (MQ=255)
 gcgGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTc                                                                                                                                                      <  1:6809348/150‑1 (MQ=255)
    gTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                                   <  1:2389151/150‑1 (MQ=255)
    gTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCGTGCCCCCGcccg                                                                                                                                                   >  1:6415147/1‑149 (MQ=255)
    gTCTTCGAGCCAATTGACGTATGACGGGTCGGGCTGGGACATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCca                                                                                                                                                   <  1:5717911/150‑1 (MQ=255)
      cTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                                 <  1:5815216/150‑1 (MQ=255)
       ttCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTTCCCGCTCCCACAg                                                                                                                                                >  2:3993643/1‑150 (MQ=255)
       ttCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTTCCCCCGCCCACAg                                                                                                                                                >  1:1265366/1‑150 (MQ=255)
       ttCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTTCCCCCACCcacac                                                                                                                                                >  1:5160534/1‑149 (MQ=255)
       ttCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCCCAg                                                                                                                                                >  2:5946874/1‑150 (MQ=255)
           aGCCATTTGACTTATGACTGGTCTTGATGGGTCAATTCGCGACCTCTTTGCGTGTCGGCACCGGCTCAGGATATGACTGCGCTGCCGCCTGAACTTTAGGCTCAAGATGGTTGTTGGTTTCTTCGCGGGCTTGCCAGATCACACAGGccc                                                                                                                                            >  2:4179668/1‑148 (MQ=255)
            gCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGa                                                                                                                                           <  1:1792036/150‑1 (MQ=255)
            gCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGa                                                                                                                                           <  2:3250920/150‑1 (MQ=255)
            gCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCCCTCCCACAGGCggc                                                                                                                                           >  1:6039970/1‑149 (MQ=255)
                  tGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCCCAGGCGGACATTcc                                                                                                                                     >  2:5054284/1‑150 (MQ=255)
                  tGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATAcc                                                                                                                                     <  2:6883386/150‑1 (MQ=255)
                    aCGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGc                                                                                                                                   <  1:3378560/150‑1 (MQ=255)
                     cGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCt                                                                                                                                  <  1:1036830/150‑1 (MQ=255)
                      gTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTg                                                                                                                                 >  2:4723383/1‑150 (MQ=255)
                      gTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCGTGCCCGCTCCCACAGGCGGACAGTCCGcgg                                                                                                                                 >  1:7298569/1‑148 (MQ=255)
                        aTGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGaa                                                                                                                               <  2:1971351/150‑1 (MQ=255)
                        aTGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGaa                                                                                                                               <  1:1499916/150‑1 (MQ=255)
                         tGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCACGCGCCCACAGGCGGACAGACCGCGGAAg                                                                                                                              >  2:342181/1‑150 (MQ=255)
                               gTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTgc                                                                                                                        <  2:2416286/150‑1 (MQ=255)
                                  gggCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCt                                                                                                                     <  2:2840472/150‑1 (MQ=255)
                                  gggCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCt                                                                                                                     <  2:7186757/150‑1 (MQ=255)
                                  gggCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCCCTCCCACAGGCGGCCAAACCGCGGAAGATTTgcgcg                                                                                                                     >  2:3175758/1‑149 (MQ=255)
                                   ggCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                                 >  2:1259333/1‑121 (MQ=255)
                                     cTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCa                                                                                                                  <  1:7221898/150‑1 (MQ=255)
                                     cTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCGCCCACACGCGGACAGACCGCGGAAGGTTTGCGCTGCa                                                                                                                  >  2:1647427/1‑150 (MQ=255)
                                          tCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCtgt                                                                                                             <  2:3460311/150‑1 (MQ=255)
                                          tCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCtgt                                                                                                             <  1:6212173/150‑1 (MQ=255)
                                            aTTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTgg                                                                                                           <  1:3450146/150‑1 (MQ=255)
                                            aTTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTgg                                                                                                           <  2:5176007/150‑1 (MQ=255)
                                             tttCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGTCATACCGCTGAAGATTTGCGCTGCACCTGgggg                                                                                                          >  2:4595823/1‑150 (MQ=255)
                                              ttCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCCCACGCGGACCAACCGCGGAAGCTTTGCGCTGCACCAGTGGGa                                                                                                         >  2:6023922/1‑150 (MQ=255)
                                              cggtgcgcaattacttccatgatgcgtttctggcaccgcagcagttgtcggcccggcgctgatcctgcctttgcaggcccggccTCTTCGCGGGCGCGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGa                                                                                                         >  1:1524141‑M1/85‑150 (MQ=255)
                                               tCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAg                                                                                                        <  1:3090055/150‑1 (MQ=255)
                                               ggtgcgcaattacttccatgatgcgtttctggcaccgcagcagttgtcggcccggcgctgatcctgcctttgcaggcccggccTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAg                                                                                                        >  2:5794720‑M1/84‑150 (MQ=255)
                                                gtgcgcaattacttccatgatgcgtttctggcaccgcagcagttgtcggcccggcgctgatcctgcctttgcaggcccggccTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGc                                                                                                       >  2:7376382‑M1/83‑150 (MQ=255)
                                                 tgcgcaattacttccatgatgcgtttctggcaccgcagcagttgtcggcccggcgctgatcctgcctttgcaggcccggccTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCg                                                                                                      >  1:987212‑M1/82‑150 (MQ=255)
                                                 tgcgcaattacttccatgatgcgtttctggcaccgcagcagttgtcggcccggcgctgatcctgcctttgcaggcccggccTCTTCGCGGGCGCGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCg                                                                                                      >  2:4495113‑M1/82‑150 (MQ=255)
                                                  cgTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCgg                                                                                                     <  1:6422731/150‑1 (MQ=255)
                                                  cgTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCgg                                                                                                     <  2:3000581/150‑1 (MQ=255)
                                                   gTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCGCCCACAGGCGGAAATTACGCAGAAGATTTGCGCGGCACCGGTGGGAGCggg                                                                                                    >  2:2101617/1‑150 (MQ=255)
                                                    tCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCCCTGCACCGGTGGGAGCGGgc                                                                                                   >  2:6381062/1‑150 (MQ=255)
                                                        ctTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCATCGCGGGCTTGCCCCCACCCACCGGCCGACATACCGCAGAAGAGATGCGCTGCACCGGTGGGACCGGgcgcgc                                                                                               >  2:5749125/1‑150 (MQ=255)
                                                        ctTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGATCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCCCAGGCGGCCATAACGCGTAAGGATTGGGCTGCACCTGTGGGAGCGGgcgcgc                                                                                               >  2:4075183/1‑150 (MQ=255)
                                                          ttGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACCTACCGCTGAAGATTTGCGCTGCCCCGGGGGGAGCGGGCGCGccc                                                                                             >  2:6703035/1‑150 (MQ=255)
                                                            gCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGGTTGCCCGCGCCCCCAAGCCGACATAACGCGGAAAAATTTGGCGGGACCAGTGGGTGCGGGGCCGccccc                                                                                           >  1:4639427/1‑148 (MQ=255)
                                                            gCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCGTGCCCGCGCCCCCAGGCGGCCATACCGCGGCAGAGGTGCGCGGCCCCGGGGGGCGCGGGCGCGCCcgc                                                                                           >  2:6437670/1‑150 (MQ=255)
                                                              gtgtCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCATGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAATATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCcgctc                                                                                         <  1:5413159/150‑3 (MQ=255)
                                                               tgtCGGCACCGGCTCAGGTTAGGACGGCGCTGCGGCCTGACCGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                        <  2:5632022/150‑1 (MQ=255)
                                                                 tCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAaga                                                                                      <  2:3710601/150‑1 (MQ=255)
                                                                   ggCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTTGGCGCGCCCGAAGGCGGGCAGAGCGCGGAATTTTTTGGGTGCAGCAGTGGGGGCGGGGGGGGCCGCGAAGAgg                                                                                    >  1:396788/1‑150 (MQ=255)
                                                                   ggCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCGGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAgg                                                                                    >  2:7247359/1‑150 (MQ=255)
                                                                    gCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGc                                                                                   <  1:2967367/150‑1 (MQ=255)
                                                                    gCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCGCCCACAGGCGGACATACCGCAGAAGATTTGCGCTGCACCAGTGGGAGCGGGCGCGCCCGCGCAGAGGc                                                                                   >  2:2562194/1‑150 (MQ=255)
                                                                    gCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCATTCCCGCTCCCCCAGGCGGGCAAACCGCGGTAGGTTTGCGGGGCACCAGTGGGAGCGGGCGCGCCCGCGAAGGGGc                                                                                   >  2:2555041/1‑150 (MQ=255)
                                                                    gCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGGTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGTCATTCCGCTGAAGTTTTGGGCGGCGCCAGTGGGAGCGGGCGCGCCCCCGAAGTGGc                                                                                   >  2:6606306/1‑150 (MQ=255)
                                                                     cACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCCCTCCCCCCGGCGGACATACCGCGGAAGATTTGCGCTGCACCAGTGGGCGCGGGCGCGCCCGCGAACAGGcc                                                                                  >  2:2668612/1‑150 (MQ=255)
                                                                       ccGGCTCTGGTTTTGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTTTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCGGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCgg                                                                                <  2:6945993/150‑1 (MQ=255)
                                                                                 ttttGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGCTGGGTGTTGGTTTCTTCGCGGGCTGGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAg                                                                      <  1:4285063/147‑1 (MQ=255)
                                                                                         gcgcTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCGGACGTTTTGCGCGGCACCGGGGGGAGCGGGCGCGCCCGCGCAGCGGCCGGGCCTGCACAGCCTGGATc                                                              >  2:6056113/1‑150 (MQ=255)
                                                                                            cTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAgc                                                           <  1:1143712/150‑1 (MQ=255)
                                                                                             tGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCGCCCCCACGCCGACATGCCGCGGAAGGTTTGCGGTGCACCAGGGGGAGCGGGCGCGCCCGCGGAGAGGCCGGGCCGGCCAAGGCAGGATCAgcg                                                          >  2:1110284/1‑150 (MQ=255)
                                                                                                   cTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGc                                                    <  2:4829460/150‑1 (MQ=255)
                                                                                                    tGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCGTTCCCGCTCCCACAGGCGGACATAACCCTGAAGATTTTCGCTGCACCGGTGGGAGCGGGCGCGCCCGCGCAGAGGCCGGGCCAGCAAAGGCACGATAAGCGCCCGGcc                                                   >  2:1255244/1‑150 (MQ=38)
                                                                                                      aaCGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAGGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGa                                                 <  1:3622382/150‑1 (MQ=255)
                                                                                                          ttCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCGCCCACAGGCGGCCATCCCGCGGAAGAGTTGCGCTGCACCAGTGGGAGCGGGCGCGCCCGCGCACACGCCGGGCCAGCCAACGCAGGATCACCGCCCGGCCGCCAAc                                             >  2:7132291/1‑150 (MQ=255)
                                                                                                            cGGCTCACGTGGGGTGTTGGTTTCTGCGCGGGGTGCCCCGCCCCCACTGGCGGACTCCCCGCTGACGTTTGGCGCGGCACGGGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAActg                                           <  1:67004/150‑1 (MQ=255)
                                                                                                                tCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACtgctgc                                       <  2:3590441/150‑1 (MQ=255)
                                                                                                                   agtggggTGGTGGTTTCTTCGCGGGCTGGCCCGCCCCCACAGGCGGACATCCCGCTTAAGATTGGCGCGGCACCTGTGGGCGGGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGt                                    <  1:5935175/146‑1 (MQ=255)
                                                                                                                       gggTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCa                                <  2:4031202/150‑1 (MQ=255)
                                                                                                                       gggTGTTGGTTTCTTCGCGGGCTTGCCCCCTCCCACACGCCGACATACCCCTGAAGATTTGCGCTGCACCAGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCAGCCAACGCAGGATCAGCGCCGGGCCCACAACTGCTGCGGGGCCa                                >  2:5894030/1‑150 (MQ=255)
                                                                                                                          tgtTGGTTTCTTCGCGGGCTTGCCCCCGCCCACAGGCCGACATGCCGCAGAATTTTTGCGCTGCACCTTTGGGAGCGGGCGCGCCCGCGCAGACGCCGGGCAAGCAAAGGCAGGATCAGCGCCCCGCCGACAACACCAGCCCTGCCAGaa                             >  2:5548157/1‑150 (MQ=255)
                                                                                                                            ttGGTTTCTTCGCGGGCTTTCCCGCTCCCACACGCGGACCAACCCCGGAAGATTTGCGCAGCAACAGTGGGAGCGGGCGCGCCCCCGAAGAGGCCGGGCCAGCAAACGCAGGATCACCGCCCGGCCGACAACTCCAGCGGTGCCAGAAAc                           >  2:3892353/1‑150 (MQ=255)
                                                                                                                            ttGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCGGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCACGATCACCGCCCGGCCGACAACAGCAGCGGGGCCACAAAc                           >  2:5292934/1‑150 (MQ=38)
                                                                                                                             tGGTTTCTTCGCGGGCTTGCCCGCTCCCACCGGCGGACATACCGCTGAAGATTTGCGCGGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACg                          <  2:153088/150‑1 (MQ=255)
                                                                                                                               gTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGca                        <  1:3993535/150‑1 (MQ=255)
                                                                                                                               gTTTCTTCGCGGGCTTGCCCCCGCCCCCAGCCCGACATACCGCCCAAGGGTTGCGCGGCACCAGGGGGTGCGGGCCCGCCCCCCGACCCGCCCGGCCCGCCAAGCCACGATCCCCGCCCGGCCGACAACCGCACCCCAGCCACAAACGca                        >  2:4268490/1‑150 (MQ=255)
                                                                                                                                 ttcttcGGGGGCGTGCCCGCTCCCACTGGCGGACATACCGCTTAAGATTTGCGCGGCACGTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGcatc                      <  1:4737937/150‑1 (MQ=255)
                                                                                                                                  TCTTCGCGGGCTTGCCCGCTCCCCCAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGcatca                     >  2:6473754‑M1/1‑150 (MQ=255)
                                                                                                                                  TCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATCCCGCTGAAGATTTGCGCGGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGcatca                     >  2:683520‑M1/1‑150 (MQ=255)
                                                                                                                                       ggGGGGTTGCCCGCTCCCACTGGGGTACATACCGCTGAAGTTTTGCGCTGCACGTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa                >  1:1405535‑M1/3‑150 (MQ=255)
                                                                                                                                       GCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa                <  2:5257020‑M1/150‑1 (MQ=255)
                                                                                                                                       GCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa                <  2:838143‑M1/150‑1 (MQ=255)
                                                                                                                                       GCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa                >  2:6236810‑M1/1‑150 (MQ=255)
                                                                                                                                       GCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGCGCCGACAACTGCTGCGGTGCCAGAAACGCATCAAGGaa                <  1:4147791‑M1/150‑1 (MQ=255)
                                                                                                                                       GCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCAGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCAAGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTCCGGAGCCAGAAACGCCTCAGGGaa                <  2:4872452‑M1/150‑1 (MQ=255)
                                                                                                                                             ttGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAAtt          >  2:3360369/1‑150 (MQ=255)
                                                                                                                                              ggCCCGCGCCCACCGGCGGACACACCGCTGAAGATTTGCGCGGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTg         <  1:5764139/149‑1 (MQ=255)
                                                                                                                                                cccGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGACGTAATTgcg       >  2:1621515/1‑150 (MQ=255)
                                                                                                                                                cccGCGCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTgcg       <  1:5069526/150‑1 (MQ=255)
                                                                                                                                                  cGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTgcgcc     >  1:6240319/1‑149 (MQ=255)
                                                                                                                                                  cGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGAGCCAGAAACGCATCATGGAAGTAATTGCGCa     >  1:2490637/1‑150 (MQ=255)
                                                                                                                                                   gCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAaga                                                                                      >  2:2781265/1‑68 (MQ=255)
                                                                                                                                                   gCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAaga                                                                                      <  1:2781265/68‑1 (MQ=255)
                                                                                                                                                   gCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTGCGCAc    >  1:3845602/1‑150 (MQ=255)
                                                                                                                                                     tCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTGCGCACCg  >  1:4450946/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGGTCTTCGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTGCGCACCG  >  NC_002947/3323983‑3324281

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