Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 381,856 T→A 10.9% D230V (GAC→GTC)  PP_0317 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947381,8560TA10.9% 170.1 / 10.1 64D230V (GAC→GTC) PP_0317methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (24/33);  new base A (4/3);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

TTCTGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTGGC  >  NC_002947/381711‑382001
                                                                                                                                                 |                                                                                                                                                 
ttctGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCgg                                                                                                                                                >  8:188653/1‑149 (MQ=255)
  ctGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCTAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGc                                                                                                                                               <  6:96748/148‑1 (MQ=255)
     ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCg                                                                                                                                              <  3:264426/146‑1 (MQ=255)
      cgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACtt                                                                                                                                          <  1:47219/149‑1 (MQ=255)
       gccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTg                                                                                                                                         <  3:122198/149‑1 (MQ=255)
       gccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTg                                                                                                                                         >  5:277607/1‑149 (MQ=255)
             gCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGgc                                                                                                                                   <  8:142241/149‑1 (MQ=255)
                aGCTGGGGGGCGGAGCCGGCGATCTGCTGCAAGGTAGCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCa                                                                                                                                <  7:21123/149‑1 (MQ=255)
                     ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                             <  1:282069/147‑1 (MQ=255)
                     ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCAttt                                                                                                                             <  7:188653/147‑1 (MQ=255)
                     ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGa                                                                                                                           <  6:163185/149‑1 (MQ=255)
                       tGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGt                                                                                                                         <  6:228016/149‑1 (MQ=255)
                        ggCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTg                                                                                                                        >  2:105204/1‑149 (MQ=255)
                         gCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGc                                                                                                                       <  2:209160/149‑1 (MQ=255)
                         gCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGc                                                                                                                       <  6:91180/149‑1 (MQ=255)
                         gCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGc                                                                                                                       <  6:95502/149‑1 (MQ=255)
                         gCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGc                                                                                                                       <  4:71134/149‑1 (MQ=255)
                              gCCGGCGATCTGCTCCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGcttctt                                                                                                                  <  4:227016/149‑1 (MQ=255)
                               ccGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCttcttc                                                                                                                 <  7:156583/149‑1 (MQ=255)
                                cGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCttcttc                                                                                                                 >  1:78718/1‑148 (MQ=255)
                                   cGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATc                                                                                                             <  6:277607/149‑1 (MQ=255)
                                   cGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATc                                                                                                             <  1:119382/149‑1 (MQ=255)
                                           tgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGAc                                                                                                     >  3:206020/1‑149 (MQ=255)
                                             cAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGa                                                                                                   >  8:72755/1‑149 (MQ=255)
                                                ggTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGACCTg                                                                                                                                         >  1:196580/1‑108 (MQ=255)
                                                ggTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGACCTg                                                                                                                                         <  2:196580/108‑1 (MQ=255)
                                                ggTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGAAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATAc                                                                                                >  2:183097/1‑149 (MQ=255)
                                                 gTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACt                                                                                               >  3:269482/1‑149 (MQ=255)
                                                    tCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGc                                                                                             <  1:105204/148‑1 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCggg                                                                                          >  5:138421/1‑149 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGAAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCgg                                                                                           <  1:183097/148‑1 (MQ=255)
                                                                tCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAgg                                                                                  >  3:91232/1‑147 (MQ=255)
                                                                 ggctgttGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAg                                                                               >  3:60073/7‑149 (MQ=255)
                                                                  cTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAg                                                                               <  5:245807/148‑1 (MQ=255)
                                                                     cATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGcac                                                                           <  7:72755/149‑1 (MQ=255)
                                                                     cATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGcac                                                                           <  4:91232/149‑1 (MQ=255)
                                                                      aTGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGca                                                                            >  6:214816/1‑147 (MQ=255)
                                                                        gcgcgcCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGCAGTGGGTCAGGTCGCCGTAGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTGGATCGGCTTGACGATACTCCGGGTCAGCAGGAAGGCACAgg                                                                        >  3:70251/1‑149 (MQ=255)
                                                                          gcgcCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACTGGCg                                                                      >  3:146844/1‑149 (MQ=255)
                                                                              cATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAgg                                                                  >  5:177590/1‑149 (MQ=255)
                                                                              cATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAgg                                                                  >  2:291092/1‑149 (MQ=255)
                                                                                    gcgcAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTc                                                                                        >  8:58157/1‑121 (MQ=255)
                                                                                    gcgcAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTc                                                                                        <  7:58157/121‑1 (MQ=255)
                                                                                        aaCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGt                                                        >  6:106918/1‑149 (MQ=255)
                                                                                        aaCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGAGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGt                                                        >  1:35633/1‑149 (MQ=255)
                                                                                           aGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATCACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTggc                                                     >  2:19091/1‑149 (MQ=255)
                                                                                                 gcggcTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAAc                                               <  2:35633/149‑1 (MQ=255)
                                                                                                 gcggcTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAAc                                               <  6:287659/149‑1 (MQ=255)
                                                                                                    gcTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACcag                                            <  2:78718/149‑1 (MQ=255)
                                                                                                      ttCGTCGGTGCCCTCTGCGCGGATGACGTGTGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGAAGGAAGGCACAGGCGAAGGTAAGTACGGTGGCTGCAACcag                                            >  4:170410/1‑147 (MQ=255)
                                                                                                         gtcgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGc                                       <  4:269482/149‑1 (MQ=255)
                                                                                                          tcgGTGCCGTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAAGCAGCAGGcc                                      <  4:206020/149‑1 (MQ=255)
                                                                                                              tGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATc                                  <  3:82328/149‑1 (MQ=255)
                                                                                                                  ctctGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCAc                                >  6:156268/1‑147 (MQ=255)
                                                                                                                          ggATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAgg                            <  5:235852/143‑1 (MQ=255)
                                                                                                                          ggATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGaaa                      <  5:156268/149‑1 (MQ=255)
                                                                                                                          ggATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGAACGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAgg                            >  6:235852/1‑143 (MQ=255)
                                                                                                                                    gagcaagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGc                          >  3:145189/7‑135 (MQ=255)
                                                                                                                                    gagcaagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGc                          <  4:145189/129‑1 (MQ=255)
                                                                                                                                    gagcaagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCa             >  1:82882/7‑148 (MQ=255)
                                                                                                                                        aagTCGCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTg        >  3:294642/3‑149 (MQ=255)
                                                                                                                                            cGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTg    >  6:271373/1‑149 (MQ=255)
                                                                                                                                             gCCGACGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTg    <  6:138693/148‑1 (MQ=255)
                                                                                                                                              ccGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCTGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTggc  <  5:214816/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCTGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTGGC  >  NC_002947/381711‑382001

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: ‑