Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,143 A→T 8.8% intergenic (+22/+46) ynfM → / ← PP_5073 efflux transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,1430AT 8.8% 191.4 / 8.4 68intergenic (+22/+46)ynfM/PP_5073efflux transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base T (2/4);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GGCGGGCGGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAC  >  NC_002947/5789998‑5790285
                                                                                                                                                  |                                                                                                                                               
ggcgggcgGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCAccc                                                                                                                                              <  7:101165/149‑1 (MQ=255)
  cgggcgggGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCca                                                                                                                                             >  3:6393/1‑148 (MQ=255)
   gggcgggGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCcac                                                                                                                                            >  7:259528/1‑148 (MQ=255)
         ggTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtgt                                                                                                                                     <  7:127785/149‑1 (MQ=255)
         ggTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtgt                                                                                                                                     <  1:86835/149‑1 (MQ=255)
          gtgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtgtg                                                                                                                                    <  4:6393/149‑1 (MQ=255)
           tgtTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTgg                                                                                                                                   <  2:41007/149‑1 (MQ=255)
                  gCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCgg                                                                                                                            <  2:149121/148‑1 (MQ=255)
                   cACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACACCCGCTCCCACAGGTTGTGGGTGCgggc                                                                                                                          <  3:114267/149‑2 (MQ=39)
                     ccAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                        <  7:280579/149‑1 (MQ=255)
                     ccAGTGGCGCTTGAACGGGATCTGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTTCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGgtgt                                                                                                                        <  2:239344/149‑1 (MQ=255)
                         tGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGAcc                                                                                                                    <  8:224562/149‑1 (MQ=255)
                           gggCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCa                                                                                                                   <  2:226701/148‑1 (MQ=255)
                            ggCTGGAACGGGATCGGGCTGTTCATTGGCAGGCTTGATGCCGGTGGCGTTGATGTTGGCTTTGCACCTGAGCAATTTGCCCCCCAAAACGGTGTGATGAATTCTTGGGCACCCCCCAACCCACAACCTGTGGGAGCGGGtgtgccaa                                                                                                                   >  1:31885/1‑144 (MQ=255)
                               tGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCTCCCACAGGTTGTgggtg                                                                                                                                >  3:67504/1‑130 (MQ=255)
                               tGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCTCCCACAGGTTGTgggtg                                                                                                                                <  4:67504/132‑3 (MQ=255)
                                 gAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGaa                                                                                                               >  3:71813/1‑146 (MQ=255)
                                 gAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAagag                                                                                                            >  2:175422/1‑149 (MQ=255)
                                 gAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAaga                                                                                                             <  4:71813/148‑1 (MQ=255)
                                  aaCGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAgg                                                                                                           <  1:72173/149‑1 (MQ=255)
                                       gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCat                                                                                                       >  5:83629/1‑148 (MQ=255)
                                       gATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCat                                                                                                       >  1:241760/1‑148 (MQ=255)
                                          cGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCat                                                                                                       >  3:31506/1‑145 (MQ=255)
                                            ggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtg                                                                                                                                      >  4:289355/1‑113 (MQ=255)
                                            ggCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCtg                                                                                                                                      <  3:289355/113‑1 (MQ=255)
                                            ggCTGTTCATTGGCAGCTTGCTGGCGGGGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGGGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCCGGGTGTGACCACGAAGAGGCCatcagat                                                                                                  <  1:175422/149‑1 (MQ=255)
                                             gCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatca                                                                                                <  7:88149/149‑1 (MQ=255)
                                               tGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGt                                                                                              >  7:155475/1‑149 (MQ=255)
                                                 ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                      >  3:124200/1‑139 (MQ=255)
                                                 ttCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatc                                                                                                      <  4:124200/139‑1 (MQ=255)
                                                         cAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGAtctc                                                                                    >  1:228556/1‑149 (MQ=255)
                                                            gCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCTATCAGATCAGTTGATGATCTCGa                                                                                  <  2:179990/149‑1 (MQ=255)
                                                            gCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGa                                                                                  <  8:155475/149‑1 (MQ=255)
                                                                 cTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAata                                                                             <  8:259528/149‑1 (MQ=255)
                                                                  tGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAatat                                                                            >  5:167891/1‑149 (MQ=255)
                                                                          gCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGGGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGTGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATatccacatcc                                                                    <  8:133648/149‑1 (MQ=255)
                                                                           ccTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCACTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatcag                                                                                               <  6:284800/121‑1 (MQ=255)
                                                                           ccTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCACTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCatcagatcag                                                                                               >  5:284800/1‑121 (MQ=255)
                                                                                   gtggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCgg                                                           >  5:309148/1‑149 (MQ=255)
                                                                                   gtggCGTTGCACCTGAGCAAGATGCCACCCAAAACGGTGTGATCCATTCGTCGGCACGCCCGCACCCACAACCTGAGGCAGCGGGTGTCACCACGAAGTGGCCAACAGATCAGTTGATGATCACGACTATATCCACATCCACTAcccag                                                           >  3:81991/1‑147 (MQ=255)
                                                                                    tggCGTTGCGCCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGGGACCACGAAGAGGACATCAGATCAGGTGATGAGCTCGACAATATCCACATCCACTTCCCGGc                                                          >  4:38546/1‑149 (MQ=255)
                                                                                     ggCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGTCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCt                                                         <  6:272900/149‑1 (MQ=255)
                                                                                      gCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCt                                                         <  2:228556/148‑1 (MQ=255)
                                                                                         ttGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                     <  1:210247/149‑1 (MQ=255)
                                                                                            cACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAgg                                                  >  1:178449/1‑149 (MQ=255)
                                                                                             aCCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                     >  8:27242/1‑145 (MQ=255)
                                                                                             aCCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCtcatc                                                     <  7:27242/145‑1 (MQ=255)
                                                                                               cTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCggg                                                                                                                           <  4:78458/73‑1 (MQ=21)
                                                                                               cTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCTCCCACAGGTTGTGGGTGCggg                                                                                                                           >  3:78458/1‑73 (MQ=21)
                                                                                               cTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGc                                               >  3:250495/1‑149 (MQ=255)
                                                                                                 caGCAAGCTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCtt                                             <  8:253082/148‑1 (MQ=255)
                                                                                                  aGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCGCGGCTCATCAGGTGCttt                                            >  3:42621/1‑149 (MQ=255)
                                                                                                   gCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCtttt                                           <  2:241760/149‑1 (MQ=255)
                                                                                                      aGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCa                                         >  5:124633/1‑148 (MQ=255)
                                                                                                       gTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAAc                                       <  2:31885/149‑1 (MQ=255)
                                                                                                            cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAc                                  >  8:97030/1‑149 (MQ=255)
                                                                                                            cACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCACGGCTCATCAGGTGCTTTTCAACCTCAc                                  <  4:42621/149‑1 (MQ=255)
                                                                                                                cAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAc                                  >  4:199761/1‑145 (MQ=255)
                                                                                                                cAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCAc                                  <  3:199761/145‑1 (MQ=255)
                                                                                                                          tgATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCAcc                    >  3:103133/1‑149 (MQ=255)
                                                                                                                          tgATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCAc                     <  6:309148/148‑1 (MQ=255)
                                                                                                                            aTGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCtt                  >  3:5386/1‑149 (MQ=255)
                                                                                                                             tGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                 <  3:135750/149‑1 (MQ=255)
                                                                                                                             tGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTg                 <  4:5386/149‑1 (MQ=255)
                                                                                                                                aTTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTc              <  4:31506/149‑1 (MQ=255)
                                                                                                                                      gggCACGCCCGCACCAACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCAACATCCACTTCCCGGCTCATCAGTTGCTTTTCAACATCACCGGTAAGCTTCAAATTGGTCTTGTCg        >  6:259856/1‑149 (MQ=255)
                                                                                                                                      gggCACGAACGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAAAATCCACATCCACTTACCGGCTCATCAGGAGCTTTTCAACCTCACCGGACAGCTTCACCTTGGTCTtctct        >  8:14772/1‑148 (MQ=255)
                                                                                                                                            gCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAc  <  6:167591/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGCGGGCGGGGTGTTCTGGCACCAGTGGGGCTGGAACGGGATCGGGCTGTTCATTGGCAGCTTGCTGGCGGTGGCCTTGATGGTGGCGTTGCACCTGAGCAAGTTGCCACCCAAAACGGTGTGATGCATTCGTGGGCACGCCCGCACCCACAACCTGTGGGAGCGGGTGTGACCACGAAGAGGCCATCAGATCAGTTGATGATCTCGACAATATCCACATCCACTTCCCGGCTCATCAGGTGCTTTTCAACCTCACCGGTCAGCTTCACCTTGGTCTTGTCGTTGAAC  >  NC_002947/5789998‑5790285

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: