Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,162,208 G→C 7.4% A326P (GCT→CCT) ‡ PP_4544 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,162,2080GC 7.4% 257.4 / 9.0 81A326P (GCT→CCT) ‡PP_4544hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (3/3);  total (42/39)
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.52e-01

GCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTCGGTTACG  >  NC_002947/5162060‑5162355
                                                                                                                                                     |                                                                                                                                                    
gCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctgg                                                                                                                                                      >  4:200237/1‑149 (MQ=255)
 cAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGTTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctggc                                                                                                                                                     <  8:101254/149‑1 (MQ=255)
 cAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctgg                                                                                                                                                      <  3:326686/148‑1 (MQ=255)
  aGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctggct                                                                                                                                                    >  3:220492/1‑149 (MQ=255)
    tCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctggctgg                                                                                                                                                  <  4:293316/149‑1 (MQ=255)
      gTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGcc                                                                                                                                                <  8:362479/149‑1 (MQ=255)
         cGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGGCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTg                                                                                                                                              <  2:156331/148‑1 (MQ=255)
          ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGa                                                                                                                                             >  4:299969/1‑148 (MQ=255)
          ggCAACACCACCAACGAGTTCAATCTGGTGGAAGCCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGATGGCCTGa                                                                                                                                             >  3:194554/1‑148 (MQ=255)
           gCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaaa                                                                                                                                           <  7:398231/149‑1 (MQ=255)
           gCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaa                                                                                                                                            >  6:36394/1‑148 (MQ=255)
           gCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaa                                                                                                                                            >  8:192757/1‑148 (MQ=255)
                accaccAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCa                                                                                                                                       >  4:419692/1‑148 (MQ=255)
                accaccAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCa                                                                                                                                       >  8:212613/1‑148 (MQ=255)
                  caccaACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCaaaa                                                                                                                                    >  1:85994/1‑149 (MQ=255)
                  caccaACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCAGCCAGCCAGGCCTGatcttcg                                                                                                                                       >  7:84604/1‑140 (MQ=255)
                  caccaACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCAGCCAGCCAGGCCTGatcttcg                                                                                                                                       <  8:84604/146‑7 (MQ=255)
                          aGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcc                                                                                                                            <  7:212613/149‑1 (MQ=255)
                           gTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGTTCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcct                                                                                                                           >  5:373887/1‑149 (MQ=255)
                                   tggtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCTTTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                                   <  6:415971/149‑1 (MQ=255)
                                   tggtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                                   >  2:192265/1‑149 (MQ=255)
                                     gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGATGGCCTGAAAACCAAAATGGGTgcctggagcctgga                                                                                                                 <  4:194554/149‑1 (MQ=255)
                                     gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACATGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctgga                                                                                                                 <  2:4760/149‑1 (MQ=255)
                                     gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTAGGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcct                                                                                                                    <  5:419224/146‑1 (MQ=255)
                                      tggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTGTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctgga                                                                                                                 >  1:90706/1‑149 (MQ=255)
                                        gAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTa                                                                                                             >  2:429566/1‑149 (MQ=255)
                                                   ccAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgc                                                                                                   >  7:108822/1‑148 (MQ=255)
                                                   ccAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgc                                                                                                   >  6:314470/1‑148 (MQ=255)
                                                    cAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGa                                                                                                 >  8:229984/1‑149 (MQ=255)
                                                    cAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGa                                                                                                 >  7:71072/1‑149 (MQ=255)
                                                     aGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgcg                                                                                                  >  2:218667/1‑147 (MQ=255)
                                                     aGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgcg                                                                                                  <  1:218667/147‑1 (MQ=255)
                                                     aGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                                >  1:41489/1‑149 (MQ=255)
                                                     aGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                                <  4:220492/149‑1 (MQ=255)
                                                       aTCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAtgt                                                                                              <  4:37803/149‑1 (MQ=255)
                                                        tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAtgtg                                                                                             >  7:136961/1‑149 (MQ=255)
                                                        tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAtgtg                                                                                             <  8:108822/149‑1 (MQ=255)
                                                        tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATAGGTGCCTGGAGCCTGGACTACAACTACCGCGAtgtg                                                                                             <  1:413648/149‑1 (MQ=255)
                                                          gACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCa                                                                                           <  8:71072/149‑1 (MQ=255)
                                                           aCTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAg                                                                                          >  6:240469/1‑149 (MQ=255)
                                                                    gCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAAcgcg                                                                                 >  1:202721/1‑149 (MQ=255)
                                                                    gCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAAcg                                                                                   <  2:85994/147‑1 (MQ=255)
                                                                        ggCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtg                                                                              <  3:200237/148‑1 (MQ=255)
                                                                         gCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtggt                                                                            <  2:41283/149‑1 (MQ=255)
                                                                         gCCATTCCGTTGTCGGCCTATGGCCAGTACCTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtggt                                                                            >  1:41283/1‑149 (MQ=255)
                                                                               ccGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCAGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGcc                                                                       <  5:36394/148‑1 (MQ=255)
                                                                                cGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGAGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTg                                                                     >  3:135815/1‑149 (MQ=255)
                                                                                  ttGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGtt                                                                   <  6:373887/149‑1 (MQ=255)
                                                                                     gTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTc                                                                  <  3:49811/148‑1 (MQ=255)
                                                                                      tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                                <  5:150941/149‑1 (MQ=255)
                                                                                         gCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtgg                                                                             <  7:399702/133‑1 (MQ=255)
                                                                                         gCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtgg                                                                             >  8:399702/1‑133 (MQ=255)
                                                                                             caagCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                         <  8:192989/146‑1 (MQ=255)
                                                                                             aTGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                         >  1:55748/1‑149 (MQ=255)
                                                                                              tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                         <  2:274359/148‑1 (MQ=255)
                                                                                                gCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGa                                                       <  7:229984/148‑1 (MQ=255)
                                                                                                gCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGAc                                                      <  2:90706/149‑1 (MQ=255)
                                                                                                   aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                                   >  2:8976/1‑149 (MQ=255)
                                                                                                    gTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                         <  4:225238/142‑1 (MQ=255)
                                                                                                    gTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                         >  3:225238/1‑142 (MQ=255)
                                                                                                        gTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                                >  4:312760/1‑131 (MQ=255)
                                                                                                        gTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                                <  3:312760/131‑1 (MQ=255)
                                                                                                         tGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAAc                                             >  2:259659/1‑149 (MQ=255)
                                                                                                          gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACg                                            >  1:118934/1‑149 (MQ=255)
                                                                                                          gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACg                                            >  4:234097/1‑149 (MQ=255)
                                                                                                               acacCGACAGCGCCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTac                                       <  2:123847/149‑1 (MQ=255)
                                                                                                                 acCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTacac                                     >  5:434779/1‑149 (MQ=255)
                                                                                                                  ccGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACAcc                                    >  7:27552/1‑149 (MQ=255)
                                                                                                                      cAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGtt                                <  3:419692/149‑1 (MQ=255)
                                                                                                                          aCCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCg                                   >  6:300658/1‑142 (MQ=255)
                                                                                                                          aCCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCg                                   <  5:300658/142‑1 (MQ=255)
                                                                                                                           ccGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTcgcgc                           <  2:332702/149‑1 (MQ=255)
                                                                                                                                ggCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTAACAAATACCGCGATTTGCAGCGAAACGCGTTGGTCAGCGTGTTCACCCATTCGGAATTCGGCAACGGCTAAACCGGTTCGCGCGGt                        >  1:124891/1‑148 (MQ=255)
                                                                                                                                 gCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTcac                     <  7:192757/149‑1 (MQ=255)
                                                                                                                                   gAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTcac                     >  3:133460/1‑147 (MQ=255)
                                                                                                                                   gAAGATCAGGCCTGACTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTcac                     >  4:140973/1‑147 (MQ=255)
                                                                                                                                      tttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCg                                                  <  8:178718/113‑1 (MQ=255)
                                                                                                                                      tttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCg                                                  >  7:178718/3‑115 (MQ=255)
                                                                                                                                      tttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAAc                                             <  2:121931/118‑1 (MQ=255)
                                                                                                                                      tttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAAc                                             >  1:121931/3‑120 (MQ=255)
                                                                                                                                        tCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCgg                                                        <  1:352430/107‑1 (MQ=255)
                                                                                                                                        tCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCgg                                                        >  2:352430/1‑107 (MQ=255)
                                                                                                                                          aGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAg            <  6:308167/149‑1 (MQ=255)
                                                                                                                                            gCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGt          >  6:373656/1‑149 (MQ=255)
                                                                                                                                            gCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGt          >  8:53044/1‑149 (MQ=255)
                                                                                                                                             ccTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTc         <  8:27552/149‑1 (MQ=255)
                                                                                                                                              ctggctggctggCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTCg        <  5:314470/149‑1 (MQ=255)
                                                                                                                                               tggctggctggcCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTCgg       >  1:98717/1‑149 (MQ=255)
                                                                                                                                                     gctggcCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTCGGTTAc   <  2:55748/147‑1 (MQ=255)
                                                                                                                                                     gctggcCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTCGGTTACg  <  5:240469/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGTCGGTTACG  >  NC_002947/5162060‑5162355

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: