Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	740,893	G→A	A463T (GCA→ACA)	phr →	deoxyribodipyrimidine photolyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	740,893	0	G	A	94.7%	128.7 / NA	38	A463T (GCA→ACA)	phr	deoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/23);  minor base C (2/0);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGGCCGCTTCAACCGGAACCGATGCAGCGCCGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGGCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCA  >  NC_000913/740660‑741133                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 gggCCGCTTCAACCGGAACCGATGCAGCGCCGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAggcggc                                                                                                                                                                                                                                  <  1:345229/250‑1 (MQ=255)          cAACCGGAACCGATGCAGCGCCGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAgggg                                                                                                                                                                                                                        <  1:83922/251‑1 (MQ=255)           aaCCGGAACCGATGCAGCGCCGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAgggg                                                                                                                                                                                                                        >  1:238858/1‑250 (MQ=255)                            cgcCGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAAGAACAGACGCGTCAg                                                                                                                                                                                                      <  1:308814/251‑1 (MQ=255)                              ccGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑CGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGc                                                                                                                                                                                                    >  1:404777/1‑251 (MQ=255)                                   ttttcgtattttcAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCg                                                                                                                                                                                               >  1:21308/1‑251 (MQ=255)                                     ttcgtattttcAACCCGACAACCCAGGGCGAGCAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGGGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCg                                                                                                                                                                                               >  2:283929/1‑249 (MQ=255)                                       cgtattttcAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGc                                                                                                                                                                                            <  1:86080/250‑1 (MQ=255)                                        gtagttgcAACCCCACAACCCAGGGGGAGAAATTT‑AGAATGGAGGCCGATTTTATCCGTCAGTGGCTACCTGAACGGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTTACGCTGGATTAGCCGCAACCGATAGGCGAGCACAAAGAAGCGAGAGTACAATCGTTGACAGCGTATGAGGCGGCGGGGAAGGGGAAATAACAGACGCGTCAGGCAATCGATccc                                                                                                                                                                                          <  2:161782/244‑1 (MQ=255)                                            tttcAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGa                                                                                                                                                                                       <  2:169732/250‑1 (MQ=255)                                                aaCCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCg                                                                                                                                                                                  >  2:391010/1‑251 (MQ=255)                                                   ccgccctcccAGGGCAAGAAATTGGATCAGG‑AGGCCGAGTTTTTCCGCCAGTGGCTACCGGAACGGCGCTATGTGCCAGGGAAAGGGGTGCATGAGCCGTGGATGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAAACGATAGTCTAGCACAAAGAAGCGAGAGTACAAACTTTGACAGCGTATGAGGCGGGGCGTAAGGGGAAATAACAGACGCTTCAGGCAATCGAGCCCAGATGCCGGAt                                                                                                                                                                               <  2:241341/244‑1 (MQ=255)                                                         aCCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATgcggcg                                                                                                                                                                         <  1:410448/251‑1 (MQ=255)                                                            cAGGGCGAGAAATTTGATCATG‑AGGGCGGGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGa                                                                                                                                                                      <  2:376911/251‑1 (MQ=255)                                                            cAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTg                                                                                                                                                                       >  1:207851/1‑250 (MQ=255)                                                                     aaaTTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCAGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                             >  1:260474/1‑251 (MQ=255)                                                                          tGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAg                                                                                                                                                         <  1:299570/250‑1 (MQ=255)                                                                                     ggCGAGTGTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTGTTCCAGCCTACAAAttt                                                                                                                                              <  2:233450/251‑1 (MQ=255)                                                                                             tATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAAt                                                                                                                                      >  2:92548/1‑251 (MQ=255)                                                                                                   ccAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGCGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATa                                                                                                                                <  2:326864/251‑1 (MQ=255)                                                                                                                  aCTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATAc                                                                                                                 <  1:6335/251‑1 (MQ=255)                                                                                                                        cgATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGa                                                                                                                      >  1:285534/1‑240 (MQ=255)                                                                                                                        cgATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGa                                                                                                                      <  2:285534/240‑1 (MQ=255)                                                                                                                         gATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGcc                                                                                                          >  2:300064/1‑251 (MQ=255)                                                                                                                            gtgCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTg                                                                                                        <  1:41682/250‑1 (MQ=255)                                                                                                                                   ggAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGaaa                                                                                                                    <  2:180552/231‑1 (MQ=255)                                                                                                                                   ggAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGaaa                                                                                                                    >  1:180552/1‑231 (MQ=255)                                                                                                                                        gtggtgCATGAGCCGTGGAAGTGGGCGCAGAAAGAAGGCGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTc                                                                                            >  1:367438/1‑250 (MQ=255)                                                                                                                                               aTGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCAGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAAAAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCg                                                                                    <  2:260474/251‑1 (MQ=255)                                                                                                                                                 gAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGGAGGCCTGATAAGACGCGTCAAGCGTCGc                                                                                   <  1:232900/250‑1 (MQ=255)                                                                                                                                                  aGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGCGAGTACAAACGTTGCCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACACAATTGTGCAAATTCAATAAATTGCCGGAAATACGTACGCCTGATAAAACGCGTCAAGCGTCGCAt                                                                                 >  1:194868/1‑251 (MQ=39)                                                                                                                                                           aGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCg                                                                        <  1:49461/251‑1 (MQ=255)                                                                                                                                                                                    gCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTccagcg                                               <  1:92548/251‑1 (MQ=255)                                                                                                                                                                                                  aaCCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAGGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGAGGACGGATGCGGCGTGAACGGCTTATCCAGCCTACAAATTGGTGCAAATTCAATAAATTGCAGGAAAGACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAATACtggagggccttcgctcggttt                                 <  2:61648/251‑21 (MQ=255)                                                                                                                                                                                                      gATAGTCGAGCACAAAGAAGCGAGCGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCACATTCAATGAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGGCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTg                              >  2:359993/1‑250 (MQ=255)                                                                                                                                                                                                                  caAAGAAGCGAGAGTACAAACGTTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAAGTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGcca                 <  2:367438/251‑2 (MQ=255)                                                                                                                                                                                                                             gagTACAAACGGTGACAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGCGCCCAGATGCCGGATGCGGCTTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCaa      <  2:364995/251‑1 (MQ=255)                                                                                                                                                                                                                                 aCAAACGTTGCCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTACGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCCAatcc  >  2:377369/1‑250 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 GGGCCGCTTCAACCGGAACCGATGCAGCGCCGTATTTTCGTATTTTCAACCCGACAACCCAGGGCGAGAAATTTGATCATG‑AGGGCGAGTTTATCCGCCAGTGGCTACCGGAACTGCGCGATGTGCCAGGGAAAGTGGTGCATGAGCCGTGGAAGTGGGCGCAGAAAGCAGGTGTGACGCTGGATTATCCGCAACCGATAGTCGAGCACAAAGAAGCGAGAGTACAAACGTTGGCAGCGTATGAGGCGGCGCGGAAGGGGAAATAACAGACGCGTCAGGCAATCGAGCCCAGATGCCGGATGCGGCGTGAACGCCTTATCCAGCCTACAAATTTGTGCAAATTCAATAAATTGCAGGAAATACGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCATCGGTGCTGGCCTATTAAGACTCCAGCGCCAGCGCGCGGTTTCTGAATTTCAGCGCCTGATACAGCCAAATCA  >  NC_000913/740660‑741133

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑