Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,356,271 A→G 5.6% Y193H (TAC→CAC)  PP_2952 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,356,2710AG 5.6% 208.6 / 2.2 71Y193H (TAC→CAC) PP_2952LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base G (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GGACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCA  >  NC_002947/3356134‑3356413
                                                                                                                                         |                                                                                                                                              
ggACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGTTGCATGTACAGCAGGCACt                                                                                                                                     >  5:248269/1‑149 (MQ=255)
 gACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    <  8:209052/149‑1 (MQ=255)
    ttGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGTACATGCAGCAggggc                                                                                                                                      >  5:341447/1‑141 (MQ=255)
    ttGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGTACATGCAGCAggggc                                                                                                                                      <  6:341447/144‑4 (MQ=255)
    ttGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGgtg                                                                                                                                 >  7:42330/1‑149 (MQ=255)
       tAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTc                                                                                                                              >  3:70797/1‑149 (MQ=255)
       tAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTc                                                                                                                              >  5:289056/1‑149 (MQ=255)
         ggTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGTACATGCAGCAggggc                                                                                                                                      >  8:236681/1‑136 (MQ=255)
         ggTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGTACATGCAGCAggggc                                                                                                                                      <  7:236681/139‑4 (MQ=255)
         ggTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGa                                                                                                                            <  8:42330/149‑1 (MQ=255)
                     gCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    <  1:8916/129‑1 (MQ=255)
                     gCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    >  2:8916/1‑129 (MQ=255)
                      ccGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAgg                                                                                                               >  4:260084/1‑149 (MQ=255)
                        ggCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGgcg                                                                                                             >  4:142364/1‑149 (MQ=255)
                          ccACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCt                                                                                                           >  4:5315/1‑149 (MQ=255)
                           cACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCtt                                                                                                          <  8:90450/149‑1 (MQ=255)
                            aCGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    <  3:134983/122‑1 (MQ=255)
                            aCGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    >  4:134983/1‑122 (MQ=255)
                                cccAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACg                                                                                                     >  6:64930/1‑149 (MQ=255)
                                 ccAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAg                                                                                                                        <  6:243755/129‑1 (MQ=255)
                                 ccAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAg                                                                                                                        >  5:243755/1‑129 (MQ=255)
                                 ccAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACgg                                                                                                    <  5:63133/149‑1 (MQ=255)
                                   aGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTg                                                                                                  <  5:64930/149‑1 (MQ=255)
                                        cGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTg                                                                                             >  4:122060/1‑149 (MQ=255)
                                               tatCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGc                                                                                                                          >  5:301047/1‑113 (MQ=255)
                                               tatCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGc                                                                                                                          <  6:301047/113‑1 (MQ=255)
                                                atCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCgcgc                                                                                       >  4:27469/1‑147 (MQ=255)
                                                  cGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAgg                                                                                   >  7:273958/1‑149 (MQ=255)
                                                    gCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTa                                                                                 >  2:312512/1‑149 (MQ=255)
                                                      gtcgtcGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGTTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACt                                                                               <  6:248269/149‑1 (MQ=255)
                                                      gtcgtcGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACt                                                                               >  1:148405/1‑149 (MQ=255)
                                                       tcgtcgCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTc                                                                              <  6:191430/149‑1 (MQ=255)
                                                              agagCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAg                                                                       <  6:289056/149‑1 (MQ=255)
                                                                  cAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGCCGGTGCATTGCCGCGCCCGGTACTCGGGCGCGCcca                                                                   >  5:280112/1‑149 (MQ=255)
                                                                    gggCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGccac                                                                  >  5:305901/1‑148 (MQ=255)
                                                                          cACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCAcc                                                           >  1:267531/1‑149 (MQ=255)
                                                                           aCGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCgg                                                                            >  5:329648/1‑131 (MQ=255)
                                                                           aCGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCgg                                                                            <  6:329648/131‑1 (MQ=255)
                                                                           aCGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCt                                                          >  4:246478/1‑149 (MQ=255)
                                                                               aCTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGccacc                                                                 >  4:8594/1‑138 (MQ=255)
                                                                               aCTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGccacc                                                                 <  3:8594/138‑1 (MQ=255)
                                                                                cTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCgg                                                     <  4:70797/149‑1 (MQ=255)
                                                                                   ggACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGtt                                                   <  2:8438/148‑1 (MQ=255)
                                                                                       ggTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAg                                              <  2:267531/149‑1 (MQ=255)
                                                                                           tGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAg                                                                       >  5:296516/1‑120 (MQ=255)
                                                                                           tGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAg                                                                       <  6:296516/120‑1 (MQ=255)
                                                                                            gCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTg                                                                                                         <  8:273905/85‑1 (MQ=255)
                                                                                            gCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTg                                                                                                         >  7:273905/1‑85 (MQ=255)
                                                                                              gTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAg                                       <  3:5315/149‑1 (MQ=255)
                                                                                                      ttGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCtt                                                                                              >  1:51608/1‑86 (MQ=255)
                                                                                                      ttGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCtt                                                                                              <  2:51608/86‑1 (MQ=255)
                                                                                                        tccACTTGGCCTACCCCCGGCCCTGCTGCATGTAAACGAGGCACTGGTGTTCGCGCAGTTCTCCACGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGc                              <  6:280112/147‑1 (MQ=37)
                                                                                                        gCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGc                              <  8:273958/148‑1 (MQ=255)
                                                                                                         ccACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCta                            >  6:26521/1‑149 (MQ=255)
                                                                                                          cACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGcg                                                             <  4:29978/115‑2 (MQ=255)
                                                                                                          cACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGcg                                                             >  3:29978/1‑114 (MQ=255)
                                                                                                                tCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGtt                     <  3:246478/149‑1 (MQ=255)
                                                                                                                tCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGtt                     <  3:27469/149‑1 (MQ=255)
                                                                                                                   tAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCtata                          >  3:342811/1‑141 (MQ=255)
                                                                                                                   tAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCtata                          <  4:342811/141‑1 (MQ=255)
                                                                                                                   tAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGc                  <  3:122060/149‑1 (MQ=255)
                                                                                                                   tAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGc                  <  5:26521/149‑1 (MQ=255)
                                                                                                                    aGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCg                 >  3:298203/1‑149 (MQ=255)
                                                                                                                        ccGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACTGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCAACAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTg                    <  3:5012/142‑1 (MQ=255)
                                                                                                                        ccGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCgg             >  2:214427/1‑149 (MQ=255)
                                                                                                                        ccGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCAACAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTg                    >  4:5012/1‑142 (MQ=255)
                                                                                                                             ccTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCt                                                                                               >  4:320706/1‑62 (MQ=255)
                                                                                                                             ccTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCt                                                                                               <  3:320706/62‑1 (MQ=255)
                                                                                                                                    cATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCtatata                        >  6:292259/1‑126 (MQ=255)
                                                                                                                                    cATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCtatata                        <  5:292259/126‑1 (MQ=255)
                                                                                                                                    cATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCa  >  4:87048/1‑148 (MQ=255)
                                                                                                                                         |                                                                                                                                              
GGACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCA  >  NC_002947/3356134‑3356413

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: