Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	2,251,616	A→G	R94R (CGT→CGC)	LY180_10745 ←	major tail tube protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	2,251,616	0	A	G	93.5%	41.4 / ‑5.0	32	R94R (CGT→CGC)	LY180_10745	major tail tube protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (12/17);  minor base T (0/2);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

AGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTCAGGC  >  CP006584/2251410‑2251906                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                   gcgagggcgaggtcGTTCTCCTTCTGAATCAGGGTGATGGTTTCGTGTTCGGACAGTTCCTTACCTGCCATCGGTGGCCGGAAATAGGTGTAGGCCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGGTTGCCCTCGCCGGTGTCGTTTTCTTTGTGGCGTCCACGCAGGACCACCTCGGCGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGa                                                                                                                                                                                                         <  2:334471/285‑1 (MQ=11)     ggTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGa                                                                                                                                                                                                     >  1:421388/1‑298 (MQ=14)                ttCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTg                                                                                                                                                                                          <  1:543836/298‑1 (MQ=14)                                tGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGtcatc                                                                                                                                                                          >  1:406808/1‑298 (MQ=14)                                            cGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTc                                                                                                                                                              <  1:329868/298‑1 (MQ=14)                                             gATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCt                                                                                                                                                             >  1:380144/1‑298 (MQ=14)                                                  cccccAGCCCCTTTCCGTCCCTCGTCCGCCGGAAATCGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCGGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGCCGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCgg                                                                                                                                                        <  2:391041/295‑1 (MQ=11)                                                   gcccAGCTCCTTCCCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCTGAGGTCTACCGGt                                                                                                                                                       <  1:588155/296‑1 (MQ=14)                                                       aGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCAAGAGAACCGAATCCGGGAAGCCACCGCGCGACCACCCCATTGCCAGCCCACCAACACACAGACCGCGGTCTACCGGtcgcc                                                                                                                                                   >  2:462754/1‑294 (MQ=11)                                                        gCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGCATCGTCCCGCTGGTAAGAGCCTGCAAAACGAATCCGCCCAGCATCCACACCGGGTCCGGCGGAAAGCCCCCAGCAACACGAATCCGGGAAGCCACCGAGCGACCAACCCATTGACAGCGCACCGCCATCCAGCCCACGGCCTACCGCTccgcg                                                                                                                                                  >  2:562264/1‑297 (MQ=11)                                                          tCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGc                                                                                                                                                <  1:32773/298‑1 (MQ=14)                                                                 cGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCCCCCAGATAACCGAACCCGGGAAGCCACCGAGCGACCCCTCCATTGACAGCGCCTCGTCATCCAGGCCGAGGGCTACCGGTCCGCTGCCGTTCac                                                                                                                                         >  2:478448/1‑297 (MQ=255)                                                                   tcgaagctgtgCCGGAAAGATGTGATGACTATGGGGATTTTCGATTCGGTGTTTTTTCCCTGGTTGCCCTTGCCGGTGTCGATGTCCTGGTGACGTCCTCGCATGACCACCTCGTCGGCCACCTTTTGTGAGCTTTCGTCTCGCCGGTAAGGGCCTGCAAAAAGAAGCGGCAATGCAACCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAAGCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATcc                                                                                                                                       <  2:345461/288‑1 (MQ=11)                                                                   cccgctgccggcgggAAATGGGTGCGGTCCACGGGGATTTTCGTTTCGGGTTTTTCTCCCTGTTTCCCCCTGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACGTGGACGGCAACCGTTTCGACGGTAGCGTCGCGCTGTGAAGAGCCTGCAAAACGAATCTGCACAGCAGCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATcc                                                                                                                                       <  2:421388/285‑1 (MQ=11)                                                                                                   ggAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTc                                                                                                       <  1:386735/298‑1 (MQ=14)                                                                                                     gggttcgcGACTTGGTTTGTTCTCCGTGTTTGCCCTCGCCGTTGTCATTTTCTTTCTGGCGTCCACGCATGACCACCTGGA‑GCGCACCGGTTCGCCGGTCTCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGt                                                                                                    <  2:89026/291‑1 (MQ=11)                                                                                                       tttttCGGCTCGGTGTCTCTCCCCTGTTTGGCCTCGCCGGGGTCGATTTCTGTCTGACGTCCGCGCATGACCACCTCGACGGCCACCGTGTCGCCGGTATCGTCGCGCTGGGAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGtt                                                                                                   <  2:109115/297‑1 (MQ=11)                                                                                                            cGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGc                                                                                              <  1:632208/298‑1 (MQ=14)                                                                                                             gaccctgtgtgtACTCCCTGTTGGCGCTCGCGCGTGTCGATTGCTGTGTGACGTGCACGCCGGGCCGCCCAGAGGGACACCGGGTCGCCGGTAGCGTCGCGCGGGTAAGAGGCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGAGAACCGAAGCCGGGAAGCCACCGGGCGGCCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCa                                                                                             <  2:656801/292‑1 (MQ=11)                                                                                                               cTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGc                                                                                           >  1:589356/1‑298 (MQ=14)                                                                                                                              cTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCa                                                                            >  1:312210/1‑298 (MQ=14)                                                                                                                                  ttGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGc                                                                        <  1:200530/298‑1 (MQ=14)                                                                                                                                   ggaccTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGcc                                                                       <  1:462754/295‑1 (MQ=14)                                                                                                                                      ccacGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCCGCAAAACGAATCGGCCCAGCATCCACACCGGTTGCGGCGTCAAACTCCCAGATAACCGAATCCGGGAAGCCGCCGAGCGACCAGTCCATTGCCGGCGCAGCGTCACCAAGGCCGAGGTCTACCGGTCCGCTCCCCATCAACCCTGTCCCGCCATTGCTGTCGAACTCACGGGTCCGTCTCGGCAGCCTCTCGGACCTTGCAACCCCCTg                                                                    >  2:32305/4‑298 (MQ=11)                                                                                                                                            cGGTGTCGATTTCTTTCTGAGGCCCCGGGATCGCCACCTGGACGGCGCCGGTTTCGCCGGGATCGTCGCGGTGGTAAGAGCCGGCAACACGGATCGGCACAGCGTCCACACCGGTTGCGGCGTAAAGCTCCCAGATAGCCGAATCCGGGAAGCCACCGAGCGACCACTCCATGGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCt                                                              <  2:389253/298‑1 (MQ=11)                                                                                                                                                    aTTTCTTTCTGACGTCCACGCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGt                                                      >  1:506698/1‑298 (MQ=14)                                                                                                                                                                       gCATGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATCGCTGTAGCCGTTCAGAAAACCGTTCATTaa                                   >  2:296373/1‑298 (MQ=9)                                                                                                                                                                         aTGACCACCTCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTCCGGGTCCGTTTTGGCAGCGTGACGGACTTTGCCACGCCCTGATCGCTGTAGCCGTTCAGAAAGGCGTTCGTTAACt                                 >  2:272367/1‑298 (MQ=9)                                                                                                                                                                                 ctCGACGGCCACCGTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTg                         >  1:197446/2‑298 (MQ=14)                                                                                                                                                                                     ccGGCCACCGTTGCGCCGGTAGCGGCGGGCGGGTAAGGGCCGGAAAAACGAATCGGAACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGcgcg                     <  2:197446/297‑1 (MQ=11)                                                                                                                                                                                              gTTTCGCCGGTATCGTCGCGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCa            <  1:607458/298‑1 (MQ=14)                                                                                                                                                                                                        tagcGTCTCGCTGGTAAGCGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGATCGACCACGCCATTGACAGCGCATCGTCGTCAAGGCCGAGGTCTACCGGGGCGCTGCCGTTCATCCCCTCACCGCGATAGTTGTCGAGCTTACTGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTCAGGc  <  1:183819/295‑1 (MQ=11)                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                   AGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTTTCGACTCGGTGTCTTCTCCCTGTTTGCCCTCGCCGGTGTCGATTTCTTTCTGACGTCCACGCATGACCACTTCGACGGCCACCGTTTCGCCGGTATCGTCACGCTGGTAAGAGCCTGCAAAACGAATCGGCACAGCATCCACACCGGTTGCGGCGTAAAGCTCCCAGATAACCGAATCCGGGAAGCCACCGAGCGACCACTCCATTGACAGCGCATCGTCATCAAGGCCGAGGTCTACCGGTGCGCTGCCGTTCATCCCCGCACCGCGATAGTTTTCGAGCTTACGGGTCAGTTTTGGCAGCGTGACGGACTTTGCAACGCCCTGATAGCTGTAGCCGTTCAGAAAGACGTTCATTAACTTGAGTTTGCGCGGCATTGCCATCGGTCAGGC  >  CP006584/2251410‑2251906

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