Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,141,208 C→A 6.5% V118V (GTG→GTT)  PP_4525 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,141,2080CA 6.5% 199.5 / 3.7 62V118V (GTG→GTT) PP_4525LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base A (2/2);  total (30/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 = 6.47e-01

GATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCG  >  NC_002947/5141064‑5141346
                                                                                                                                                |                                                                                                                                          
gATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCt                                                                                                                                        >  2:234599/1‑149 (MQ=255)
 aTCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTGAACCTg                                                                                                                                       <  8:221684/149‑1 (MQ=255)
    aGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCg                                                                                                                                     <  3:234805/148‑1 (MQ=255)
     gTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTc                                                                                                                                            >  5:97985/1‑140 (MQ=255)
     gTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTc                                                                                                                                            <  6:97985/140‑1 (MQ=255)
     gTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCggg                                                                                                                                   <  5:100923/149‑1 (MQ=255)
      tACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCggg                                                                                                                                   >  5:138334/1‑148 (MQ=255)
              gCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTgg                                                                                                                          >  8:226930/1‑149 (MQ=255)
               cTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTgg                                                                                                                          <  6:274103/148‑1 (MQ=255)
               cTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGt                                                                                                                         <  7:226930/149‑1 (MQ=255)
                  aTCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTg                                                                                                                        <  6:99969/147‑1 (MQ=255)
                      cacaGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAgcg                                                                                                                  >  7:53921/1‑149 (MQ=255)
                       acaGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAgcgc                                                                                                                 >  5:95560/1‑149 (MQ=255)
                              cACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGc                                                                                                           >  1:44116/1‑148 (MQ=255)
                                 gCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACg                                                                                                       >  3:156714/1‑149 (MQ=255)
                                   cATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGcc                                                                                                     >  3:25107/1‑149 (MQ=255)
                                   cATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGcc                                                                                                     <  4:25107/149‑1 (MQ=255)
                                   cATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGAGCGCCAGCAAAGcc                                                                                                     >  6:2155/1‑149 (MQ=255)
                                    aTCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCg                                                                                                    <  8:112919/149‑1 (MQ=255)
                                       cAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGc                                                                                                  >  5:128552/1‑148 (MQ=255)
                                        aGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAgc                                                                                                                   <  3:98314/130‑1 (MQ=255)
                                        aGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAgc                                                                                                                   >  4:98314/1‑130 (MQ=255)
                                          cccGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGt                                                                                              >  4:305988/1‑149 (MQ=255)
                                            cGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTgg                                                                                            >  1:166054/1‑149 (MQ=255)
                                            cGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTgg                                                                                            <  6:95560/149‑1 (MQ=255)
                                               cAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTc                                                                                         <  7:31439/149‑1 (MQ=255)
                                                   gCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTcgccg                                                                                     >  8:6155/1‑149 (MQ=255)
                                                   gCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTcgccg                                                                                     >  5:95158/1‑149 (MQ=255)
                                                         tGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGc                                                                               <  4:156714/149‑1 (MQ=255)
                                                            ccagTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCATCGCCGGCAAGTGGGTCGCCGCCGTGCTGt                                                                            <  5:2155/146‑1 (MQ=255)
                                                                 ccACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTGAACCTGcaag                                                                                                                                   <  8:266337/89‑4 (MQ=255)
                                                                 ccACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTGAACCTGcaag                                                                                                                                   >  7:266337/1‑86 (MQ=255)
                                                                  cacaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATg                                                                       <  2:314677/148‑1 (MQ=255)
                                                                   acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTg                                                                                >  5:91902/1‑138 (MQ=255)
                                                                   acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTg                                                                                <  6:91902/138‑1 (MQ=255)
                                                                   acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATgg                                                                      >  7:138507/1‑148 (MQ=255)
                                                                   acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATgg                                                                      <  8:138507/148‑1 (MQ=255)
                                                                   acaAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGt                                                                     <  6:138334/149‑1 (MQ=255)
                                                                     aaGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGc                                                                               <  2:197957/137‑1 (MQ=255)
                                                                     aaGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGc                                                                               >  1:197957/1‑137 (MQ=255)
                                                                      agagCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTccc                                                                  <  7:6155/149‑1 (MQ=255)
                                                                      agagCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTccc                                                                  <  1:209244/149‑1 (MQ=255)
                                                                      agagCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTccc                                                                  <  3:305988/149‑1 (MQ=255)
                                                                        agcagcATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCa                                                                <  7:118871/149‑1 (MQ=255)
                                                                        agcagcATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCa                                                                <  5:302093/149‑1 (MQ=255)
                                                                                  aGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTc                                                      >  6:297746/1‑149 (MQ=255)
                                                                                    ggTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTccc                                                    <  2:44116/149‑1 (MQ=255)
                                                                                      tgtgGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCa                                                  <  6:128552/149‑1 (MQ=255)
                                                                                             cAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCt                                           >  1:127718/1‑149 (MQ=255)
                                                                                                    ccTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGc                                     >  6:19178/1‑148 (MQ=255)
                                                                                                        cAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCcac                                  <  6:95158/147‑1 (MQ=255)
                                                                                                           gCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCcacgcacg                             >  5:287120/1‑149 (MQ=255)
                                                                                                                   ccgcTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAg                     <  2:107219/149‑1 (MQ=255)
                                                                                                                           tgctgcCCGCACTTGCAGGTTCCCCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGGCGCCGCCGGGCTGTACATGGGCCCCAGCACGAAGTCCCCAGCCGCCTGACCGCCACGCACGgcg                          <  2:193475/136‑1 (MQ=255)
                                                                                                                           tgctgcACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTtcgt             >  3:62570/1‑149 (MQ=255)
                                                                                                                           tgctgcACGCACTTGCAGGTTCACCTGCGAGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGgcg                          >  1:193475/1‑136 (MQ=255)
                                                                                                                                cacgcacTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTcgcg         >  7:96005/1‑148 (MQ=255)
                                                                                                                                     acTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTccc                                                                  <  6:214994/86‑1 (MQ=255)
                                                                                                                                     acTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTccc                                                                  >  5:214994/1‑86 (MQ=255)
                                                                                                                                      cTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCg  <  2:127718/149‑1 (MQ=255)
                                                                                                                                       ttGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTcgc          >  6:155390/1‑140 (MQ=255)
                                                                                                                                       ttGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTcgc          <  5:155390/140‑1 (MQ=255)
                                                                                                                                       ttGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCg  >  7:11370/1‑148 (MQ=255)
                                                                                                                                       ccgCAGGTGAACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCCGCACGAAGTCCCCAGCCGGCTGACCGCCCCGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCg  >  6:283581/3‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCG  >  NC_002947/5141064‑5141346

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: ‑