Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,462,394 G→C 8.7% 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 8.7% 233.3 / 7.4 77intergenic (+17/+17)PP_4798/PP_4799membrane‑bound lytic murein transglycosylase/muramoyltetrapeptide carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base C (3/4);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCA  >  NC_002947/5462246‑5462538
                                                                                                                                                      |                                                                                                                                                
gatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTg                                                                                                                                                  >  8:219217/1‑149 (MQ=255)
   cccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc                                                                                                                                               >  4:138461/6‑149 (MQ=255)
    ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg                                                                                                                                              >  1:375146/5‑149 (MQ=255)
       gTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGcc                                                                                                                                           >  1:258827/1‑149 (MQ=255)
       gTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGcc                                                                                                                                           <  2:258827/149‑1 (MQ=255)
          aTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCg                                                                                                                                        <  7:144734/149‑1 (MQ=255)
          aTGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGCCCCg                                                                                                                                        <  6:18759/149‑1 (MQ=255)
             ccGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGt                                                                                                                                      <  6:54129/148‑1 (MQ=255)
              cGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCGCAGCGGCCCCggaat                                                                                                                                    <  7:345731/149‑4 (MQ=255)
                     aTTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACg                                                                                                                             <  3:138461/149‑1 (MQ=255)
                      ttATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                            <  1:129237/149‑1 (MQ=255)
                      ttATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGc                                                                                                                            >  6:94955/1‑149 (MQ=255)
                         tCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCt                                                                                                                         >  7:333180/1‑149 (MQ=255)
                         tCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCt                                                                                                                         <  5:346309/149‑1 (MQ=255)
                            gggCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                      <  2:55007/149‑1 (MQ=255)
                               ccAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                      >  6:120019/1‑146 (MQ=255)
                               ccAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGCATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                      >  6:272016/1‑146 (MQ=255)
                                cAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAAtc                                                                                                                                  <  7:361522/133‑1 (MQ=255)
                                 aaTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                      >  7:267896/1‑144 (MQ=255)
                                 aaTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGaa                                                                                                                      >  3:23363/1‑144 (MQ=255)
                                 aaTACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTCCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCGGTAATCTCACGCCCTGAACAGATc                                                                                                                <  8:6752/149‑1 (MQ=255)
                                  ttACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATc                                                                                                                <  1:274710/148‑1 (MQ=255)
                                      ccGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGc                                                                                                            >  2:9281/1‑149 (MQ=255)
                                        gCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTc                                                                                                          >  4:211819/1‑149 (MQ=255)
                                        gCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTc                                                                                                          >  5:174041/1‑149 (MQ=255)
                                               aaTCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCAcc                                                                                                   <  5:27658/149‑1 (MQ=255)
                                                      aaTGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGttt                                                                                            <  6:380244/149‑1 (MQ=255)
                                                       aTGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTg                                                                                           <  3:107739/149‑1 (MQ=255)
                                                         gCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCt                                                                                         >  6:366052/1‑149 (MQ=255)
                                                           ccAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTg                                                                                        >  7:183856/1‑148 (MQ=255)
                                                            cAGCAAGGGCTGGGATGGCCGGCGGATGGATACCCGACCCACAAGCTGCTTGAGAGACTCGACGCCAGCAGTAAGCATTCCGGGCCCGCTGCGCGGCCCCTGTAATCTCACGCCCAGAACAGATCCTGCTCCAGCACCACTGTTTGCt                                                                                         >  1:111752/1‑148 (MQ=255)
                                                                     cTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTg                                                                                        <  6:175991/138‑1 (MQ=255)
                                                                     cTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTg                                                                                        >  5:175991/1‑138 (MQ=255)
                                                                     cTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCa                                                                             <  4:23363/149‑1 (MQ=255)
                                                                         gATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCt                                                                         >  3:323204/1‑149 (MQ=255)
                                                                            ggCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCCTTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGcc                                                                          <  7:203613/145‑1 (MQ=255)
                                                                            ggCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGcc                                                                          >  8:203613/1‑145 (MQ=255)
                                                                               cggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGTGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCAccc                                                                    >  5:78809/1‑148 (MQ=255)
                                                                                ggcggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGa                                                                  <  5:272016/149‑1 (MQ=255)
                                                                                   ggATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGcc                                                               <  6:174041/149‑1 (MQ=255)
                                                                                        gCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAg                                                          <  5:120019/149‑1 (MQ=255)
                                                                                            cccGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg                                                        <  7:67526/147‑1 (MQ=255)
                                                                                              cGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCaa                                                    <  8:209570/149‑1 (MQ=255)
                                                                                                 cccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCa                                                 <  8:333180/149‑1 (MQ=255)
                                                                                                 cccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCa                                                 <  8:267896/149‑1 (MQ=255)
                                                                                                  ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCa                                                 >  2:210798/1‑148 (MQ=255)
                                                                                                  ccACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCAt                                                <  2:375146/149‑1 (MQ=255)
                                                                                                   cacaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                               <  5:366052/149‑1 (MQ=255)
                                                                                                    acaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTGACGCCCTGAACAGATCCTGCGCCAGCACCAGTGTTTGCTGAACACTGTGCAGCCTCACCCCAGCACCCAGGGGCAAGCATa                                               >  4:214474/1‑148 (MQ=255)
                                                                                                    acaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCCTa                                               >  1:245983/1‑148 (MQ=255)
                                                                                                    acaAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATa                                               >  3:48627/1‑148 (MQ=255)
                                                                                                     caAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACa                                             <  4:323204/149‑1 (MQ=255)
                                                                                                        gctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTc                                          <  3:211819/149‑1 (MQ=255)
                                                                                                        gctgctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTc                                          >  2:9154/1‑149 (MQ=255)
                                                                                                            ctTGAGAGCCTGCTCCAGCATTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                        <  5:61546/147‑1 (MQ=255)
                                                                                                            ctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                        <  8:15364/147‑1 (MQ=255)
                                                                                                            ctTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                        <  8:183856/147‑1 (MQ=255)
                                                                                                             ttGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACt                                                                                  <  2:224969/104‑1 (MQ=255)
                                                                                                             ttGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACt                                                                                  >  1:224969/1‑104 (MQ=255)
                                                                                                                    ccTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGc                              <  4:48627/149‑1 (MQ=255)
                                                                                                                     cTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                        >  1:208437/1‑138 (MQ=255)
                                                                                                                     cTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCgg                                        <  2:208437/138‑1 (MQ=255)
                                                                                                                     cTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGcc                             >  5:82263/1‑149 (MQ=255)
                                                                                                                             cagcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCAt                                                <  7:105849/124‑1 (MQ=255)
                                                                                                                             cagcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCAt                                                >  8:105849/1‑124 (MQ=255)
                                                                                                                               gcagTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTgcgc                     <  7:219217/149‑1 (MQ=255)
                                                                                                                                   tAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg                                                        <  6:108258/110‑1 (MQ=255)
                                                                                                                                   tAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAggg                                                        >  5:108258/1‑110 (MQ=255)
                                                                                                                                    cccaccgccGGGGCCGCCCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAAc                 <  3:15458/141‑1 (MQ=255)
                                                                                                                                      gCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCg              >  1:293160/1‑149 (MQ=255)
                                                                                                                                           ccGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGc         >  6:94069/1‑149 (MQ=255)
                                                                                                                                            cGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGc                                                  >  5:137546/1‑107 (MQ=255)
                                                                                                                                            cGGGGCCGCGCAGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGc                                                  <  6:137546/107‑1 (MQ=255)
                                                                                                                                            cGGGGCCGCGCAGCGGCCCCGGGAATCTAACGCCCTGAACAGATCCTGCTCCAGCACCACACTTTGGTGAACAATGTCCACCGACACCCCCCCCCCCAGGGGCAACCATACAATCGGGTCAAAGTGCCCGCTGCCCCAAc                 >  4:15458/1‑140 (MQ=255)
                                                                                                                                                gCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg    >  4:364950/1‑149 (MQ=255)
                                                                                                                                                gCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCgg    >  5:281788/1‑149 (MQ=255)
                                                                                                                                                   gCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGCGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGca  >  5:250904/1‑148 (MQ=255)
                                                                                                                                                   gCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTcaca                                  <  5:335712/116‑1 (MQ=255)
                                                                                                                                                   gCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTcaca                                  >  6:335712/1‑116 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGGTAATCTCACGCCCTGAACAGATCCTGCTCCAGCACCAGTGTTTGCTGAACACTGTCCAGCCTCACCCGAGCCCCCAGGGGCAAGCATACATTCGGGTCACAGTGCCCGCTGCGCCAACCCGCCAGCACCGGCA  >  NC_002947/5462246‑5462538

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: