Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,049,913 C→T 15.2% intergenic (+46/‑55) quiA → / → oprB‑III quinate dehydrogenase/carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,049,9130CT15.2% 134.1 / 15.7 59intergenic (+46/‑55)quiA/oprB‑IIIquinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (5/4);  total (31/28)
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.77e-01

TCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAC  >  NC_002947/4049769‑4050059
                                                                                                                                                |                                                                                                                                                    
tcGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTg                                                                                                                                                    >  6:174338/1‑147 (MQ=255)
 cGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACCAGGCCCCTGTAg                                                                                                                                                 >  7:38028/1‑149 (MQ=255)
  gCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTa                                                                                                                                                  <  6:158763/147‑1 (MQ=255)
   ccGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               <  6:279362/149‑1 (MQ=255)
     gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGAGc                                                                                                                                             <  2:35346/149‑1 (MQ=255)
     gAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAATCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGc                                                                                                                                             >  8:217624/1‑149 (MQ=255)
      aaGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCg                                                                                                                                            <  3:20315/149‑1 (MQ=255)
        gACCGGCAAGCAGTCCGTGGTGATCACTGCTGGCGGCGACCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                                                          <  4:377388/149‑1 (MQ=255)
            ggCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCACGCCCGCTCCAACAAGGCCCCGGTAGGACCGGGC‑GTg                                                                                                                                     >  4:255315/1‑149 (MQ=255)
                 gCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCACCTGTTGGAGAGGGC‑TTGGCcgc                                                                                                                                >  4:284652/1‑149 (MQ=255)
                             gATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                                                          <  4:245109/128‑1 (MQ=255)
                             gATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCggg                                                                                                                                          >  3:245109/1‑128 (MQ=255)
                              aTCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCa                                                                                                                    >  2:31865/1‑148 (MQ=255)
                               tCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCaaa                                                                                                                  <  2:117321/149‑1 (MQ=255)
                                cACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTTGCCTGCTTCGCGGGCAAGCCCCCTCCGACAAAGCCCCCGTAAGCGCCGGG‑CTGCCCGCGAAATCCCCGCaaa                                                                                                                  >  2:193498/1‑148 (MQ=255)
                                 aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               >  8:204590/1‑119 (MQ=255)
                                 aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                               <  7:204590/119‑1 (MQ=255)
                                   tgctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAATCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACa                                                                                                              >  1:264223/1‑149 (MQ=255)
                                     ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAATCCCGCTCCTACAACGCCCCTGTAGGAGCGGGCTTTGCCCGCGAATGCCCCTCAAAGACaa                                                                                                             <  1:208520/149‑1 (MQ=255)
                                     ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGc                                                                                                                                    <  4:235075/125‑1 (MQ=255)
                                     ctgGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGc                                                                                                                                    >  3:235075/1‑125 (MQ=255)
                                               ccgccAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCGACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGCATGCCCCTCAAAGACAACCCGAGAGTCt                                                                                                  >  5:304633/1‑149 (MQ=255)
                                                   cAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTTGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTaa                                                                                             <  1:193498/149‑1 (MQ=255)
                                                    aGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATCCCCCTCAAAGACAACCCGAGATTCT‑GTAAt                                                                                            <  1:31865/149‑1 (MQ=255)
                                                                      aCTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCAt                                                                          <  3:284652/149‑1 (MQ=255)
                                                                                ttCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGCGCGGGC‑TTGCCCGCGCATGCCCCTCAAAGACAACCCGAGCTTCT‑GTAATGCCAGCCGCAATTCGCAGCACCCCCAcc                                                                >  6:276485/1‑149 (MQ=39)
                                                                                   ttACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGCCCTCAAGGCCCCGGTTGGAGCGGGC‑GAGCCCCCGGATGCCCCTCAAACACAAACCGCCAATGTGGTAAT‑CCAACCGCAATTACCCACACccccccccat                                                             >  8:31495/1‑147 (MQ=255)
                                                                                     aCGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCctg                                                             <  8:38028/147‑1 (MQ=255)
                                                                                       gCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGGAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCctgctgc                                                         >  6:353646/1‑149 (MQ=255)
                                                                                                  aaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCATACAAGCCCCCTGTAGTAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctggcc                                              <  1:356092/149‑1 (MQ=255)
                                                                                                       gcgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcac                                         <  2:264223/149‑1 (MQ=255)
                                                                                                          tCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagcag                                      <  1:85172/149‑2 (MQ=255)
                                                                                                             ccTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCGGTAGGAGCGGGC‑GTGGCCGCGGATGCCCCGCAAAGACAACCCGAGATTCT‑GTAATGCCAACCGCAATTGGCCTCACCCCCACCCAGCAGCTGGCCCTGGGCACCCCCACGGc                                   >  8:187594/1‑149 (MQ=37)
                                                                                                             ccTGCTTCGCGGGCAAGCCCCCTACTACAAGGCCCCTGTAGGAGCGGGC‑GTTCGCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTTATGCCATGCGCAATTAGCCTGACCCCCCCCCAGCTTCTGGCCCGGGGCAGCCCCCCTGc                                   >  8:73907/1‑149 (MQ=255)
                                                                                                              cTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcc                                  <  7:187594/149‑1 (MQ=255)
                                                                                                                        ggCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCCCCACTGCCCTTGCCAGcg                        >  4:95734/1‑149 (MQ=255)
                                                                                                                          cAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgac                      <  2:222584/149‑1 (MQ=255)
                                                                                                                               ccGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATTCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc                  >  3:210167/1‑148 (MQ=255)
                                                                                                                               ccGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc                  >  8:12696/1‑148 (MQ=255)
                                                                                                                               ccGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtc                  <  7:12696/148‑1 (MQ=255)
                                                                                                                                cGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                <  5:174338/149‑1 (MQ=255)
                                                                                                                                cGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCA‑CCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGCTCTCa               >  4:375466/1‑149 (MQ=255)
                                                                                                                                 gCTCCTACAGGGGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCGCCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa               >  7:249383/1‑149 (MQ=255)
                                                                                                                                 gCTCCTACAGGGGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GGAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa               >  3:362176/1‑149 (MQ=255)
                                                                                                                                 gCTCCTACAAGGCCCCTGTAGGACCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAAGGCCAACCGCAATTCCCATCACCCCCCCCCTCCTCATGGCCCTGGCCAGCACCAATCCACTTCCAACCGAGGATCTCa               >  4:31492/1‑149 (MQ=255)
                                                                                                                                   tCCTACAGGGGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCa                           >  5:189473/1‑135 (MQ=255)
                                                                                                                                   tCCTACAGGGGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCa                           <  6:189473/135‑1 (MQ=255)
                                                                                                                                     cTACAGGGGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GAAATGCCATCCGCAATTCGCATACCCCCCACCATCCTGCTGGCCATGGCCAGCCCCCCTGCCCTTGCAAGCCACGATCACAtacc           >  4:243512/1‑146 (MQ=25)
                                                                                                                                     cTACAAGGCCCCTGTAGGAGCGGTC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAc           <  6:304633/149‑1 (MQ=255)
                                                                                                                                     cTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                               <  4:115563/129‑1 (MQ=255)
                                                                                                                                     cTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                               <  4:120962/129‑1 (MQ=255)
                                                                                                                                     cTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                               >  3:120962/1‑129 (MQ=255)
                                                                                                                                     cTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                               >  3:115563/1‑129 (MQ=255)
                                                                                                                                      tACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                             >  2:317417/1‑130 (MQ=255)
                                                                                                                                      tACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                             <  1:317417/130‑1 (MQ=255)
                                                                                                                                        cAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATTCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                       >  6:119525/1‑134 (MQ=255)
                                                                                                                                        cAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                       <  5:119525/134‑1 (MQ=255)
                                                                                                                                            gCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                             <  4:140642/124‑1 (MQ=255)
                                                                                                                                            gCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                             >  3:140642/1‑124 (MQ=255)
                                                                                                                                              cccTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAc  >  4:22603/1‑149 (MQ=255)
                                                                                                                                              cccTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAc  >  4:277788/1‑149 (MQ=255)
                                                                                                                                              cccTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAc  >  5:159768/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAC  >  NC_002947/4049769‑4050059

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