Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,566,958 G→C 6.5% intergenic (‑7/+91) pctB ← / ← PP_2250 methyl‑accepting chemotaxis protein PctB/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,566,9580GC 6.5% 249.8 / 4.5 77intergenic (‑7/+91)pctB/PP_2250methyl‑accepting chemotaxis protein PctB/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base C (3/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTT  >  NC_002947/2566814‑2567104
                                                                                                                                                 |                                                                                                                                                   
cAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCg                                                                                                                                                 >  7:152587/1‑149 (MQ=255)
cAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCACCGCGTCATGAGCGCGCTGc                                                                                                                                                  <  1:371506/148‑1 (MQ=255)
  ggTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCggg                                                                                                                                               <  1:235304/149‑1 (MQ=255)
  ggTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCggg                                                                                                                                               <  2:196246/149‑1 (MQ=255)
   gTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGAGCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCgggg                                                                                                                                              >  4:6522/1‑149 (MQ=255)
       tCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCGTa                                                                                                                                          >  6:250040/1‑149 (MQ=255)
       tCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTa                                                                                                                                          >  5:71342/1‑149 (MQ=255)
                 gCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGAcc                                                                                                                                    >  3:247153/1‑145 (MQ=255)
                    tGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGatat                                                                                                                             >  8:129606/1‑149 (MQ=255)
                     gAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATc                                                                                                                            >  2:115857/1‑149 (MQ=255)
                     gAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATc                                                                                                                            >  5:355303/1‑149 (MQ=255)
                      aaCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCg                                                                                                                           >  1:173255/1‑149 (MQ=255)
                       aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCgg                                                                                                                          >  5:251375/1‑149 (MQ=255)
                       aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCgg                                                                                                                          <  8:344189/149‑1 (MQ=255)
                       aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCgg                                                                                                                          >  3:141244/1‑149 (MQ=255)
                            gCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGa                                                                                                                     <  6:251375/149‑1 (MQ=255)
                              aaaTGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAgg                                                                                                                   <  4:247153/149‑1 (MQ=255)
                                   ccAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTaa                                                                                                              <  7:232783/149‑1 (MQ=255)
                                    cAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                                  >  7:229662/1‑144 (MQ=255)
                                    cAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                                  <  8:229662/144‑1 (MQ=255)
                                    cAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTaaa                                                                                                             <  8:152587/149‑1 (MQ=255)
                                    cAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTaa                                                                                                              >  6:304501/1‑148 (MQ=255)
                                    cAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTaa                                                                                                              >  8:234151/1‑148 (MQ=255)
                                     aaaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCga                                                                                                                                  <  2:205425/127‑1 (MQ=255)
                                     aaaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCga                                                                                                                                  >  1:205425/1‑127 (MQ=255)
                                           acgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGa                                                                                                      <  7:129606/149‑1 (MQ=255)
                                           acgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGa                                                                                                      <  1:115857/149‑1 (MQ=255)
                                            cgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGaa                                                                                                     <  7:234151/149‑1 (MQ=255)
                                             gacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                                  <  8:31044/135‑1 (MQ=255)
                                             gacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                                  >  7:31044/1‑135 (MQ=255)
                                             gacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAg                                                                                                    >  8:179184/1‑149 (MQ=255)
                                                gccgcAGGCGCTGAGCCGGATCTTGTGGCGAAATTTCAGGGTCAGGTTTCGAGACCTTGAGCCACTGAGGGGAGGGCGCAGCGCGCGCATGACGCGGTGCGGGGCCTTTCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTg                                                                                                  <  8:181483/148‑1 (MQ=255)
                                                 ccgcAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCCGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCgag                                                                                                                                 >  2:130780/1‑116 (MQ=255)
                                                  cgcAGGCGCTGAGCAGGATCTTGTGTCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCtt                                                                                               <  5:304501/149‑1 (MQ=255)
                                                   gcAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTa                                                                                              <  8:103365/149‑1 (MQ=255)
                                                    cAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAg                                                                                             <  1:194672/149‑1 (MQ=255)
                                                    cAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAg                                                                                             <  1:107109/149‑1 (MQ=255)
                                                         gcTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGaa                                                                                         >  7:221111/1‑148 (MQ=255)
                                                           tGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGaa                                                                                                     >  5:129141/1‑134 (MQ=255)
                                                           tGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGaa                                                                                                     <  6:129141/134‑1 (MQ=255)
                                                                aGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTa                                                                                 <  8:221111/149‑1 (MQ=255)
                                                                 ggATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTa                                                                                 >  5:146645/1‑148 (MQ=255)
                                                                    tCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGa                                                                             >  3:280900/1‑149 (MQ=255)
                                                                        gtgGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGtt                                                                         >  3:131883/1‑149 (MQ=255)
                                                                         tgGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGttt                                                                        <  4:141244/149‑1 (MQ=255)
                                                                         tgGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGtt                                                                         <  6:355303/148‑1 (MQ=255)
                                                                           gCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGTGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTc                                                                       >  4:239481/1‑148 (MQ=255)
                                                                           gCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGttt                                                                        <  4:110747/147‑1 (MQ=255)
                                                                           gCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGttt                                                                        >  3:110747/1‑147 (MQ=255)
                                                                            cGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGAAGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCga                                                                     >  5:266007/1‑149 (MQ=255)
                                                                                 tttCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGAGCTTAGGGAAAATTCTTACAGACGTTTCgacggac                                                                >  6:268905/1‑149 (MQ=255)
                                                                                        ttCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGt                                                          >  6:15175/1‑148 (MQ=255)
                                                                                           aTGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCAc                                                       <  6:71342/148‑1 (MQ=255)
                                                                                              tttCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGAtt                                                                                                        <  3:125830/96‑1 (MQ=255)
                                                                                              tttCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGAtt                                                                                                        >  4:125830/1‑96 (MQ=255)
                                                                                              tttCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCgacggac                                                                >  7:208206/1‑136 (MQ=255)
                                                                                              tttCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCgacggac                                                                <  8:208206/136‑1 (MQ=255)
                                                                                                  gagaCCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATcagc                                               <  7:179184/149‑1 (MQ=255)
                                                                                                             ccAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGa                                    >  4:233659/1‑149 (MQ=255)
                                                                                                             ccAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGa                                    >  5:122696/1‑149 (MQ=255)
                                                                                                              cAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGGCACGGATCAGCCAGCTGCTCGAt                                   <  5:268905/149‑1 (MQ=255)
                                                                                                                   aggggagggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATcagca                              <  5:15175/149‑1 (MQ=255)
                                                                                                                         ccccgcACCGCGTCATGAGCGCGCTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTg                        >  8:338584/4‑149 (MQ=255)
                                                                                                                          ggCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTg                        >  4:34528/1‑148 (MQ=255)
                                                                                                                          cccgcACCGCGTCATGAGCGCGCTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTAc                                                                                <  1:331338/90‑1 (MQ=37)
                                                                                                                          cccgcACCGCGTCATGAGCGCGCTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTAc                                                                                >  2:331338/3‑92 (MQ=37)
                                                                                                                           gctctgcgcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGGTTGAGGGGTAAAAGATTGAAGTTATTAGGGAAAATTCTTACAGACGTATCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCa                      <  3:6522/144‑1 (MQ=255)
                                                                                                                           gcgcAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCGTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTg                        >  6:328943/1‑147 (MQ=255)
                                                                                                                            cgcACCGCGTCATGAGCGCGCTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAg                    >  1:44200/1‑149 (MQ=255)
                                                                                                                                gcgcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctg                >  5:389874/1‑149 (MQ=255)
                                                                                                                                 cgcgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctgc               <  2:173255/149‑1 (MQ=255)
                                                                                                                                   cgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCg                                                                      >  4:178140/1‑93 (MQ=255)
                                                                                                                                   cgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCg                                                                      <  3:178140/93‑1 (MQ=255)
                                                                                                                                   cgcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCctgctg             >  5:51852/1‑149 (MQ=255)
                                                                                                                                    gcTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCTGCACGCctgctga            <  6:122696/149‑2 (MQ=255)
                                                                                                                                             cgcgGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtg    <  3:239481/149‑1 (MQ=255)
                                                                                                                                               cgGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCtgtt  <  4:280900/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGCAGCGGCTGTT  >  NC_002947/2566814‑2567104

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