Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,955,263 C→T 5.9% A38T (GCC→ACC)  fliJ ← flagellar protein FliJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,955,2630CT 5.9% 208.4 / 2.2 68A38T (GCC→ACC) fliJflagellar protein FliJ
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GGCCGTTTCCAGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCGTGAAAATCGCCC  >  NC_002947/4955129‑4955409
                                                                                                                                       |                                                                                                                                                  
ggCCGTTTCCAGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAACGCGCTCCAGCTCGGCCAGCTTGGCCGGGGCCGTGGCCAgctgct                                                                                                                                      >  4:194264/1‑149 (MQ=255)
   cGTTTCCAGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAActgctgct                                                                                                                                   <  8:160898/149‑1 (MQ=255)
     tttCCAGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAActgctgctg                                                                                                                                  >  1:55977/1‑148 (MQ=255)
         cAGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAAt                                                                                                                             <  5:309965/149‑1 (MQ=255)
          aGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATg                                                                                                                            <  5:195996/149‑1 (MQ=255)
          aGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATg                                                                                                                            >  7:60045/1‑149 (MQ=255)
           gCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATgg                                                                                                                           <  5:156356/149‑1 (MQ=255)
               ccccAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGcc                                                                                                                         >  6:200593/1‑147 (MQ=255)
                          gcTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAgcc                                                                                                            >  7:161986/1‑149 (MQ=255)
                            tGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAgccgc                                                                                                          >  4:149517/1‑149 (MQ=255)
                                 aTTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCATATGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTg                                                                                                      >  1:232937/1‑149 (MQ=255)
                                 aTTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTTCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTTGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCATTgc                                                                                                     >  7:276112/1‑149 (MQ=255)
                                 aTTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTg                                                                                                      <  3:149517/148‑1 (MQ=255)
                                     aCCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTc                                                                                                 <  8:161986/149‑1 (MQ=255)
                                      ccagccagCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAgccgc                                                                                                          >  7:299720/1‑139 (MQ=255)
                                      ccagccagCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAgccgc                                                                                                          <  8:299720/139‑1 (MQ=255)
                                      ccagccagCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCg                                                                                                >  5:103922/1‑149 (MQ=255)
                                      ccagccagCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCg                                                                                                <  1:38951/149‑1 (MQ=255)
                                        agccagcTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGc                                                                                              >  4:215970/1‑149 (MQ=255)
                                         gccagcGTCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGcc                                                                                             <  3:194264/149‑1 (MQ=255)
                                                 ccATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGc                                                                                     <  6:103922/149‑1 (MQ=255)
                                                  cATTGACCCCTTGCCCGCCGCGGGTGATCCACTGCAACTGGTTATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGcc                                                                                    <  6:36378/149‑1 (MQ=255)
                                                   aTTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGTCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCa                                                                                   >  5:136538/1‑149 (MQ=255)
                                                   aTTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGTCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCa                                                                                   >  5:131737/1‑149 (MQ=255)
                                                       aCCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTc                                                                               >  3:221652/1‑149 (MQ=255)
                                                           cTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGa                                                                                                               <  5:226216/113‑1 (MQ=255)
                                                           cTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTccacc                                                                           <  3:85619/149‑1 (MQ=255)
                                                                  gccgcGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGgcg                                                                    >  5:165604/1‑149 (MQ=255)
                                                                            aTCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGgcgcgc                                                          >  5:236041/1‑149 (MQ=255)
                                                                               cACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTCGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCgg                                                        >  6:19904/1‑148 (MQ=255)
                                                                                aCTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGCCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCa                                                      >  6:88840/1‑149 (MQ=255)
                                                                                  cTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCAc                                                     >  4:259697/1‑149 (MQ=255)
                                                                                   tGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACg                                                    >  3:271348/1‑149 (MQ=255)
                                                                                   tGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACg                                                    >  4:299423/1‑149 (MQ=255)
                                                                                       aCTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCg                                                <  7:250996/149‑1 (MQ=255)
                                                                                       aCTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCg                                                <  2:232937/149‑1 (MQ=255)
                                                                                        cTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCgg                                               >  2:302762/1‑149 (MQ=255)
                                                                                           gTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCggg                                              >  6:164800/1‑147 (MQ=255)
                                                                                           gTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCggg                                              >  4:203061/1‑147 (MQ=255)
                                                                                             aaTCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCgg                                               <  6:321952/144‑1 (MQ=255)
                                                                                             aaTCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCgg                                               >  5:321952/1‑144 (MQ=255)
                                                                                              aTCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAgc                                          <  3:215970/148‑1 (MQ=255)
                                                                                                ccTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAAc                                                                                                                                          >  3:120902/1‑50 (MQ=255)
                                                                                                ccTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAAc                                                                                                                                          <  4:120902/50‑1 (MQ=255)
                                                                                                 cTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCggg                                              >  8:294298/1‑141 (MQ=255)
                                                                                                 cTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCggg                                              <  7:294298/141‑1 (MQ=255)
                                                                                                  tCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCAt                                     >  8:298791/1‑149 (MQ=255)
                                                                                                    gcgAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATgg                                   >  2:58627/1‑149 (MQ=255)
                                                                                                    gcgAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATgg                                   <  5:88840/149‑1 (MQ=255)
                                                                                                    gcgAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATgg                                   <  2:278058/149‑1 (MQ=255)
                                                                                                    gcgAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATg                                    <  7:298791/148‑1 (MQ=255)
                                                                                                      gAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCt                                 <  6:236041/149‑1 (MQ=255)
                                                                                                       aaaGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCtt                                <  8:60045/149‑1 (MQ=255)
                                                                                                         aGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTAt                              >  5:95519/1‑149 (MQ=255)
                                                                                                         aGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTAt                              <  8:30667/149‑1 (MQ=255)
                                                                                                             cTCCAGCTCGGCCAGCTTGGCCTGGGTCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAAc                          <  6:136538/149‑1 (MQ=255)
                                                                                                             cTCCAGCTCGGCCAGCTTGGCCTGGGTCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAAc                          <  6:131737/149‑1 (MQ=255)
                                                                                                             cTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAAc                          >  8:188688/1‑149 (MQ=255)
                                                                                                              tCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAAcc                         >  7:97779/1‑149 (MQ=255)
                                                                                                                   cTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCg                    <  5:19904/149‑1 (MQ=255)
                                                                                                                      ggCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCgggg                 <  4:221652/149‑1 (MQ=255)
                                                                                                                      ggCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCgggg                 >  2:121548/1‑149 (MQ=255)
                                                                                                                      ggCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCgggg                 <  1:302762/149‑1 (MQ=255)
                                                                                                                         cAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCg              >  8:131252/1‑149 (MQ=255)
                                                                                                                                gCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCGTGAAAAt       <  5:164800/149‑1 (MQ=255)
                                                                                                                                 ccTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCGTGAAAATc      <  3:299423/149‑1 (MQ=255)
                                                                                                                                     ggCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCGTGAATATCGc    <  6:165604/147‑1 (MQ=255)
                                                                                                                                     ggCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCGTGAAAATCGccc  <  7:131252/149‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                  
GGCCGTTTCCAGCTGCCCCAGAAAGCGCTGGTAATTGACCAGCCAGCTGCCATTGACCCCTTGCCCGCCGCGGTTGATCCACTGCAACTGGTAATCCTCGCGAAAGCGCTCCAGCTCGGCCAGCTTGGCCTGGGCCGTGGCCAACTGCTGCTGGAAATGGCCAAGGCGCTGAGCCGCCTTGCGCTCGGCCTCTTCGGCCATGTCCACCACCGGCGCCAGGCGCGCGGCACGTCCGGGCAGCGCCATGGCTTATCAACCTGCCGGGGGCGTGAAAATCGCCC  >  NC_002947/4955129‑4955409

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: