Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,872,133 T→C 10.3% Q83R (CAA→CGA)  PP_3420 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,872,1330TC10.3% 181.2 / 9.6 68Q83R (CAA→CGA) PP_3420sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (4/3);  total (35/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.47e-01

ATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATC  >  NC_002947/3871988‑3872279
                                                                                                                                                 |                                                                                                                                                  
aTTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGc                                                                                                                                                  >  3:303980/1‑148 (MQ=255)
aTTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGc                                                                                                                                                  >  8:272434/1‑148 (MQ=255)
     ccAGCAAGGTGTCCGGCAGGTGCTCGACATTATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGACGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGc                                                                                                                                            <  7:17479/149‑1 (MQ=255)
      cAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGcc                                                                                                                                           <  7:267641/149‑1 (MQ=255)
          aaGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCATTCGCAATTGCCGcct                                                                                                                                       >  4:175730/1‑147 (MQ=255)
           aGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCATTCGCAATTGCCGcctg                                                                                                                                      <  3:175730/149‑4 (MQ=255)
                            tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGc                                                                                                                           <  3:174187/143‑1 (MQ=255)
                            tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGc                                                                                                                           >  4:174187/1‑143 (MQ=255)
                            tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                     >  4:39432/1‑149 (MQ=255)
                            tCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                     >  8:195911/1‑149 (MQ=255)
                              gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                     <  3:39432/147‑1 (MQ=255)
                              gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                     <  8:32910/147‑1 (MQ=255)
                              gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGc                                                                                                                     <  1:183283/147‑1 (MQ=255)
                              gACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCt                                                                                                                   <  7:272434/149‑1 (MQ=255)
                                 aTGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCtgaa                                                                                                                <  1:190455/149‑1 (MQ=255)
                                 aTGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCtgaa                                                                                                                <  1:147619/149‑1 (MQ=255)
                                  tGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCtgaat                                                                                                               <  8:86609/149‑1 (MQ=255)
                                              tCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTg                                                                                                         >  2:75591/1‑143 (MQ=255)
                                              tCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTg                                                                                                         <  1:75591/143‑1 (MQ=255)
                                                gATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGt                                                                                                   <  4:303980/147‑1 (MQ=255)
                                                gATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGt                                                                                                   >  4:124108/1‑147 (MQ=255)
                                                    cGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCAt                                                                                             >  1:80867/1‑149 (MQ=255)
                                                     gggTCGATTCGCCCACCGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCCCCGCTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGCCCATTCCCTGAATGAACTGCGCGGTAAGCATg                                                                                            >  4:165941/1‑149 (MQ=255)
                                                               ggCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGACAGGCGGCAATTGCGAATGCCGCACTTCGGACTGGGACATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGctct                                                                                  >  8:177139/1‑149 (MQ=255)
                                                                gCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCAtt                                                                                                                       <  6:28298/111‑1 (MQ=255)
                                                                gCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCAtt                                                                                                                       >  5:28298/1‑111 (MQ=255)
                                                                 ccAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTcc                                                                                <  2:245240/149‑1 (MQ=255)
                                                                      gCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCGGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCt                                                                           <  7:195911/149‑1 (MQ=255)
                                                                              gACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGca                                                                   >  2:132843/1‑149 (MQ=255)
                                                                                 cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCgcgc                                                                    >  3:353205/1‑145 (MQ=255)
                                                                                 cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCgcgc                                                                    <  4:353205/145‑1 (MQ=255)
                                                                                 cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcccca                                                                >  6:341414/1‑149 (MQ=255)
                                                                                 cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                >  8:190083/1‑149 (MQ=255)
                                                                                 cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                >  2:377430/1‑149 (MQ=255)
                                                                                 cTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacca                                                                >  7:98083/1‑149 (MQ=255)
                                                                                    cAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGGGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACg                                                              >  1:133416/1‑148 (MQ=255)
                                                                                      aGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAg                                                           <  8:47239/149‑1 (MQ=255)
                                                                                      aGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAg                                                           <  7:177139/149‑1 (MQ=255)
                                                                                       gTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGcacacagtat                                                           >  3:271481/1‑141 (MQ=255)
                                                                                            gcAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCTCGCTGTCGCGCACCACGAATCGCGTc                                                     >  8:173358/1‑149 (MQ=255)
                                                                                                 tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCGAGGCTGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATTAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                <  6:237534/149‑1 (MQ=255)
                                                                                                 tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                <  1:132843/149‑1 (MQ=255)
                                                                                                 tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                <  1:188373/149‑1 (MQ=255)
                                                                                                 tCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAAc                                                <  7:190083/149‑1 (MQ=255)
                                                                                                        gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAgcg                                                         >  4:208623/1‑133 (MQ=255)
                                                                                                        gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAgcg                                                         <  3:208623/133‑1 (MQ=255)
                                                                                                        gcgGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAAcc                                          >  3:237006/1‑148 (MQ=255)
                                                                                                                aGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCAcc                                 >  3:110921/1‑149 (MQ=255)
                                                                                                                aGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCAcc                                 >  5:329130/1‑149 (MQ=255)
                                                                                                                    aCTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAg                                       <  3:33983/139‑1 (MQ=255)
                                                                                                                    aCTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAg                                       >  4:33983/1‑139 (MQ=255)
                                                                                                                         cccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCAcc                          >  1:25985/1‑147 (MQ=255)
                                                                                                                          ccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGt                                                      <  6:86536/118‑1 (MQ=255)
                                                                                                                          ccATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGt                                                      >  5:86536/1‑118 (MQ=255)
                                                                                                                               gtggtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTaa                   <  4:345452/148‑1 (MQ=255)
                                                                                                                                tggtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATg                 <  4:237006/149‑1 (MQ=255)
                                                                                                                                  gtgCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGAt               <  4:110921/149‑1 (MQ=255)
                                                                                                                                   cgaagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGAt               >  4:82833/7‑148 (MQ=255)
                                                                                                                                     ccAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa            <  2:25985/149‑1 (MQ=255)
                                                                                                                                     aagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGa                <  1:294575/141‑1 (MQ=255)
                                                                                                                                     aagtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGGCGAAACGTAACCCAGCCCCCCTGCCACCGCGGTAATGa                >  2:294575/5‑145 (MQ=255)
                                                                                                                                       gtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAgg                      <  3:193776/135‑1 (MQ=255)
                                                                                                                                       gtgcggcATTCGCAATTGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAgg                      >  4:193776/3‑137 (MQ=255)
                                                                                                                                       aGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa            >  8:241557/1‑147 (MQ=255)
                                                                                                                                       aGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCaa            >  5:92026/1‑147 (MQ=255)
                                                                                                                                         gcgggAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGAGCGACGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAgg        >  7:158551/1‑149 (MQ=255)
                                                                                                                                              aaTTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGAATGGCATCGCGCTCGACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGAt   >  5:259243/1‑149 (MQ=255)
                                                                                                                                               aTTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATc  <  3:124108/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATTGGCCAGCAAGGTGTCCGGCAGGTGCTCGACATGATCGAGGAACTCGATCCGGGTCGATTCGGCCAGGGCCGGGGTGACCTGCAAGTGCTGCAGGTCCAGCCGCGGGTCCAGCAACTCGCCCATGGTGGTGCCAGGCGGCAATTGCGAATGCCGCACTTCGGCCTGGGCCATTGCCTGAATGAACTGCGCGGTAAGCATGGCATCGCGCTCCACGCTGTCGCGCACCACGAAGCGCGTCGAAACGTAACCCAGCCCACCTGCCACCGAGGTAATGATCAACAGGCTGATC  >  NC_002947/3871988‑3872279

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