Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,049,910	C→G	15.1%	intergenic (+43/‑58)	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_002947	4,049,910	0	C	G	15.1%	140.7 / 17.3	59	intergenic (+43/‑58)	quiA/oprB‑III	quinate dehydrogenase/carbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base G (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 = 2.59e-01

CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAC  >  NC_002947/4049762‑4050050                                                                                                                                                     |                                                                                                                                               cTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGc                                                                                                                                                >  4:20315/1‑149 (MQ=255)    cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGccc                                                                                                                                              <  3:358553/148‑1 (MQ=255)    cGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGccc                                                                                                                                              <  4:48851/148‑1 (MQ=255)     gTCTCGCCCAACACCTGCAAGCAGTACGTCGTGATAACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGgcctt                                                                                                                                            <  8:313435/149‑3 (MQ=38)        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                                                                                                                                      <  7:204590/119‑1 (MQ=255)                                         aCTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGa                                                                                                                                      >  8:204590/1‑119 (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                  <  6:189473/135‑1 (MQ=255)                                                                                                                                           tCCTACAGGGGCCTTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCa                  >  5:189473/1‑135 (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                      >  3:115563/1‑129 (MQ=255)                                                                                                                                             cTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                      >  3:120962/1‑129 (MQ=255)                                                                                                                                             cTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                      <  4:120962/129‑1 (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                    >  3:140642/1‑124 (MQ=255)                                                                                                                                                    gCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                    <  4:140642/124‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               CTACGTCTCGCCGAAGACCGGCAAGCAGTACGTGGTGATCACTGCTGGCGGCGCCCGCCAGTCGACTGACCGTGGCGACTATGTGATTTCTTACGCTCTGCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCT‑GTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAC  >  NC_002947/4049762‑4050050

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: ‑