Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,566,942 G→C 7.0% L4V (CTG→GTG)  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,9420GC 7.0% 227.0 / 5.3 71L4V (CTG→GTG) pctBmethyl‑accepting chemotaxis protein PctB
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (3/2);  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 = 5.40e-01

TGGCGGATGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTG  >  NC_002947/2566794‑2567090
                                                                                                                                                    |                                                                                                                                                     
tGGCGGATGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAg                                                                                                                                                       >  1:196246/1‑149 (MQ=255)
tGGCGGATGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAg                                                                                                                                                       >  2:107109/1‑149 (MQ=255)
  gCGGATGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAgcg                                                                                                                                                     >  2:194672/1‑149 (MQ=255)
       tGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTc                                                                                                                                                <  2:193177/149‑1 (MQ=255)
                    cAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCg                                                                                                                                  >  7:152587/1‑149 (MQ=255)
                    cAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCACCGCGTCATGAGCGCGCTGc                                                                                                                                   <  1:371506/148‑1 (MQ=255)
                      ggTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCggg                                                                                                                                <  2:196246/149‑1 (MQ=255)
                      ggTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCggg                                                                                                                                <  1:235304/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                                                                                                             >  5:355303/1‑149 (MQ=255)
                                         gAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATc                                                                                                             >  2:115857/1‑149 (MQ=255)
                                          aaCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCg                                                                                                            >  1:173255/1‑149 (MQ=255)
                                           aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCgg                                                                                                           >  5:251375/1‑149 (MQ=255)
                                           aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCgg                                                                                                           >  3:141244/1‑149 (MQ=255)
                                           aCAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCgg                                                                                                           <  8:344189/149‑1 (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                                                                                               >  8:234151/1‑148 (MQ=255)
                                                        cAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTaa                                                                                               >  6:304501/1‑148 (MQ=255)
                                                         aaaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCga                                                                                                                   <  2:205425/127‑1 (MQ=255)
                                                         aaaaCCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCga                                                                                                                   >  1:205425/1‑127 (MQ=255)
                                                              cACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGTCATGa                                                                                                                                            >  3:346212/1‑97 (MQ=255)
                                                              cACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGTCATGa                                                                                                                                            <  4:346212/97‑1 (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:107109/149‑1 (MQ=255)
                                                                        cAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAg                                                                              <  1:194672/149‑1 (MQ=255)
                                                                             gcTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGaa                                                                          >  7:221111/1‑148 (MQ=255)
                                                                               tGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGaa                                                                                      <  6:129141/134‑1 (MQ=255)
                                                                               tGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGaa                                                                                      >  5:129141/1‑134 (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                                                                                         >  4:125830/1‑96 (MQ=255)
                                                                                                                  tttCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGAtt                                                                                         <  3:125830/96‑1 (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)
                                                                                                                                                    |                                                                                                                                                     
TGGCGGATGGTATTGCGTTGCAGGTAGTCGTTGTAGAGCGTGAACAGGGCAAATGCCAAAACCACGACGCCGCAGGCGCTGAGCAGGATCTTGTGGCGAAATTTCAGGTTCATGTTTCGAGACCTTGAGCCAATGAGGGGAGGGCGCAGCGCGCTCATGACGCGGTGCGGGGCCTATCGACCGAGATATCGGCTTGAGGGGTAAAAGATTGAAGTGCTTAGGGAAAATTCTTACAGACGTTTCGACGGACGGCAGTCACGGATCAGCCAGCTGCTCGATCAGCAGGCGTGCAGCCTG  >  NC_002947/2566794‑2567090

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: