Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,735,797	A→T	10.2%	L550Q (CTG→CAG)	hutU ←	urocanate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,735,797	0	A	T	10.2%	189.4 / 9.2	68	L550Q (CTG→CAG)	hutU	urocanate hydratase
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (3/4);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCC  >  NC_002947/5735658‑5735940                                                                                                                                            |                                                                                                                                                 ggTCACGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGc                                                                                                                                          <  5:217595/148‑1 (MQ=255) ggTCACGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGc                                                                                                                                          <  4:161890/148‑1 (MQ=255)   tCACGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTcc                                                                                                                                       >  1:216471/1‑149 (MQ=255)   tCACGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTcc                                                                                                                                       <  1:81701/149‑1 (MQ=255)      cGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTg                                                                                                                                    <  7:230025/149‑1 (MQ=255)       gTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGGCCTGGACCTGccgatgat                                                                                                                                   <  4:26039/149‑8 (MQ=255)            tactCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAg                                                                                                                              >  4:342/1‑149 (MQ=255)            tactCCAAATCCCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAg                                                                                                                              >  6:149531/1‑149 (MQ=255)               tCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑Tc                                                                                                                           >  4:316441/1‑148 (MQ=255)               tCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCAACCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑Tc                                                                                                                           >  4:109503/1‑148 (MQ=255)                ccAACTCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGa                                                                                                                         <  3:259896/149‑1 (MQ=255)                 cAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGAt                                                                                                                        >  2:3568/1‑149 (MQ=255)                  aaaTCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGAGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATg                                                                                                                       <  1:284454/149‑1 (MQ=255)                   actCGCGTGTAGCCGCAGGCCCGGGTGGGGGTTCTTGCGCGCTCCGTCTTGTTGTTCAGTATCTCCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCTTCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATgg                                                                                                                      <  8:62804/147‑1 (MQ=255)                               ccGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAACGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAAc                                                                                                          >  2:49642/1‑149 (MQ=255)                                     gCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCAt                                                                                                    >  7:190098/1‑149 (MQ=255)                                       ccGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCg                                                                                                  >  3:266884/1‑149 (MQ=255)                                        cGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGAGCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCgg                                                                                                 >  7:168690/1‑149 (MQ=255)                                             gtgGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAGCTCGATGGCGATGTCATAACCGGCATCGGCg                                                                                               <  1:29394/147‑1 (MQ=255)                                             gtgGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGAGCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGt                                                                                              <  2:216471/147‑1 (MQ=255)                                              tgGGTACCTGAGCGCTCCTTCTTGTTGGTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCTATGTCATAACCGGCATCGGCGTGa                                                                                            <  7:232099/148‑1 (MQ=255)                                                   aCCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCAt                                                                                       >  5:303940/1‑148 (MQ=255)                                                         gcgcTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCa                                                                                <  7:25597/149‑1 (MQ=255)                                                               cttcttGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGAt                                                                                                                  <  2:11831/109‑1 (MQ=255)                                                               cttcttGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGAt                                                                                                                  >  1:11831/1‑109 (MQ=255)                                                               cttcttGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCt                                                                          >  7:199029/1‑149 (MQ=255)                                                               cttcttGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCt                                                                          <  8:168690/149‑1 (MQ=255)                                                                        gttCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGtt                                                                  <  5:149531/148‑1 (MQ=255)                                                                           cAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGt                                                               <  3:342/148‑1 (MQ=255)                                                                           cAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGt                                                               <  1:3568/148‑1 (MQ=255)                                                                              tATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAg                                                            >  6:38166/1‑148 (MQ=255)                                                                              tATCACCCCTTGCAGGGTGCGATCCAAAGCGT‑GTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTa                                                          >  4:303267/1‑149 (MQ=255)                                                                                tCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCa                                                                                >  6:33969/1‑126 (MQ=255)                                                                                tCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCa                                                                                <  5:33969/126‑1 (MQ=255)                                                                                tCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTa                                                          <  4:266884/148‑1 (MQ=255)                                                                                  aCCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTAccc                                                       >  2:64293/1‑149 (MQ=255)                                                                                   ccccTTGCAGGGTGCGATCCAAGGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCa                                                             >  5:169650/1‑142 (MQ=255)                                                                                   ccccTTGCAGGGTGCGATCCAAGGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCa                                                             <  6:169650/142‑1 (MQ=255)                                                                                   ccccTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTa                                                          >  3:192712/1‑145 (MQ=255)                                                                                   ccccTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTa                                                          <  4:192712/145‑1 (MQ=255)                                                                                    cccTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCgg                                                     >  4:193364/1‑149 (MQ=255)                                                                                    cccTTGCAGGGTGCGATCCAAAACGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCgg                                                     <  1:49642/149‑1 (MQ=255)                                                                                                  gATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcg                                         >  7:271988/1‑147 (MQ=255)                                                                                                     ccAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCggcggc                                    >  4:127193/1‑149 (MQ=255)                                                                                                      cAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCggcggc                                    <  6:303940/148‑1 (MQ=255)                                                                                                       aaaGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGcc                                   >  3:158712/1‑148 (MQ=255)                                                                                                           cGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATc                              <  8:190098/149‑1 (MQ=255)                                                                                                            gtggtgATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCg                             >  6:78599/1‑149 (MQ=255)                                                                                                                tgATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTAc                         <  3:109503/149‑1 (MQ=255)                                                                                                                          gTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGAc               <  3:193364/149‑1 (MQ=255)                                                                                                                             atcatcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCggcggc                                    >  7:100775/1‑125 (MQ=255)                                                                                                                             atcatcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCggcggc                                    <  8:100775/125‑1 (MQ=255)                                                                                                                             atcatcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGAt            >  4:231021/1‑149 (MQ=255)                                                                                                                              tcatcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATg           <  5:38166/149‑1 (MQ=255)                                                                                                                              caaggagCAGGGCCTGGACCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcggcgg                                     <  6:103869/117‑1 (MQ=255)                                                                                                                              caaggagCAGGGCCTGGACCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcggcgg                                     >  5:103869/7‑123 (MQ=255)                                                                                                                              caaggagCAGGGCCTGGACCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcggcgg                                     <  8:147040/117‑1 (MQ=255)                                                                                                                              caaggagCAGGGCCTGGACCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcggcgg                                     >  7:147040/7‑123 (MQ=255)                                                                                                                                atcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcggcg                                      <  1:210480/120‑1 (MQ=255)                                                                                                                                atcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTcggcg                                      >  2:210480/1‑120 (MQ=255)                                                                                                                                atcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGa          <  3:303267/148‑1 (MQ=255)                                                                                                                                atcGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCCGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGa          <  2:308142/148‑1 (MQ=255)                                                                                                                                  gagCAGGGCCTGGACCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGt                                                               <  7:119079/91‑1 (MQ=255)                                                                                                                                  gagCAGGGCCTGGACCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGt                                                               >  8:119079/3‑93 (MQ=255)                                                                                                                                     cAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGc    <  5:78599/149‑1 (MQ=255)                                                                                                                                        gTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGccc  >  5:22565/1‑148 (MQ=255)                                                                                                                                         tCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGccc  >  6:138062/1‑147 (MQ=255)                                                                                                                                         tCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCAGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGAcc        >  3:220429/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                                 GGTCACGTTACTACTCCAAATCGCGTGGAGCCGCAGGCCCGCTTGTGGGTACCTGCGCGCTCCTTCTTGTTGTTCAGTATCACCCCTTGCAGGGTGCGATCCAAAGCGTGGTGATCAGCCGGTGATCATCGGCAGGTCCAGGCCCTGCTCCTTGGCGCAG‑TCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCC  >  NC_002947/5735658‑5735940

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