Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,883	+C	coding (667/687 nt)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,883	1	.	C	97.2%	103.8 / ‑2.7	36	coding (667/687 nt)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (18/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑C‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCT  >  NC_000913/3815636‑3816122                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc‑cctc                                                                                                                                                                                                                                                <  1:172078/251‑4 (MQ=255) cgcgCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc‑cctc                                                                                                                                                                                                                                                <  1:73145/251‑4 (MQ=255)           aTTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGc                                                                                                                                                                                                                                      >  2:273252/1‑251 (MQ=255)               aaCAGTCCTGCGGTGAAGAAATAGTGGCTTTTGCGCCCGGACTTCACCGTAAACTCGCCAATCGTTAACCACGCCTTGCTAATCTCAAAGTCAATAAACTGTCGCGTATATGGTTTCATGCCTTCGCTCCCCATCTTCCTTTCCTACTGCCCAAAAAAAGGCGTCACATCAGTCGCCTTAAAAACCAGTGTTCCAGCGCCGCCTTCTGCGTCGCTACAATGGACACGATTCCCAC‑TCGGGCCAGATCcaac                                                                                                                                                                                                                                  <  2:76760/251‑1 (MQ=255)                                                              cgtgaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                                                                                                   >  1:358508/5‑251 (MQ=255)                                                                              tttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCATTCGCTCCTCATCTTACTTTTCTACAGCCAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGCTTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                                   >  1:171959/1‑251 (MQ=255)                                                                                    acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCACCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTACAATTCAATGATGCGcccc                                                                                                                                                             >  2:261160/1‑250 (MQ=255)                                                                                               aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                                                                                                                                       <  2:325593/182‑1 (MQ=255)                                                                                               aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAg                                                                                                                                                                                                                       >  1:325593/1‑182 (MQ=255)                                                                                                aGCGCAAATTCAAGAAACCGGAGCTGATATGGTTTCAGGCCTTCGCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCAACAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                                                                 <  2:365351/251‑1 (MQ=255)                                                                                                  cgcAAATTCAATAAACTGGCGCTGATAGGGTTGCATGCCTGCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTGCGGTCATCAc                                                                                                                                               <  2:130000/251‑1 (MQ=255)                                                                                                          cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTAATCACTACGTTc                                                                                                                                        <  1:273252/250‑1 (MQ=255)                                                                                                           aaTAAACTGGCGCGTAGATAGTTACATGCCTTCGCTCCTCATCGTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTATAAAACAGGTTACCCGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTTGCCTTCTGCCGTCCCCTGCACGTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCAt                                                                                                                                      <  2:151129/251‑1 (MQ=255)                                                                                                                    gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGtctc                                                                                                                             >  1:38518/1‑251 (MQ=255)                                                                                                                                   cATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTc                                                                                                              <  1:73552/251‑1 (MQ=255)                                                                                                                                       ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACg                                                                                                         >  1:325715/1‑251 (MQ=255)                                                                                                                                        cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCGGGGTCTTCA‑ACgg                                                                                                        >  2:234991/1‑250 (MQ=255)                                                                                                                                            gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAACAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTc                                                                                                    <  1:261160/251‑1 (MQ=255)                                                                                                                                                                 cAGACAAAAAAAAGGCGACTCACCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGcc                                                                               <  2:325715/251‑1 (MQ=255)                                                                                                                                                                             aGGCGAATCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTACATGTCGGTCTCTGCGGCAGAGTCTTCACAGGT‑TACCAGACCGCAAACCGCTTCGCCGTTCAAAATTcg                                                                   >  2:263583/1‑250 (MQ=255)                                                                                                                                                                                    tcatcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCTTCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGa                                                             >  2:153578/1‑250 (MQ=255)                                                                                                                                                                                       ccaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACt                                                         <  1:208975/250‑1 (MQ=255)                                                                                                                                                                                         aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                                                                  <  1:319035/161‑1 (MQ=255)                                                                                                                                                                                         aGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcat                                                                                                                                                  >  2:319035/1‑161 (MQ=255)                                                                                                                                                                                          gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCACCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTg                                                        >  1:49479/1‑249 (MQ=255)                                                                                                                                                                                             gCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCt                                                   <  1:133100/251‑1 (MQ=255)                                                                                                                                                                                              ccTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTa                                                  <  1:42400/251‑1 (MQ=255)                                                                                                                                                                                                  aaaaaTCAGTTTGCCAGCGCCGCCTTCTGAGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCAt                                              >  1:288129/1‑251 (MQ=255)                                                                                                                                                                                                           tttGCCAGAACCGCCTTCTGCGTCGCTAGAATGGATTCGATTCCC‑CCTCGGGCCCGAGCAAACAAGCTGAGAAGCGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCGCCTGCACTTCAATGTTGCACCCGTCTTCGGTCATCACTATGTTCATGTCGATATCTGCGGCCGAGGCGTCA‑ACGTATTCCATATCGCTATCCGCTTCGCCGCTCACAATTCCGACAGAAACTGCGGCTAccgtcgttttgct                                     >  2:37430/1‑240 (MQ=255)                                                                                                                                                                                                                cAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGAt                                <  2:358508/251‑1 (MQ=255)                                                                                                                                                                                                                    gccgccTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTAACCATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTggg                            >  2:403103/1‑250 (MQ=255)                                                                                                                                                                                                                     ccgccTTCTGCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCCACCATCCCTTTCATCGGATTGGtt                           >  2:186316/1‑251 (MQ=255)                                                                                                                                                                                                                              gCGTCGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCtt                   <  1:165778/250‑1 (MQ=255)                                                                                                                                                                                                                                  cGCTACAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTAGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGt              >  2:258485/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCAATGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTGCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGttt            <  1:403103/249‑1 (MQ=255)                                                                                                                                                                                                                                          tGGATTCGATTCCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTccacc      >  2:341129/1‑251 (MQ=255)                                                                                                                                                                                                                                            gattcgattcCC‑CCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAg    <  2:38518/251‑1 (MQ=255)                                                                                                                                                                                                                                               tcgattcCC‑CCGCGGGCCAGATCCATCAATATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCt  <  1:23095/250‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 CGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCC‑C‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCA‑ACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCT  >  NC_000913/3815636‑3816122

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