Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,153,506	T→A	8.7%	intergenic (‑308/‑61)	trmA ← / → btuB	tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase/vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,153,506	0	T	A	8.7%	139.2 / 5.1	46	intergenic (‑308/‑61)	trmA/btuB	tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase/vitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base A (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

GCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGCTACGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGT  >  CP009273/4153242‑4153770                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                         gCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGCTACGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCACGAAGACCTGCCGGCCAACGTCGCATCAGGTACTCATCATcgc                                                                                                                                                                                                                                                                 >  1:454411/1‑274 (MQ=255)                                  aTTTGTGAGCTACGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCAGCTCTTAGTATCTTAGATACCACTCAACGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCACGGAATATTAATGAAACATACGGCATCCTtcctct                                                                                                                                                                                                                               >  1:295092/1‑275 (MQ=255)                                                         tAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGGTTTACAATg                                                                                                                                                                                                           >  1:195694/1‑272 (MQ=255)                                                                     aCCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAACCTGGAGCTGTCGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTTTCATCATcgcg                                                                                                                                                                                                                                                                <  1:301872/207‑1 (MQ=255)                                                                          tAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCAGGCCAACGTCGCATCTGTTTCTCAACATCGCGTAACATTGATGAAACATTCGGCATCCATCTCCTATTGTGGATGCTTTTCACGGATTaaaaggacttatgtgct                                                                                                                                                                                       >  2:712133/1‑262 (MQ=255)                                                                               tCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCACGAAGACCTGCCGGCCAACGTCGCATCTGGATCTCATCATCGCGTCATATTGATGCAACCTGCGGCATCCTTCTTCTAtt                                                                                                                                                                                                                            >  1:320801/1‑233 (MQ=255)                                                                                                                                       cgcAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTgg                                                                                                                                                                                                                        <  1:256357/181‑1 (MQ=255)                                                                                                                                                      aGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTtctc                                                                                                                                                                                                                                                                         >  2:270803/1‑117 (MQ=255)                                                                                                                                                      aGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTCCAATGATTAAAAAAGCTTCCATGCGTACGGCGTGTTCCGTCACGACATTTTACGCTTGGGCCGCGGATACCAGCCGGGCTACTTGcgtcgtgccgggcca                                                                                                            >  2:261383/1‑265 (MQ=255)                                                                                                                                                          gCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCtcatca                                                                                                                                                                                                                                                                     <  1:270807/117‑1 (MQ=255)                                                                                                                                                              ttatTGCGGTATTCGGACACTGCCATTCGGTGTGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATTCTCTCGTCGTTACTGCTAACCGTTTTg                                                                                                    <  1:483172/272‑1 (MQ=255)                                                                                                                                                                 ttGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATCGATGAGAGCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAAc                                                                                                 <  1:712141/274‑1 (MQ=255)                                                                                                                                                                  tGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTACGTCACGGCATTTTCCGCTTGGGCACAGGCTACCAGCCCGGAAACGCCCGTCGTTCCTGCAAATGTTTTTGAACat                                                                                               >  2:799070/1‑274 (MQ=255)                                                                                                                                                                   gCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCtt                                                                                                                                                                                                                  >  1:168071/1‑159 (MQ=255)                                                                                                                                                                   gCGTTAGGCGGCCACTGCCATTAGCTGGGAACTCATCCTCTCTTTGTATCTTTGATTCACCCCCACGCCCGAAGACCTGCCGGCCAACCTCGACTCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCGTCTATTGTGGATGCTTTACAATGATGAAAAAAGCGTCGCTGCTGGCGGCGTGTTCCGTCACGGCATTTGCCGCTTGGGCACAGGATACCAGCCCGGATTCTCTCGTCGTTACTGCTATCCGTCTTGAACAg                                                                                               <  2:195691/274‑1 (MQ=255)                                                                                                                                                                       tATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACa                                                                                                                                                                                                              <  2:168071/159‑1 (MQ=255)                                                                                                                                                                             gACACTGCCATTCGGTGGGAAGTTATCATCTCTTTGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGACGCATCTGGTTCTCATCATCGCGGTATATTGATGAAACCTGCGGCATCCTTCTTCTATTTTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGACGCGCAGCAc                                                                                    <  1:799079/275‑1 (MQ=255)                                                                                                                                                                                   gCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCCTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGc                                                                               <  1:449135/274‑1 (MQ=255)                                                                                                                                                                                            tGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGc                                                                                                                                                                                                                                          >  2:411821/1‑110 (MQ=255)                                                                                                                                                                                            tGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGa                                                                                                   >  1:293912/1‑245 (MQ=255)                                                                                                                                                                                                aaGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATcc                                                                                                                                                                                                                                      <  1:411826/110‑1 (MQ=255)                                                                                                                                                                                                aaGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAg                                                                                               <  2:293908/245‑1 (MQ=255)                                                                                                                                                                                                          ctctTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATCCTCTCGTCGTTACTGCTAACCGTTTTGa                                                                                                   >  2:228052/1‑231 (MQ=255)                                                                                                                                                                                                          ctctTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGcc                                                                                             >  2:265178/1‑237 (MQ=255)                                                                                                                                                                                                              tAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAg                                                                                               <  1:228056/231‑1 (MQ=255)                                                                                                                                                                                                              tAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCgcgc                                                                                         <  1:265182/237‑1 (MQ=255)                                                                                                                                                                                                               aGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAACCCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGCACAGCCgcgc                                                                                         >  1:140383/1‑236 (MQ=255)                                                                                                                                                                                                               aGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCct                                                                                                                                                                                                                                     >  1:45193/1‑96 (MQ=255)                                                                                                                                                                                                                   tCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCcttctt                                                                                                                                                                                                                                 <  2:45193/96‑1 (MQ=255)                                                                                                                                                                                                                         aTACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCtcg                                                                                                                         >  2:50906/1‑194 (MQ=255)                                                                                                                                                                                                                            ccccTCCAAGCCCGAAGACCTGCTAGCCCACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGANACTCtcg                                                                                                                         >  1:573361/1‑191 (MQ=255)                                                                                                                                                                                                                             cccTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAatat                                                                                                                                                                                                                                                          >  1:63488/1‑61 (MQ=255)                                                                                                                                                                                                                             cccTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCgctgct                                                                                                                                                                                       >  1:400441/1‑128 (MQ=255)                                                                                                                                                                                                                             cccTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTcgtcgt                                                                                                                     <  1:50906/194‑1 (MQ=255)                                                                                                                                                                                                                               cTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTg                                                                                                                                                      >  1:684832/1‑159 (MQ=255)                                                                                                                                                                                                                                tCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTcgtcgt                                                                                                                     <  2:573356/191‑1 (MQ=255)                                                                                                                                                                                                                                 ccAAGCCCGAAGACCTGCCGTCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACg                                                                                                                                                                                   <  2:400436/128‑1 (MQ=255)                                                                                                                                                                                                                                 ccAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATtgat                                                                                                                                                                                                                                                      <  2:63488/61‑1 (MQ=255)                                                                                                                                                                                                                                 ccAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTAtt                                                                                                                                                                                                                            >  2:612763/1‑87 (MQ=255)                                                                                                                                                                                                                                   aaGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGca                                                                                                                                                  <  2:684825/159‑1 (MQ=255)                                                                                                                                                                                                                                    aGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCACATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTggg                                                                                                                                                    >  2:626266/1‑156 (MQ=255)                                                                                                                                                                                                                                     gCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCAACCTTCTTCTATTGTgg                                                                                                                                                                                                                        <  1:612768/87‑1 (MQ=255)                                                                                                                                                                                                                                        cGAAGACCTGCCGGCCAACGTCGCATCTGGTTCACATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGcaca                                                                                                                                                <  1:626273/156‑1 (MQ=255)                                                                                                                                                                                                                                          aaGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTCCCCCAACCACCGTGGTGACCCGTCCGTATACCGACCGCNGTCCGTCGTCCTcc                       >  2:600933/1‑274 (MQ=255)                                                                                                                                                                                                                                                 gCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGAACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCCCCCACCACCGTTGTGACCCGTCAGGATATCAACCGCTGGCAGTCGCCCTCGGTCact                 >  1:811696/1‑272 (MQ=255)                                                                                                                                                                                                                                                          cGTCGCATCTGGTTCACATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCt                                                                                                                                                        >  2:535295/1‑130 (MQ=255)                                                                                                                                                                                                                                                              gCATCTGGTTCACATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTTACGGCGTGTTCCGTCACGGCATTTTCCGCTTggg                                                                                                                                                    <  1:535300/130‑1 (MQ=255)                                                                                                                                                                                                                                                               cATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTNCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGt  <  1:600938/275‑1 (MQ=255)                                                                                                                                                                                                                                                                        ctcATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGa                                                                                                   >  1:257545/1‑169 (MQ=255)                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                         GCTCATCAGATGTCCAGATCTTGATGAATTCCTATTTGTGAGCTACGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGT  >  CP009273/4153242‑4153770

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