Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,089,216	A→G	5.9%	Q26R (CAG→CGG)	fumC‑I →	class 2 fumarate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,089,216	0	A	G	5.9%	200.8 / 2.1	68	Q26R (CAG→CGG)	fumC‑I	class 2 fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base G (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

CGCCTTGCAGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTAT  >  NC_002947/1089072‑1089360                                                                                                                                                    |                                                                                                                                                 cGCCTTGCAGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGAcc                                                                                                                                              <  8:78006/149‑1 (MQ=255)    cTTGCAGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAg                                                                                                                                           >  4:84147/1‑149 (MQ=255)    cTTGCAGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAg                                                                                                                                           >  2:148290/1‑149 (MQ=255)         aGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGa                                                                                                                                                >  1:65016/1‑139 (MQ=255)         aGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGa                                                                                                                                                <  2:65016/139‑1 (MQ=255)           ttGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCgg                                                                                                                                    >  1:56814/1‑149 (MQ=255)              gCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATTAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTc                                                                                                                                  <  8:69848/148‑1 (MQ=255)                  aTTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACaa                                                                                                                             <  2:300536/149‑1 (MQ=255)                     gCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGTGGGGTACCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACaa                                                                                                                             >  1:115346/1‑149 (MQ=255)                                ccTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGCCAACTTCCCGATCAGCg                                                                                                               >  8:113678/1‑149 (MQ=255)                                      aTTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCg                                                                                                               <  8:45073/143‑1 (MQ=255)                                      aTTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCg                                                                                                               >  7:45073/1‑143 (MQ=255)                                          ctcCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTa                                                                                                      <  5:291315/148‑1 (MQ=255)                                                aTGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGc                                                                                               <  2:56814/149‑1 (MQ=255)                                                          gCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTcc                                                                                                                        >  5:22323/1‑114 (MQ=255)                                                          gCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTcc                                                                                                                        <  6:22323/114‑1 (MQ=255)                                                          gCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCAtt                                                                                     >  6:29815/1‑149 (MQ=255)                                                           cAAGGATTTGATGATGTGTGATTCCCGTATGGAGCGTGACAGGATGGGTGAGCTGCAGGTGCCGG‑‑‑GCCAGGCCCTCTAGGGCGCCCAGACCCAGCGGGCGTTCGACAACTTCCCGATCAGc                                                                                                                <  2:249495/121‑1 (MQ=255)                                                           cAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCGGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTc                                                                                    >  1:219767/1‑149 (MQ=255)                                                           cAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCGGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTc                                                                                    <  2:219767/149‑1 (MQ=255)                                                           cAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCGGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTc                                                                                    >  7:264774/1‑149 (MQ=255)                                                           cAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCGGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTc                                                                                    <  8:264774/149‑1 (MQ=255)                                                           cAAGGATTTGATGATGAGTGAAACCCGAATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGCCGCCCAGACCCAGCGCGCGGACGACAACTTCCCGCTAAGc                                                                                                                >  1:249495/1‑121 (MQ=255)                                                             aGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCg                                                                                                               <  2:66783/120‑1 (MQ=255)                                                             aGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCg                                                                                                               >  1:66783/1‑120 (MQ=255)                                                             aGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCgg                                                                                                >  6:24421/1‑135 (MQ=255)                                                             aGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCgg                                                                                                <  5:24421/135‑1 (MQ=255)                                                              ggATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTg                                                                                 <  3:84147/149‑1 (MQ=255)                                                                     atgatgaGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCttgt                                                                           <  5:147781/148‑1 (MQ=255)                                                                       gatgaGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGc                                                                        >  3:244574/1‑149 (MQ=255)                                                                          gaGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTgg                                                                     >  5:321604/1‑149 (MQ=255)                                                                           aGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGc                                                                    >  2:131480/1‑149 (MQ=255)                                                                             tGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCttgt                                                                           >  2:186488/1‑140 (MQ=255)                                                                             tGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCttgt                                                                           <  1:186488/140‑1 (MQ=255)                                                                               aTACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTg                                                                                 <  6:317262/132‑1 (MQ=255)                                                                               aTACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTg                                                                                 >  5:317262/1‑132 (MQ=255)                                                                               aTACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCaa                                                                 <  7:167531/148‑1 (MQ=255)                                                                               aTACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAg                                                                >  8:148401/1‑149 (MQ=255)                                                                                tACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAgg                                                               <  2:115346/149‑1 (MQ=255)                                                                                 aCCCGTATCGAGCGTGACAGCATCGGTGAGCTGCAGGAGCCGG‑‑‑CCCAGGCCCTCTACGGCCCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGc                                                              >  5:28877/1‑149 (MQ=255)                                                                                  cccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGcc                                                             >  6:147409/1‑149 (MQ=255)                                                                                   ccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCg                                                            >  6:312067/1‑149 (MQ=255)                                                                                   ccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCg                                                            >  2:229858/1‑149 (MQ=255)                                                                                   ccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCg                                                            >  8:73017/1‑149 (MQ=255)                                                                                   ccGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCg                                                            >  4:145552/1‑149 (MQ=255)                                                                                        tCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCg                                                       >  6:123076/1‑149 (MQ=255)                                                                                        tCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCg                                                       >  4:233777/1‑149 (MQ=255)                                                                                          gAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCgg                                                      <  5:29815/148‑1 (MQ=255)                                                                                              gTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggccaag                                                 <  7:148401/149‑1 (MQ=255)                                                                                                 aCAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggccaaggcc                                              <  1:148290/149‑1 (MQ=255)                                                                                                         ggTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCg                                       <  5:123076/148‑1 (MQ=255)                                                                                                             aGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACt                                   <  5:147409/148‑1 (MQ=255)                                                                                                               cTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTTTGCCGGCCCAGTTCATTCGTGcc                                                                               <  3:77946/102‑1 (MQ=255)                                                                                                                tGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGt                                                                                          <  5:4053/90‑1 (MQ=255)                                                                                                                tGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGt                                                                                          >  6:4053/1‑90 (MQ=255)                                                                                                                  cAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggccaaggc                                               >  6:124760/1‑131 (MQ=255)                                                                                                                  cAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCggccaaggc                                               <  5:124760/131‑1 (MQ=255)                                                                                                                  cAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAg                               <  6:321604/147‑1 (MQ=255)                                                                                                                   aGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCaa                            <  3:233777/149‑1 (MQ=255)                                                                                                                     gTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCGCAGCggccttggc                                               <  6:307084/128‑6 (MQ=255)                                                                                                                     gTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCTTGGCCGCAGCggccttggc                                               >  5:307084/1‑123 (MQ=255)                                                                                                                        ccGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGt                        <  1:131480/148‑1 (MQ=255)                                                                                                                           g‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGc                    >  1:34055/1‑149 (MQ=255)                                                                                                                                ccAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATTCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACCTCGAACTGGAGCAAATGTCAGCCg                  >  4:47590/1‑149 (MQ=255)                                                                                                                                      ccTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCttgt                                                                           <  6:76629/86‑1 (MQ=255)                                                                                                                                      ccTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCttgt                                                                           >  5:76629/1‑86 (MQ=255)                                                                                                                                          tACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAgg            >  6:23542/1‑145 (MQ=255)                                                                                                                                          tACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAgg            <  5:23542/145‑1 (MQ=255)                                                                                                                                                    aGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTAt  >  7:173448/1‑145 (MQ=255)                                                                                                                                                    aGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTAt  <  8:173448/145‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CGCCTTGCAGTTGGCCAATTGCCGTCTGCTGCCTTCGATTTCTCCCGATGAACTGATGCAAGGATTTGATGATGAGTGATACCCGTATCGAGCGTGACAGCATGGGTGAGCTGCAGGTGCCGG‑‑‑CCCAGGCCCTCTACGGCGCCCAGACCCAGCGCGCGGTCGACAACTTCCCGATCAGCGGCCAGCGTATGCCGGCCCAGTTCATTCGTGCCTTGTTGCTGGCCAAGGCCGCTGCGGCCAAGGCCAACGTCGAACTGGAGCAACTGTCTGCCGGGCAGGGCCAGGCTAT  >  NC_002947/1089072‑1089360

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