Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,469,842 T→C 9.0% intergenic (+11/+24) vacJ → / ← PP_2164 lipoprotein VacJ/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,469,8420TC 9.0% 186.1 / 7.1 68intergenic (+11/+24)vacJ/PP_2164lipoprotein VacJ/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (3/3);  total (34/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 = 6.63e-01

ACCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCG  >  NC_002947/2469695‑2469985
                                                                                                                                                   |                                                                                                                                               
aCCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTg                                                                                                                                                <  6:247426/149‑1 (MQ=255)
 ccTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACGGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTgg                                                                                                                                               >  5:252213/1‑149 (MQ=255)
 ccTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTgg                                                                                                                                               <  4:47905/149‑1 (MQ=255)
        gACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGc                                                                                                                                         <  6:117456/148‑1 (MQ=255)
        gACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGct                                                                                                                                        >  7:98259/1‑148 (MQ=255)
            cccGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGc                                                                                                                                    >  1:301029/1‑149 (MQ=255)
             ccGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTgggg                                                                                                                                             <  6:109434/139‑1 (MQ=255)
             ccGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTgggg                                                                                                                                             >  5:109434/1‑139 (MQ=255)
                  ccAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCCGAGGCCGCTTTGCGGcccc                                                                                                                              <  8:98259/149‑3 (MQ=255)
                           tctcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAAt                                                                                                                     >  5:330369/1‑149 (MQ=255)
                             tcGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAAt                                                                                                                     <  7:234074/147‑1 (MQ=255)
                               ggCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCAt                                                                                                                 <  2:66769/149‑1 (MQ=255)
                                   gAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGc                                                                                                              >  4:253901/1‑148 (MQ=255)
                                   gAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGc                                                                                                              >  6:176298/1‑148 (MQ=255)
                                        gCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAAt                                                                                                        <  3:233313/149‑1 (MQ=255)
                                         cTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATc                                                                                                       >  1:37970/1‑149 (MQ=255)
                                              tCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCggg                                                                                                                           <  4:23011/124‑1 (MQ=255)
                                                aGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGAGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGcc                                                                                                >  7:200572/1‑149 (MQ=255)
                                                   gTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCg                                                                                              >  7:127472/1‑148 (MQ=255)
                                                          gTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGc                                                                                                    <  7:113101/135‑1 (MQ=255)
                                                          gTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGc                                                                                                    >  8:113101/1‑135 (MQ=255)
                                                          gTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctg                                                                                      >  6:297907/1‑149 (MQ=255)
                                                          gTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTctg                                                                                      <  8:74147/149‑1 (MQ=255)
                                                              aTCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCtgctgc                                                                                  >  4:3264/1‑149 (MQ=255)
                                                                  tCATTCGCAATGCCTACCTTCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAATAGGCCTCGGGTTAAATGAATGGCAAGAATAAGCAGCCCGAAATTCTTCTgacctc                                                                              >  5:247997/1‑147 (MQ=255)
                                                                  tCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgc                                                                              >  2:335543/1‑149 (MQ=255)
                                                                     ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGc                                                                                                    <  3:320727/124‑1 (MQ=255)
                                                                     ttCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGc                                                                                                    >  4:320727/1‑124 (MQ=255)
                                                                         aaaTGCCTAACTGAAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTAGGGTTCAATGCATGGCCAGAATCCGCATCCCGAACTTCTGCTGCCCGCCCGGCTg                                                                       <  6:247997/148‑1 (MQ=255)
                                                                          aaTGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTgg                                                                      >  1:329342/1‑149 (MQ=255)
                                                                          aaTGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTgg                                                                      >  7:31748/1‑149 (MQ=255)
                                                                           aTGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATGTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAGTGCATGGCCAGTATCCGCAGCCGGAACTTCTGCTGCCCGCCCGGCTGGt                                                                     >  6:19101/1‑149 (MQ=255)
                                                                           aTGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCTGCTGGt                                                                     <  2:301029/149‑1 (MQ=255)
                                                                              ccTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCg                                                                   <  6:330369/148‑1 (MQ=255)
                                                                                  ccTGCAAAACCGCGAGTTAAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCaa                                                               <  8:309506/148‑1 (MQ=255)
                                                                                        aaaCCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATAc                                                        >  6:282897/1‑149 (MQ=255)
                                                                                         aaCCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCAGTGAACCAGGCGCCGCAAAGCGCCCTCGGGTTCAATGCATGGCAAGAATCCGCAGCCCGAACATCAGCTGCCCGCCCGGCTGGTCGGCAATCCATaa                                                        >  7:166685/1‑147 (MQ=255)
                                                                                                   tCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGc                                             >  4:49153/1‑149 (MQ=255)
                                                                                                    cAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGTCCAGAATCCGCAGCCCGAACTTCTGCTTCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGc                                             <  2:37970/148‑1 (MQ=255)
                                                                                                     aaGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCt                                            <  7:68792/148‑1 (MQ=255)
                                                                                                     aaGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                           <  3:3264/149‑1 (MQ=255)
                                                                                                       ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGa                                                                                                          <  4:289516/84‑1 (MQ=255)
                                                                                                       ggTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGa                                                                                                          >  3:289516/1‑84 (MQ=255)
                                                                                                        gTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATa                                                         >  6:32737/1‑132 (MQ=255)
                                                                                                        gTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATa                                                         <  5:32737/132‑1 (MQ=255)
                                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                           >  4:15466/1‑141 (MQ=255)
                                                                                                             ggATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTg                                           <  3:15466/141‑1 (MQ=255)
                                                                                                                tGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAAc                                                                                           <  6:159159/90‑1 (MQ=255)
                                                                                                                tGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAAc                                                                                           >  5:159159/1‑90 (MQ=255)
                                                                                                                  gCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGcccc                               <  8:31748/148‑1 (MQ=255)
                                                                                                                   cGAGGTCGAAGACGAGTTCTGATCTGTGAACGTGGGGCCGCATAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGGCCGAACTTGTGCTGCCCGCGCGGCTGGGCGGGAATCCATACCTCTTCGGTGCTGGCATGCAGCCCCtt                             <  8:166685/149‑1 (MQ=255)
                                                                                                                      ggTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAg                          <  3:253901/149‑1 (MQ=255)
                                                                                                                      ggTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAg                          >  5:298319/1‑149 (MQ=255)
                                                                                                                            aGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCg                     >  3:323169/1‑148 (MQ=255)
                                                                                                                             gacgacTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGa                   <  5:282897/149‑1 (MQ=255)
                                                                                                                                     cTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCTAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCa                                                           <  1:128518/101‑1 (MQ=255)
                                                                                                                                     cTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCTAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCa                                                           >  2:128518/1‑101 (MQ=255)
                                                                                                                                        tgcattgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTc                            <  1:261197/124‑1 (MQ=37)
                                                                                                                                        tgcattgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTc                            >  2:261197/6‑129 (MQ=37)
                                                                                                                                        tgcattgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTc        <  2:180534/144‑1 (MQ=38)
                                                                                                                                        aTCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgc                                                                              <  1:248943/79‑1 (MQ=255)
                                                                                                                                        aTCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTgcccgc                                                                              >  2:248943/1‑79 (MQ=255)
                                                                                                                                        aTCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGa           <  5:297907/146‑1 (MQ=255)
                                                                                                                                           attgAACCCGAGGCCGCTTTGCGGCCCCAGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGAt     >  7:119783/3‑149 (MQ=38)
                                                                                                                                             tgAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGa      <  6:298319/146‑1 (MQ=255)
                                                                                                                                               aaCCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCg  >  4:163645/1‑148 (MQ=255)
                                                                                                                                                  cTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGa           <  6:190409/136‑1 (MQ=255)
                                                                                                                                                  cTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGa           >  5:190409/1‑136 (MQ=255)
                                                                                                                                                   tGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg             <  6:261471/133‑1 (MQ=255)
                                                                                                                                                   tGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACg             >  5:261471/1‑133 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCTGATCGACACCCGCGCCAGCCTGCTCTCGGCCGAAAAGCTGATTCAGGGTGACAAGTACATCTTCATTCGCAATGCCTACCTGCAAAACCGCGAGTTCAAGGTGAAGGATGGCGAGGTCGAAGACGACTTCTGATCTGTGAACCTGGGGCCGCAAAGCGGCCTCGGGTTCAATGCATGGCCAGAATCCGCAGCCCGAACTTCTGCTGCCCGCCCGGCTGGTCGGCAATCCATACCACTTCGGTGCTGGCATGCAGCCCCTTCAGCGCCGGATGGACGGAGTCGATCCG  >  NC_002947/2469695‑2469985

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: