Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,089,682	T→C	6.4%	L250P (CTG→CCG)	kdgD →	5‑dehydro‑4‑deoxyglucarate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,089,682	0	T	C	6.4%	180.3 / 2.2	62	L250P (CTG→CCG)	kdgD	5‑dehydro‑4‑deoxyglucarate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (33/25);  new base C (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CCGCCTACAAGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATG  >  NC_002947/4089535‑4089824                                                                                                                                                    |                                                                                                                                                ccgccTACAAGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCGCGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTg                                                                                                                                                <  5:303906/149‑1 (MQ=255)   gccTACAAGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGa                                                                                                                                              <  3:155293/149‑1 (MQ=255)    ccTACAAGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACGACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGAc                                                                                                                                             <  2:254592/149‑1 (MQ=255)         aaGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATccg                                                                                                                                        >  2:117372/1‑149 (MQ=255)          aGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATccgc                                                                                                                                       <  2:261276/149‑1 (MQ=255)              ccTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATcagcaacc                                                                                                                                   >  3:23025/1‑149 (MQ=255)                        ccGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCgg                                                                                                                         <  2:37438/149‑1 (MQ=255)                                 tcctccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGt                                                                                                                >  7:75866/1‑149 (MQ=255)                                   ctccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGt                                                                                                                >  6:275971/1‑147 (MQ=255)                                   ctccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGt                                                                                                                >  2:228938/1‑147 (MQ=255)                                   ctccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTg                                                                                                               <  3:4462/148‑1 (MQ=255)                                   ctccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGa                                                                                                              <  7:60650/149‑1 (MQ=255)                                   ctccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCAACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGa                                                                                                              <  1:335314/149‑1 (MQ=255)                                    tccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATccgcaacc                                                                                                                                   >  2:154238/1‑127 (MQ=255)                                    tccGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATccgcaacc                                                                                                                                   <  1:154238/127‑1 (MQ=255)                                       gCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCAt                                                                                                          >  3:109888/1‑149 (MQ=255)                                           tgtTCAACTTCGTGTTGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGt                                                                                                       <  1:117372/148‑1 (MQ=255)                                             ttCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTc                                                                                                      >  3:82711/1‑147 (MQ=255)                                                aaCTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGc                                                                                                 >  3:219906/1‑149 (MQ=255)                                                 aCTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGcc                                                                                                >  5:245306/1‑149 (MQ=255)                                                   ttCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCg‑g                                                                                             >  5:132544/1‑149 (MQ=255)                                                      gTGCCGAAAACCGCCATGGACTTATACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAATCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTAAGCATCGTCAAGGCCG‑gcgc                                                                                          <  8:300613/149‑1 (MQ=255)                                                          cGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGa                                                                                      <  3:310032/149‑1 (MQ=255)                                                               aCCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGc                                                                                 <  5:275971/149‑1 (MQ=255)                                                                ccGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGcc                                                                                <  4:219906/149‑1 (MQ=255)                                                                    cATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGt                                                                             >  6:248361/1‑148 (MQ=255)                                                                           ttCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGt                                                                             <  7:309614/141‑1 (MQ=255)                                                                           ttCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGt                                                                             >  8:309614/1‑141 (MQ=255)                                                                                taaCGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCCGGACATCCGAAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCgctc                                                                >  5:297114/2‑146 (MQ=255)                                                                                cAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGcc                                                                >  7:29444/1‑149 (MQ=255)                                                                                 aaCGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGtt                                                                            >  7:62273/1‑136 (MQ=255)                                                                                 aaCGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGtt                                                                            <  8:62273/136‑1 (MQ=255)                                                                                 aaCGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGccc                                                               >  4:314885/1‑149 (MQ=255)                                                                                 aaCGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGccc                                                               >  5:256814/1‑149 (MQ=255)                                                                                      cATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGc                                                          <  1:228938/149‑1 (MQ=255)                                                                                            ccGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTAccc                                                    >  5:244568/1‑149 (MQ=255)                                                                                              gggacgacCACGCCACCTTGGCCAAGCGGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCg                                                  <  4:23025/147‑1 (MQ=255)                                                                                                   gacCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAAATGCCGGTTACGACGCCGGCCCGGt                                                            >  6:127101/1‑134 (MQ=255)                                                                                                       aCGCCACCGTTGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCCGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTgaccgac                                         <  2:265019/149‑1 (MQ=255)                                                                                                       aCGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCGCGCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTgaccga                                          <  1:88860/149‑1 (MQ=255)                                                                                                       aCGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTgaccgac                                         <  8:29444/149‑1 (MQ=255)                                                                                                           cACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCCGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCa                                     >  8:73835/1‑149 (MQ=255)                                                                                                              cGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATTCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCAccg                                  >  8:167520/1‑149 (MQ=255)                                                                                                                   ccAAGCTGATCGACGATTTCTTCCTGCCCTACCCGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGaa                             <  7:73835/149‑1 (MQ=255)                                                                                                                    cAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCAccgccg                               >  4:213466/1‑146 (MQ=255)                                                                                                                       gCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAagag                          >  4:66085/1‑148 (MQ=255)                                                                                                                       gCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAagag                          >  4:308681/1‑148 (MQ=255)                                                                                                                           aTCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACg                      >  5:308736/1‑148 (MQ=255)                                                                                                                              gacgaTTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAAt                  <  6:132544/149‑1 (MQ=255)                                                                                                                               acgaTTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATg                 <  6:256814/149‑1 (MQ=255)                                                                                                                                  aTTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCAccg                                  <  7:206919/129‑1 (MQ=255)                                                                                                                                  aTTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCAccg                                  >  8:206919/1‑129 (MQ=255)                                                                                                                                    ttcttcCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTcgc            >  3:67116/1‑149 (MQ=255)                                                                                                                                    ttcttcCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTcgc            >  5:16316/1‑149 (MQ=255)                                                                                                                                      cttcCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTcgccg          >  8:137139/1‑149 (MQ=255)                                                                                                                                           tGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCt      >  1:285261/1‑148 (MQ=255)                                                                                                                                           tGCCCTACCTGGACATCCGCAACCACAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCt      >  3:50543/1‑148 (MQ=255)                                                                                                                                            gCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGa    >  1:293349/1‑149 (MQ=255)                                                                                                                                             cccTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGa    <  3:314885/148‑1 (MQ=255)                                                                                                                                               cTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATg  >  1:160340/1‑148 (MQ=255)                                                                                                                                                   cTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAaga                           <  3:270043/119‑1 (MQ=255)                                                                                                                                                   cTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAaga                           >  4:270043/1‑119 (MQ=255)                                                                                                                                                    |                                                                                                                                                CCGCCTACAAGGCCCTGGGCGTACCGGTGTACTCCTCCGCTGTGTTCAACTTCGTGCCGAAAACCGCCATGGACTTCTACAACGCCATCGCCCGTGACGACCACGCCACCGTGGCCAAGCTGATCGACGATTTCTTCCTGCCCTACCTGGACATCCGCAACCGCAAGGCCGGCTACGCCGTGAGCATCGTCAAGGCCG‑GCGCCAGAATTGCCGGTTACGACGCCGGCCCGGTGCGTACCCCGCTGACCGACCTCACCGCCGAAGAGTACGAAATGCTCGCCGCGCTGATG  >  NC_002947/4089535‑4089824

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