Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,607,908 G→A 7.8% L77L (CTC→CTT cyaB ← cyclolysin secretion/processing ATP‑binding protein CyaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,607,9080GA 7.8% 257.3 / 8.7 89L77L (CTC→CTTcyaBcyclolysin secretion/processing ATP‑binding protein CyaB
Reads supporting (aligned to +/- strand):  ref base G (46/36);  new base A (3/4);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CGGCGCGCGGGAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACTTTCAA  >  NC_002947/5607761‑5608053
                                                                                                                                                    |                                                                                                                                                 
cggcgCGCGGGAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGa                                                                                                                                                  <  1:357516/149‑1 (MQ=255)
 ggcgCGCGGGAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAg                                                                                                                                                 >  4:244272/1‑149 (MQ=255)
  gcgcgcGGGAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAgg                                                                                                                                                >  5:83772/1‑149 (MQ=255)
    gcgcgGGAATGTGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCTAGGGCGTCTGCTGCAGCCGGGCCGTATCGACGCTGCGGTGTTTGGCCTTGAGGcc                                                                                                                                              <  1:335241/149‑1 (MQ=255)
        gggAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGt                                                                                                                                          >  5:271577/1‑149 (MQ=255)
         ggAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGtt                                                                                                                                         <  5:249408/149‑1 (MQ=255)
          gAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTg                                                                                                                                        >  4:375576/1‑149 (MQ=255)
          gAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTCAAGGCCaaacac                                                                                                                                        <  2:126915/149‑5 (MQ=255)
          gAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTCAAGGCCaaacac                                                                                                                                        <  5:249402/149‑5 (MQ=255)
           aaTGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGc                                                                                                                                       <  4:25553/149‑1 (MQ=255)
           aaTGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCAGGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTCAAGGCCaaacacc                                                                                                                                       <  1:147590/149‑6 (MQ=255)
                 ggTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAAt                                                                                                                                 <  1:178810/149‑1 (MQ=255)
                   tCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGc                                                                                                                               <  4:182685/149‑1 (MQ=255)
                     cTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGccc                                                                                                                             >  3:199030/1‑149 (MQ=255)
                      tGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGccc                                                                                                                             >  8:244961/1‑148 (MQ=255)
                           cAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAagta                                                                                                                       <  3:244272/149‑1 (MQ=255)
                             gAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGc                                                                                                                                       <  5:391972/131‑1 (MQ=255)
                             gAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGc                                                                                                                                       >  6:391972/1‑131 (MQ=255)
                             gAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAagtagt                                                                                                                     >  7:122299/1‑149 (MQ=255)
                                 cTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGc                                                                                                                 <  6:385487/149‑1 (MQ=255)
                                    cTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGcctg                                                                                                              >  5:237612/1‑147 (MQ=255)
                                     tCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGt                                                                                                             >  2:183274/1‑149 (MQ=255)
                                       gCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGt                                                                                                             <  8:122299/147‑1 (MQ=255)
                                             cACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGaaa                                                                                                      <  1:383843/148‑1 (MQ=255)
                                             cACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAAc                                                                                                     <  3:375576/149‑1 (MQ=255)
                                                      gatgatGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAACCGCTCACCCg                                                                                            >  4:284564/1‑149 (MQ=255)
                                                      gatgatGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCg                                                                                            >  6:61800/1‑149 (MQ=255)
                                                      gatgatGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCg                                                                                            >  5:109853/1‑149 (MQ=255)
                                                             aaGTACGCACCACGCGTACCCACTGCACTAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCtgctgc                                                                                     <  3:313639/149‑1 (MQ=255)
                                                               gTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCAt                                                                                   >  3:355600/1‑149 (MQ=255)
                                                               gTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCAt                                                                                   >  2:46876/1‑149 (MQ=255)
                                                                  ataaccacGCGTATCCACTGCAAAAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGAAAGACGCTGCGGTCTTTGGCCTAGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCAACCGATGCTGCATGAt                                                                                >  4:249682/4‑149 (MQ=255)
                                                                   tCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATc                                                                               >  8:405604/1‑149 (MQ=255)
                                                                     accacGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCg                                                                              <  4:211606/148‑1 (MQ=255)
                                                                      ccacGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGtg                                                                            <  4:355600/149‑1 (MQ=255)
                                                                          gcgTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGtgctg                                                                         >  1:90747/1‑149 (MQ=255)
                                                                          gcgTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGtgctgc                                                                        <  4:339157/149‑1 (MQ=255)
                                                                          gcgTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGtgctg                                                                         >  8:227043/1‑148 (MQ=255)
                                                                             tATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAAc                                                                     <  2:90747/149‑1 (MQ=255)
                                                                                ccACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCgctg                                                                                         >  3:407728/1‑126 (MQ=255)
                                                                                ccACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCgctg                                                                                         <  4:407728/126‑1 (MQ=255)
                                                                                 cACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCt                                                                 >  8:347470/1‑149 (MQ=255)
                                                                                  aCTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTc                                                                >  8:40689/1‑149 (MQ=255)
                                                                                  aCTGCAATAGCGGGCAAGGGCGTCTGCTGAAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTc                                                                <  2:321616/149‑1 (MQ=255)
                                                                                    tGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCgg                                                              >  1:405661/1‑149 (MQ=255)
                                                                                    tGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCgg                                                              >  4:369418/1‑149 (MQ=255)
                                                                                    tGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCgg                                                              >  7:364344/1‑149 (MQ=255)
                                                                                      cAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCgg                                                              >  4:228940/1‑147 (MQ=255)
                                                                                      cAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGc                                                             <  3:369418/148‑1 (MQ=255)
                                                                                        ttAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTGACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCaga                                                          <  6:271577/148‑1 (MQ=255)
                                                                                         tAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCaa                                                                      >  1:365009/1‑136 (MQ=255)
                                                                                         tAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCaa                                                                      <  2:365009/136‑1 (MQ=255)
                                                                                         tAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCa                                                            >  8:32727/1‑146 (MQ=255)
                                                                                         tAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCa                                                            >  3:318702/1‑146 (MQ=255)
                                                                                              ggCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACgg                                                     >  5:385215/1‑148 (MQ=255)
                                                                                              ggCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGc                                                    <  2:405661/149‑1 (MQ=255)
                                                                                              ggCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGc                                                    <  3:228940/149‑1 (MQ=255)
                                                                                                 aaGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCaa                                                  >  1:57425/1‑148 (MQ=255)
                                                                                                 aaGGGCGTCTGCTGCAGCCGGGCAGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCgc                                                                                           <  7:300012/107‑1 (MQ=255)
                                                                                                 aaGGGCGTCTGCTGCAGCCGGGCAGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCgc                                                                                           >  8:300012/1‑107 (MQ=255)
                                                                                                     gCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCa                                             <  6:83772/149‑1 (MQ=255)
                                                                                                        tCTGCTGCAGCCGGGCCGGATCGACGATGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTTCAACTGCTCGGCAGAGACGGCAACACCATgg                                          >  8:6783/1‑149 (MQ=255)
                                                                                                         ctgctgCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGt                                         <  4:318702/149‑1 (MQ=255)
                                                                                                              gcAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAAcgc                                    >  3:371619/1‑149 (MQ=255)
                                                                                                               cAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAAcgcg                                   >  2:178875/1‑149 (MQ=255)
                                                                                                                aGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGTTCGGCAGAGACGGCAACACCATGGTAACgcgc                                  <  1:46876/149‑1 (MQ=255)
                                                                                                                  ccgggccggATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACgcgc                                  >  8:40519/1‑147 (MQ=255)
                                                                                                                   cgggccggATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAg                               <  7:40689/149‑1 (MQ=255)
                                                                                                                        ccggATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGt                            >  7:58240/1‑148 (MQ=255)
                                                                                                                        ccggATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTg                           <  5:61800/149‑1 (MQ=255)
                                                                                                                             tCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGAt                         >  8:287939/1‑146 (MQ=255)
                                                                                                                               gACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGa                     <  6:109853/148‑1 (MQ=255)
                                                                                                                               gACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGa                     <  4:199030/148‑1 (MQ=255)
                                                                                                                                aCGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGaca                   >  7:186228/1‑149 (MQ=255)
                                                                                                                                aCGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGaca                   >  3:306885/1‑149 (MQ=255)
                                                                                                                                    tGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGacacac                >  1:55622/1‑148 (MQ=255)
                                                                                                                                      cGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACAcc               <  4:306885/147‑1 (MQ=255)
                                                                                                                                       ggTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGc            >  2:233766/1‑149 (MQ=255)
                                                                                                                                         aacttGGCCTCAAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGCCCCACCAGcc           >  2:185464/4‑148 (MQ=255)
                                                                                                                                         aacttGGCCTCAAGGCCAAGTTGCTGAAATGCCAGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTATCGCGCCAGCGTGATATGGCCCACCAGcc           >  8:404521/4‑148 (MQ=255)
                                                                                                                                           tttGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCAc        >  5:185722/1‑149 (MQ=255)
                                                                                                                                           cttGGCCTCAAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCa         >  6:243075/2‑148 (MQ=255)
                                                                                                                                              ggCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACttt     >  3:88006/1‑149 (MQ=255)
                                                                                                                                              ggCCTCAAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACtt      <  5:243075/148‑1 (MQ=255)
                                                                                                                                               gCCTTGAGGCGAAGTTGCTGAAATGCCCGAAGTAGTTGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACttt     <  3:249682/148‑1 (MQ=255)
                                                                                                                                               gCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACttt     <  7:40519/148‑1 (MQ=255)
                                                                                                                                               gCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACttt     <  2:57425/148‑1 (MQ=255)
                                                                                                                                                ccTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACTTTc    >  3:276415/1‑148 (MQ=255)
                                                                                                                                                 cTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACTTTCaa  >  2:188466/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGGCGCGCGGGAATGGGGGTCCTGCACCAGAACCTGCTCGCCCTCCACGCGCGCGATGATGAAGTACTCACCACGCGTATCCACTGCAATAGCGGGCAAGGGCGTCTGCTGCAGCCGGGCCGGATCGACGCTGCGGTGTTTGGCCTTGAGGCCAAGTTGCTGAAATGCCCGAAGTAGTGGGCCGGTGGAGAAACGCTCACCCGCTGCTGCATGATCGTGCTGCAACTGCTCGGCAGAGACGGCAACACCATGGTAACGCGCCAGCGTGATAAGACACACCAGCCCACTTTCAA  >  NC_002947/5607761‑5608053

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: