Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,606,901 A→C 9.2% L413W (TTG→TGG)  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,606,9010AC 9.2% 177.3 / 4.4 61L413W (TTG→TGG) cyaBcyclolysin secretion/processing ATP‑binding protein CyaB
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base C (2/4);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGATACCCCTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAGGTTCCACGGCCATCGAC  >  NC_002947/5606757‑5607047
                                                                                                                                                 |                                                                                                                                                  
aTCGATACCCCTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGt                                                                                                                                                <  4:360193/149‑1 (MQ=255)
aTCGATACCCCTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGt                                                                                                                                                <  8:157533/149‑1 (MQ=255)
  cGATACCCCTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTgg                                                                                                                                              <  2:319683/149‑1 (MQ=255)
       ccccTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACCTTGTGGCtgg                                                                                                                                          <  1:350362/148‑3 (MQ=255)
        cccTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACCTTGTGGCtgggg                                                                                                                                        <  5:201760/149‑5 (MQ=255)
        cccTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGt                                                                                                                                                >  4:146573/1‑141 (MQ=255)
        cccTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGt                                                                                                                                                <  3:146573/141‑1 (MQ=255)
        cccTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCg                                                                                                                                        >  4:39951/1‑149 (MQ=255)
          cTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGt                                                                                                                                       >  4:126126/1‑148 (MQ=255)
                 ctgGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACaa                                                                                                                                  <  3:39951/146‑1 (MQ=255)
                    gAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGGGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCCCAAGGTGGCTACCGTGACAAGCTTGc                                                                                                                            >  3:210350/1‑149 (MQ=255)
                         tGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACCTTGTGGCtgg                                                                                                                                          <  1:126593/130‑3 (MQ=255)
                         tGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATc                                                                                                                       <  6:74646/149‑1 (MQ=255)
                         tGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCACCCCCCATGTGGCtgg                                                                                                                                          >  2:126593/1‑128 (MQ=255)
                                   ttGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATACAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACg                                                                                                              >  1:285608/1‑149 (MQ=255)
                                   ttGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGc                                                                                                             >  7:257453/1‑149 (MQ=255)
                                   ttGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGc                                                                                                             <  4:277477/149‑1 (MQ=255)
                                       gCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCttt                                                                                                          >  6:290886/1‑148 (MQ=255)
                                         cAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCAt                                                                                                       >  3:20261/1‑149 (MQ=255)
                                                atgatgGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTCCCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGc                                                                                                   >  5:27742/1‑146 (MQ=255)
                                                  gatgGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCgaga                                                                                               <  4:20261/148‑1 (MQ=255)
                                                   atgGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGAc                                                                                              <  3:172253/148‑1 (MQ=255)
                                                   atgGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGAc                                                                                              <  2:167414/148‑1 (MQ=255)
                                                         ggccacccgcgCCCGTCGCGCATGTTGAATGCGATCAGCTCGCCCACCGTGAGGGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAg                                                                                       <  8:257453/142‑1 (MQ=255)
                                                              cacGCGCCCTGCGCGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCt                                                                                  <  2:241969/149‑1 (MQ=255)
                                                               acGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTg                                                                                 >  6:210747/1‑149 (MQ=255)
                                                                  cgcCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGt                                                                              <  5:210747/149‑1 (MQ=255)
                                                                      cTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTg                                                                          <  5:290886/149‑1 (MQ=255)
                                                                              aTGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGc                                                                                                   <  1:324429/116‑1 (MQ=255)
                                                                              aTGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGc                                                                                                   >  2:324429/1‑116 (MQ=255)
                                                                                gTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTgg                                                                   >  6:177678/1‑146 (MQ=255)
                                                                                        ccaTCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATaa                                                        <  8:49862/147‑1 (MQ=255)
                                                                                        cGATCAGCTCGCCCACCGTCAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACCAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGCCGCTTTCATTGGCGAGACTGGATAGGTTCTGGGTGTTGAAACTGGCAGATTCATag                                                        >  8:19089/1‑148 (MQ=255)
                                                                                         gATCAGCTCGCCCAGCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTgg                                                                                >  7:314523/1‑124 (MQ=255)
                                                                                         gATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGCCCAGCCGGGCCCCCAGCCACAAGGTGGATAACGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCACATACATAAg                                                       >  7:329200/1‑149 (MQ=255)
                                                                                         gATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTgg                                                                                <  8:314523/124‑1 (MQ=255)
                                                                                         gATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGAt                                                              >  8:152918/1‑142 (MQ=255)
                                                                                         gATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGAt                                                              <  7:152918/142‑1 (MQ=255)
                                                                                         gATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAg                                                       >  8:239727/1‑149 (MQ=255)
                                                                                           tCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAg                                                                                       <  7:20067/115‑1 (MQ=255)
                                                                                           tCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAg                                                                                       >  8:20067/1‑115 (MQ=255)
                                                                                                   gCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACt                                              <  3:126126/149‑1 (MQ=255)
                                                                                                     ccACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGa                                            >  5:353829/1‑149 (MQ=255)
                                                                                                     ccACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGAGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTg                                             >  7:144622/1‑148 (MQ=255)
                                                                                                       aCCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGAt                                           >  8:349086/1‑148 (MQ=255)
                                                                                                                gggacTTCTATGACAAGACTGTCCCCCAGCCACAAGGTGGATACCGTGCCAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAgc                                  <  8:329200/144‑1 (MQ=255)
                                                                                                                    cTTCGATGACCAGCCGGGCCCCCAGCCACCAGGGGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCgg                             <  2:15420/149‑1 (MQ=255)
                                                                                                                    cTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCtt                                                                                                           >  4:75617/1‑71 (MQ=255)
                                                                                                                    cTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCtt                                                                                                           <  3:75617/71‑1 (MQ=255)
                                                                                                                    cTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGAt                                           >  2:298816/1‑135 (MQ=255)
                                                                                                                    cTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGAt                                           <  1:298816/135‑1 (MQ=255)
                                                                                                                    cTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCgg                             >  1:9127/1‑149 (MQ=255)
                                                                                                                        gATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCa                         >  6:115095/1‑149 (MQ=255)
                                                                                                                             ccAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGGCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGa                    <  2:285608/149‑1 (MQ=255)
                                                                                                                               aGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAgg                  <  6:353829/149‑1 (MQ=255)
                                                                                                                                  cGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAGGTTc               >  4:232785/1‑149 (MQ=255)
                                                                                                                                        cccAGCGAGATGGTGGCTACCGGGCCAAGCGTGCCGATCAGGCCGACGCTTGCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGAGACATTTGCGGCCATCTGATCATTCCAGCGCCGGGTCACCTGAGGTTCCACGGc         <  8:144622/149‑1 (MQ=255)
                                                                                                                                         ccAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCAGTGGCGAGACTGGATAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAGGTTCCACGGc         >  3:47067/1‑148 (MQ=255)
                                                                                                                                            gCCACAATGTGGCTACCGTGACAAGCTTGCAGATGAGGCCGATGCTGTCATTGGCGAGACTGGACAGGTTCTGGGGCGTGAACCTGGCAGATACATCAGCGGCCAACTGATCCTTCCCGCGCCGGGTCACCTGAGGTTCCACTGCCATc     <  4:210350/149‑1 (MQ=255)
                                                                                                                                            gCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAGGTTCCACGGCCATc     >  4:84590/1‑149 (MQ=255)
                                                                                                                                             ccACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGTGCCGGGTCACCTGAGGTTCCACGGCCATc     >  8:72113/1‑148 (MQ=255)
                                                                                                                                               acaAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAGGTTCCACGGCCATCGAc  >  7:135208/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATCGATACCCCTGTCTGCTGGAAGCTGGTCCAGAGTTGGGCCAGGCGGATGATGGGCTGCGACACGCGCCCTGCGAGCATGTTGAAGGCGATCAGCTCGCCCACCGTGAGTGCGCCTTCGATGACCAGCCGGGCCCCCAGCCACAAGGTGGCTACCGTGACAAGCTTGCCGATCAGGCCGACGCTTTCATTGGCGAGACTGGAAAGGTTCTGGGTCTTGAAACTGGCAGATACATAAGCGGCCAACTGATCATTCCAGCGCCGGGTCACCTGAGGTTCCACGGCCATCGAC  >  NC_002947/5606757‑5607047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: