Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,141,207 T→G 6.3% N119H (AAC→CAC)  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,2070TG 6.3% 199.8 / 4.4 64N119H (AAC→CAC) PP_4525LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base G (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

TGGTGATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCG  >  NC_002947/5141060‑5141346
                                                                                                                                                   |                                                                                                                                           
tGGTGATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGtt                                                                                                                                             >  1:314677/1‑148 (MQ=255)
    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)
                                                                                                                                                   |                                                                                                                                           
TGGTGATCAGTACCAGTTGCTCATCGCACAGCGGCACGCCATCCAGCCCGGCAAGGCTCGGTGGCCCGTCCACAAGAGCAGCATCCAGGGTGTGGTTCAACAGCCCTTCAAGCAACTCGCCGCTGGGTGCTGCACGCACTTGCAGGTTCACCTGCGGGTAAGACTGGTGATAGCGCGCCAGCAACGCCGGCAAGTGGGTCGCCGCCGTGCTGTACATGGTCCCCAGCACGAAGTCCCCAGCCGGCTGACCGCCACGCACGGCGCCCAGCGCTTCGTCGCGCAAGGCG  >  NC_002947/5141060‑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: ‑