Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,356,268 T→C 5.5% M194V (ATG→GTG) ‡ 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,2680TC 5.5% 215.4 / 2.2 72M194V (ATG→GTG) ‡PP_2952LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/32);  new base C (2/2);  total (38/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.74e-01

TCACATCCAGCCAGGACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCA  >  NC_002947/3356121‑3356413
                                                                                                                                                   |                                                                                                                                                 
tCACATCCAGCCAGGACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATg                                                                                                                                                  >  6:247818/1‑149 (MQ=255)
             ggACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGTTGCATGTACAGCAGGCACt                                                                                                                                     >  5:248269/1‑149 (MQ=255)
              gACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    <  8:209052/149‑1 (MQ=255)
                 ttGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGTACATGCAGCAggggc                                                                                                                                      <  6:341447/144‑4 (MQ=255)
                 ttGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGTACATGCAGCAggggc                                                                                                                                      >  5:341447/1‑141 (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                                                                                                                                    >  4:134983/1‑122 (MQ=255)
                                         aCGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTg                                                                                                                                    <  3:134983/122‑1 (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                                                                                                                          <  6:301047/113‑1 (MQ=255)
                                                            tatCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGc                                                                                                                          >  5:301047/1‑113 (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                                                                 <  3:8594/138‑1 (MQ=255)
                                                                                            aCTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGccacc                                                                 >  4:8594/1‑138 (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                  <  5:26521/149‑1 (MQ=255)
                                                                                                                                tAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGc                  <  3:122060/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                        <  5:292259/126‑1 (MQ=255)
                                                                                                                                                 cATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCtatata                        >  6:292259/1‑126 (MQ=255)
                                                                                                                                                 cATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCa  >  4:87048/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCACATCCAGCCAGGACTTGTAGGTTATCCCTTCGCCGGCCACGGCCCAGCGCCGGGCGATATCGGCGTCGTCGCAGAGCAGGGCGCCACGCACTTGGACGGTCTGCGTTCCCAGTTGCCACTTGTCGTAGACCCGCCCCTGCTGCATGTACAGCAGGCACTGGTGTTCGAGCAGATCTCCAGGCGCTTGGACGGTGCCTTGCCGCGCCAGGTACTCGGGCGAGGCCACCAGCACCTTGCGGTTCCAAGGTGCCAGTGGCAGGGCTATATAGTTGGCGTCGGCGTTCAGGCCA  >  NC_002947/3356121‑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: ‑