Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,566,953:1 +G 6.7% intergenic (‑2/+96) 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,9531.G 6.7% 301.1 / 11.6 75intergenic (‑2/+96)pctB/PP_2250methyl‑accepting chemotaxis protein PctB/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (37/33);  new base G (3/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGC  >  NC_002947/2566814‑2567094
                                                                                                                                            |                                                                                                                                              
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                                                                                                                        >  1:205425/1‑127 (MQ=255)
                                     aaaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCga                                                                                                                        <  2:205425/127‑1 (MQ=255)
                                          cACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGTCATGa                                                                                                                                                 <  4:346212/97‑1 (MQ=255)
                                          cACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGTCATGa                                                                                                                                                 >  3:346212/1‑97 (MQ=255)
                                           acgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGa                                                                                            <  7:129606/149‑1 (MQ=255)
                                           acgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGa                                                                                            <  1:115857/149‑1 (MQ=255)
                                            cgacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGaa                                                                                           <  7:234151/149‑1 (MQ=255)
                                             gacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                        >  7:31044/1‑135 (MQ=255)
                                             gacgCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAggg                                                                                                        <  8:31044/135‑1 (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                                                              >  3:110747/1‑147 (MQ=255)
                                                                           gCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGttt                                                              <  4:110747/147‑1 (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                                                                      >  2:331338/3‑92 (MQ=37)
                                                                                                                          cccgcACCGCGTCATGAGCGCGCTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTAc                                                                      <  1:331338/90‑1 (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)
                                                                                                                                            |                                                                                                                                              
CAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGC  >  NC_002947/2566814‑2567094

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: