Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,285,376 A→T 5.5% E594V (GAG→GTG)  PP_4658 → methyl‑accepting chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,285,3760AT 5.5% 219.4 / 2.1 72E594V (GAG→GTG) PP_4658methyl‑accepting chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base T (2/2);  total (38/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 = 6.50e-01

GGCAGCCGTGCCCTGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCA  >  NC_002947/5285228‑5285524
                                                                                                                                                    |                                                                                                                                                    
ggCAGCCGTGCCCTGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAaga                                                                                                                                                      <  6:95801/149‑1 (MQ=255)
 gCAGCCGTGCCCTGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAagag                                                                                                                                                     <  7:352416/149‑1 (MQ=255)
  cacggctGCCCTGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAgaga                                                                                                                                                    <  7:293177/143‑1 (MQ=255)
           ccTGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCg                                                                                                                                             >  4:152600/1‑147 (MQ=255)
            cTGACCGACAGCAGGGTTGATCTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTAGGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTTCGGCGGCAGAGCAGCAGAGTGCGGTGTCTGAAGAGATCAGCCGAAg                                                                                                                                          <  2:297864/149‑1 (MQ=38)
             tGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAg                                                                                                                                          >  2:121042/1‑148 (MQ=255)
                 cGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATcc                                                                                                                                     >  3:197243/1‑149 (MQ=255)
                  gACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCt                                                                                                                                    <  8:205773/149‑1 (MQ=255)
                      gcagcGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCCGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAAt                                                                                                                                <  7:219960/149‑1 (MQ=255)
                        agcGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAAtgt                                                                                                                              <  2:241858/149‑1 (MQ=255)
                              gAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGTAGAGCAGCAGAGTGCGGTGGCTGAAGTGATCAGCCGAAGCATCCTGAATGTGCGTGa                                                                                                                        >  2:271136/1‑149 (MQ=255)
                              gAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGAAGAGCAGCAGAGTGCGGTGGCTGAAGTGATCAGCCGAAGCATCCTGAATGTGCGTGa                                                                                                                        >  2:265435/1‑149 (MQ=255)
                                    gcgcGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATgtgt                                                                                                                   <  4:335881/148‑1 (MQ=255)
                                       cgAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCaga                                                                                                               <  6:34970/149‑1 (MQ=255)
                                       cgAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCaga                                                                                                               <  8:140905/149‑1 (MQ=255)
                                         aaaGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCagag                                                                                                              >  4:338909/1‑148 (MQ=255)
                                         aaaGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATAGAGCATCCAGTCGATGAACCAGCAGATTGCTGAGGCGGCAGAGCAGAAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCagag                                                                                                              >  7:353169/1‑148 (MQ=255)
                                            ggCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCaaa                                                                                                          <  6:65542/149‑1 (MQ=255)
                                             gCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAAc                                                                                                         >  7:267107/1‑149 (MQ=255)
                                                gggTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCg                                                                                                       >  2:149391/1‑148 (MQ=255)
                                                 ggTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCgca                                                                                                     <  5:154188/149‑1 (MQ=255)
                                                 ggTAATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCTATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATGAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCgca                                                                                                     <  8:353169/149‑1 (MQ=255)
                                                  gTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCgcag                                                                                                    >  5:102829/1‑149 (MQ=255)
                                                     aTCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGcagcag                                                                                                 >  5:202713/1‑149 (MQ=255)
                                                      tCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCagcagc                                                                                                >  8:2553/1‑149 (MQ=255)
                                                          tGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCCGCCGCAGCCTCCTGGATGTGCGTGATGTGGCAGCGCACGCCGCAGCAGCCAGt                                                                                            >  6:74745/1‑149 (MQ=37)
                                                           ggAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGt                                                                                            >  7:37569/1‑148 (MQ=255)
                                                            gagaGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGa                                                                                          >  3:274550/1‑149 (MQ=255)
                                                             agagCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGAc                                                                                         <  2:88151/149‑1 (MQ=255)
                                                             agagCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGAc                                                                                         >  1:271617/1‑149 (MQ=255)
                                                                 cATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGCGATCAGCCGAAGCATCCTGAATGTGCGTGATCTGTCAGAGCAAACCGCAGCAGCCAGTGACgaga                                                                                     >  1:346662/1‑149 (MQ=255)
                                                                 cATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATgtgt                                                                                                                   <  8:274156/119‑1 (MQ=255)
                                                                 cATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATgtgt                                                                                                                   >  7:274156/1‑119 (MQ=255)
                                                                 cATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACgaga                                                                                     >  5:32806/1‑149 (MQ=255)
                                                                 cATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACgaga                                                                                     >  3:168616/1‑149 (MQ=255)
                                                                 cATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACgaga                                                                                     >  3:159327/1‑149 (MQ=255)
                                                                  aTTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAACCCGCAGCAGCCAGTGACGAGAc                                                                                    >  6:234578/1‑149 (MQ=255)
                                                                    tACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTc                                                                                                                  >  1:365802/1‑117 (MQ=255)
                                                                    tACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTc                                                                                                                  <  2:365802/117‑1 (MQ=255)
                                                                    tACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACgg                                                                                  >  2:272927/1‑149 (MQ=255)
                                                                     aCCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGAAGTGTCAGCGCAAACCGCAGCAGCCAGTGACGAGACgg                                                                                  >  3:361752/1‑148 (MQ=255)
                                                                      ccAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCa                                                                                <  3:152600/149‑1 (MQ=255)
                                                                       cAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGAt                                                                                                                       <  5:317270/109‑1 (MQ=255)
                                                                       cAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGAt                                                                                                                       >  6:317270/1‑109 (MQ=255)
                                                                        aGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGc                                                                              >  6:256579/1‑149 (MQ=255)
                                                                             ggTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGa                                                                         <  2:271617/149‑1 (MQ=255)
                                                                                aTCGAGCATCCAGTCGATGAACCAGCAGATTGCTGAGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAgt                                                                       <  8:267107/148‑1 (MQ=255)
                                                                                 tCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAgtgt                                                                     <  1:149391/149‑1 (MQ=255)
                                                                                 tCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAgtg                                                                      <  4:197243/148‑1 (MQ=255)
                                                                                   gAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGtgtg                                                                    <  4:159327/148‑1 (MQ=255)
                                                                                          cAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGgc                                                            >  2:306331/1‑149 (MQ=255)
                                                                                           aGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGgcg                                                           >  2:363267/1‑149 (MQ=255)
                                                                                            gTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGgcgc                                                          <  8:37569/149‑1 (MQ=255)
                                                                                                aTGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGtt                                                      <  1:363267/149‑1 (MQ=255)
                                                                                                   aaCCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTtgg                                                    <  3:338909/148‑1 (MQ=255)
                                                                                                       agcagATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTtgggtggg                                               >  7:363930/1‑149 (MQ=255)
                                                                                                        gcagATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGa                                                                  >  6:114943/1‑129 (MQ=255)
                                                                                                        gcagATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGa                                                                  <  5:114943/129‑1 (MQ=255)
                                                                                                        gcagATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTtgggtggg                                               >  7:270656/1‑148 (MQ=255)
                                                                                                        gcagATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTtgggtggg                                               >  5:76649/1‑148 (MQ=255)
                                                                                                            aTTGCTGCGGCGGTAGAGCAGCAGAGTGCGGTGGCTGAAGTGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCTTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGtt                                          <  1:271136/149‑1 (MQ=255)
                                                                                                            aTTGCTGCGGCGGTAGAGCAGCAGAGTGCGGTGGCTGAAGTGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGtt                                          <  1:265435/149‑1 (MQ=255)
                                                                                                             ttGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTg                                         >  6:352496/1‑149 (MQ=255)
                                                                                                               gctgcGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGc                                        <  6:102829/148‑1 (MQ=255)
                                                                                                                 tgcGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGa                                     <  6:76649/149‑1 (MQ=255)
                                                                                                                   cggcggCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGa                                     <  6:202713/147‑1 (MQ=255)
                                                                                                                    agcggcAGAGCAGCAGAGGGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGGGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGGGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGacac                                  <  5:283232/148‑1 (MQ=255)
                                                                                                                       ggcAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTgg                               >  6:33625/1‑149 (MQ=255)
                                                                                                                         cAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTc                             <  1:306331/149‑1 (MQ=255)
                                                                                                                            agcagcAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGc                          <  1:272927/149‑1 (MQ=255)
                                                                                                                                 cagAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCCAGtt                     >  7:169219/1‑149 (MQ=255)
                                                                                                                                        cGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCAAGTTCCGCGTc              >  7:373477/1‑149 (MQ=255)
                                                                                                                                         ggTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCt             <  4:274550/149‑1 (MQ=255)
                                                                                                                                          gTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTCGGTGGGCAGTTGCAGACACTGGTCAGCCAGATCCGCGTCt             >  3:42684/1‑148 (MQ=255)
                                                                                                                                                 aaGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTg     <  6:32806/149‑1 (MQ=255)
                                                                                                                                                    agaTCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCa  <  5:256579/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCAGCCGTGCCCTGACCGACAGCAGCGTCGAACTGGCGCGAAAGGCCGGGTCATCGCTGGAGAGCATTACCAGTACGGTATCGAGCATCCAGTCGATGAACCAGCAGATTGCTGCGGCGGCAGAGCAGCAGAGTGCGGTGGCTGAAGAGATCAGCCGAAGCATCCTGAATGTGCGTGATGTGTCAGAGCAAACCGCAGCAGCCAGTGACGAGACGGCAGCGTCGAGTGTGGAGCTGGCGCGGTTGGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCA  >  NC_002947/5285228‑5285524

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