Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,566,947 Δ1 bp 9.9% coding (5/1932 nt) pctB ← methyl‑accepting chemotaxis protein PctB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,566,9470C. 9.9% 220.0 / 18.1 71coding (5/1932 nt)pctBmethyl‑accepting chemotaxis protein PctB
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base . (4/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGC  >  NC_002947/2566801‑2567094
                                                                                                                                                  |                                                                                                                                                     
tGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTc                                                                                                                                                     <  2:193177/149‑1 (MQ=255)
             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)
                                                                                                                                                  |                                                                                                                                                     
TGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTGCTGC  >  NC_002947/2566801‑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: