Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,462,394 G→C 6.2% intergenic (+17/+17) PP_4798 → / ← PP_4799 membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,462,3940GC 6.2% 198.6 / 3.7 64intergenic (+17/+17)PP_4798/PP_4799membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base C (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

CGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGG  >  NC_002947/5462250‑5462536
                                                                                                                                                |                                                                                                                                              
cGGGTAATGCCGATGGCATTTTCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                                                                                                                                             <  3:210374/148‑1 (MQ=255)
cGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                                                                                                                                             <  3:97818/148‑1 (MQ=255)
cGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                                                                                                                                             <  5:269393/148‑1 (MQ=255)
    tAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGccc                                                                                                                                        >  4:257826/1‑149 (MQ=255)
      aTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGCCCCg                                                                                                                                      <  5:284958/149‑1 (MQ=255)
       tGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCgg                                                                                                                                     <  5:105956/149‑1 (MQ=255)
         ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCGGCGCCAGCAGTAAGCATTCCGGGGCCGCGGCGCGGCCCCGGt                                                                                                                                    >  7:121132/1‑148 (MQ=255)
          cGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTaa                                                                                                                                  <  3:257826/149‑1 (MQ=255)
            aTGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtc                                                                                                                                >  7:207600/1‑149 (MQ=255)
            aTGGCATAATAGGGGCCAATACCCCCAAGCCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGATAGCATGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtc                                                                                                                                <  8:260982/149‑1 (MQ=255)
                cATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCAc                                                                                                                            >  6:304028/1‑149 (MQ=255)
                      cGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCt                                                                                                                       >  4:165805/1‑148 (MQ=255)
                          gCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                    >  8:228365/1‑147 (MQ=255)
                            cAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                    >  5:92273/1‑145 (MQ=255)
                               tACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATcc                                                                                                             <  8:207600/149‑1 (MQ=255)
                                    gCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTc                                                                                                        <  4:224742/149‑1 (MQ=255)
                                          cAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtc                                                                                                                                <  6:180686/119‑1 (MQ=255)
                                          cAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtc                                                                                                                                >  5:180686/1‑119 (MQ=255)
                                            aTCCGCAATGCCCAGCTAGGTCTGGGATGGCCGGCTGATGTCTACCCGACCGACAAGCTGCTTTAGAGCCTGCGCCATCAGTAAGCATTCCGGGGACGCTGCGCGGCCCCGGTAATCTGACGCCCTGAACAGATCCTGCTCCAGCACCa                                                                                                <  7:149750/149‑1 (MQ=255)
                                               cGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAgt                                                                                              <  5:304028/148‑1 (MQ=255)
                                                gCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAgtgt                                                                                            >  8:42524/1‑149 (MQ=255)
                                                     gCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGcc                                                                                                                                         >  3:287960/1‑99 (MQ=255)
                                                     gCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGcc                                                                                                                                         <  4:287960/99‑1 (MQ=255)
                                                                gCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCTCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGc                                                                                        <  1:44003/137‑1 (MQ=255)
                                                                gCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCTCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGc                                                                                        >  2:44003/1‑137 (MQ=255)
                                                                       tGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCa                                                                     <  7:228365/149‑1 (MQ=255)
                                                                         gCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAc                                                                    <  1:134799/148‑1 (MQ=255)
                                                                         gCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAc                                                                    <  7:42524/148‑1 (MQ=255)
                                                                          ccGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGa                                                                                     >  3:51891/1‑130 (MQ=255)
                                                                          ccGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGa                                                                                     <  4:51891/130‑1 (MQ=255)
                                                                            ggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGa                                                                <  6:92273/149‑1 (MQ=255)
                                                                                gATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGccc                                                            >  7:276907/1‑149 (MQ=255)
                                                                                    gCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCCGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGAGTTTGCTGAACACTCTCCAGCCTCACCCGAGCCCCCac                                                        >  2:267591/1‑148 (MQ=255)
                                                                                         ccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAgg                                                       >  7:92910/1‑145 (MQ=255)
                                                                                         ccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAgg                                                       <  8:92910/145‑1 (MQ=255)
                                                                                           gACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGcccc                                                           >  3:289503/1‑139 (MQ=255)
                                                                                           gACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGcccc                                                           <  4:289503/139‑1 (MQ=255)
                                                                                            aCCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGc                                                <  8:276907/149‑1 (MQ=255)
                                                                                              ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAACCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGcc                                                             <  3:177033/134‑1 (MQ=255)
                                                                                              ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAACCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGcc                                                             >  4:177033/1‑134 (MQ=255)
                                                                                                acaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                             >  7:252861/1‑148 (MQ=255)
                                                                                                acaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCACCCTCACCCGAGCCACCAGCGGCAAGCATa                                             >  6:79099/1‑148 (MQ=255)
                                                                                                 caAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACa                                           <  1:267591/149‑1 (MQ=255)
                                                                                                  aaGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACAt                                          >  8:283991/1‑149 (MQ=255)
                                                                                                        ctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                      <  3:165805/147‑1 (MQ=255)
                                                                                                        ctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGt                                    <  1:293136/149‑1 (MQ=255)
                                                                                                        ctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGt                                    <  4:253629/149‑1 (MQ=255)
                                                                                                            agagCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACCCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTcaca                                <  5:79099/149‑1 (MQ=255)
                                                                                                             gagCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAg                               >  6:265312/1‑149 (MQ=255)
                                                                                                                ccTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCAt                                              >  4:143564/1‑131 (MQ=255)
                                                                                                                ccTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCAt                                              <  3:143564/131‑1 (MQ=255)
                                                                                                                ccTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGc                            >  2:249988/1‑149 (MQ=255)
                                                                                                                        agcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAg                                                        >  6:105940/1‑113 (MQ=255)
                                                                                                                        agcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAg                                                        <  5:105940/113‑1 (MQ=255)
                                                                                                                        agcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCCCAGTGc                            >  3:181469/1‑141 (MQ=255)
                                                                                                                        agcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGc                            <  4:181469/141‑1 (MQ=255)
                                                                                                                               aGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAAccc             >  5:171926/1‑149 (MQ=255)
                                                                                                                               aGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAAccc             >  4:187416/1‑149 (MQ=255)
                                                                                                                                 cATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTc                                                                             <  1:115830/83‑1 (MQ=255)
                                                                                                                                 cATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                      <  6:255475/122‑1 (MQ=255)
                                                                                                                                 cATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                      >  5:255475/1‑122 (MQ=255)
                                                                                                                                 cATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGc           >  4:95462/1‑149 (MQ=255)
                                                                                                                                  aTTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGc           <  7:259985/148‑1 (MQ=255)
                                                                                                                                  aTTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGc           >  8:259985/1‑148 (MQ=255)
                                                                                                                                      cGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCa      >  1:234089/1‑149 (MQ=255)
                                                                                                                                          gCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg  <  7:283991/149‑1 (MQ=255)
                                                                                                                                           ccGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg  >  8:17519/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGG  >  NC_002947/5462250‑5462536

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