Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	793,985	A→G	N394D (AAC→GAC)	LY180_03720 →	rhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	793,985	0	A	G	98.1%	156.1 / ‑4.3	53	N394D (AAC→GAC)	LY180_03720	rhsA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (31/21);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAAAAAGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGG  >  CP006584/793778‑794267                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                            gtgtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTAGGAAAAAGACCGCATCACCATCCCCGACAGCCTGAACCGCCGCGAAGTGCTGCACACGCAGGGCGACGCCGGGCGAACGCGGGGGGTGAAAAAGGAACACCCGGACGGCAGCgacccgccgacgc                                                                                                                                                                                                   >  2:695895/1‑286 (MQ=1)     aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGCAAACGGAACACGCGGACGGCAGCGTCACAAAGAGTCAGtt                                                                                                                                                                                               >  1:586193/1‑298 (MQ=9)     aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGtt                                                                                                                                                                                               <  1:774883/298‑1 (MQ=14)           gccgcaATAAACAGGGTCGAAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTAGCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGc                                                                                                                                                                                         <  2:796663/295‑1 (MQ=11)              gcaATCCACGGGGGGGGAGCTTTACTTACGATGAGAAATACCGGGGCCGGATGGTGGCGCCCCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGt                                                                                                                                                                                      <  2:68998/298‑1 (MQ=11)                       aGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCt                                                                                                                                                                             <  1:563925/298‑1 (MQ=17)                              cAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGGGGTGAAAAAGGAACACGCGGACGGCAGCGTCCCACAGAGGCCGTTTAACGCCGTGGGCAGGCCCAGGGca                                                                                                                                                                      >  2:391225/1‑298 (MQ=11)                                gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGcaca                                                                                                                                                                    >  1:660064/1‑298 (MQ=21)                                 cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAg                                                                                                                                                                   >  1:636381/1‑298 (MQ=21)                                    tacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACg                                                                                                                                                               >  1:153316/1‑298 (MQ=21)                                       ttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCCTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGCAGTGCTGCCCACGCAGGGCGAAGCCGGGCTGAAGCGGGGGGTGAAAAAGGAACACGCGCACGGCAGCCTCCCACAGAGGCCACTTGACGCCGAGGGGAGGCGCCGGGGCCAG‑ACgcgt                                                                                                                                                            >  2:133523/1‑295 (MQ=11)                                        tacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATg                                                                                                                                                           <  1:232706/298‑1 (MQ=21)                                          cGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGAGGcc                                                                                                                                                         >  1:137181/1‑298 (MQ=19)                                           gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCg                                                                                                                                                        >  1:601602/1‑298 (MQ=21)                                                  aaTACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGCGGGCAGGCTCAGGGCACAG‑ACGGATGCCgcaggcag                                                                                                                                                 <  1:488461/298‑1 (MQ=255)                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCCGGGCGCAGCCGGGCCGAAGCGGGTGGTGAAAAAGGAACACGCGGACCGCAGCGTCACCCACCGTCAGTTTGACGCCGTGGGCCGGCGCAGGGCACAG‑ACGGCGGCCGCCGGACGGaca                                                                                                                                             >  2:179584/1‑298 (MQ=11)                                                       cGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGacaa                                                                                                                                            <  1:77205/298‑1 (MQ=25)                                                              ggATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAGGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTa                                                                                                                                     >  1:602313/1‑298 (MQ=255)                                                               gATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTAc                                                                                                                                    >  1:313476/1‑298 (MQ=255)                                                                              tCACACGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATgtggtg                                                                                                                     <  1:577030/298‑1 (MQ=255)                                                                                 cacGGGCCGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACg                                                                                                                  >  1:709050/1‑298 (MQ=255)                                                                                        cGACCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCtca                                                                                                           <  1:139446/298‑1 (MQ=255)                                                                                        cGAACGGAATTCCGTTACCGTTACGGCGGGGACGGGCGGGGGACAGAACAGCCAAACCCGGCAGGCTTAAGCTACACGTATCAGGAGGGGAAAGACCGCATCACCATCACCGACAGCCTGGACGGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCtca                                                                                                           <  2:153316/298‑1 (MQ=11)                                                                                           ccGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCtcatca                                                                                                        >  1:300545/1‑298 (MQ=255)                                                                                              gAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGGGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGGCAGTTTGAAGCCGTGGGCAGGCCCAGGGCACAGCACGCGT‑CCGCAGGCAGGACAACCAAGGACAGCCCGGAAGTGGTGACCGGCCccaaccccg                                                                                                     >  2:58959/1‑292 (MQ=11)                                                                                                       tACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGGCTCATCACGCGCATAACCa                                                                                            >  1:427249/1‑298 (MQ=255)                                                                                                          gggtgcggaggCGACGGGGGGGGGAAAGAACAGCTAAACCCGGCAGGGTTAAGCTACACGTAGCAGGATGGGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGAACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGc                                                                                         <  2:307301/288‑1 (MQ=14)                                                                                                              agaggagggcgggggggggGAAAGACCAACGAAACCCCGCAGGTTTAAGCGACACGGATCAGTAGGGGAAAGACCGCATCGCCATCGCCGAGAGCCGGGACGGCGGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGGGAAAAAGGAACACGCGGACGGCAGCGGCACACAGAGTCAGTTGGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGa                                                                                     <  2:313476/287‑1 (MQ=14)                                                                                                                 aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATgg                                                                                  >  1:369193/1‑298 (MQ=255)                                                                                                                   aGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCa                                                                                >  1:272122/1‑298 (MQ=255)                                                                                                                         gggcggtTTAAAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCAt                                                                          <  1:659355/298‑1 (MQ=255)                                                                                                                          ggcgggTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGATCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATc                                                                         >  1:642186/1‑298 (MQ=255)                                                                                                                            cgggTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCgg                                                                       <  1:306377/298‑1 (MQ=255)                                                                                                                                 gACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGttt                                                                  <  1:391225/298‑1 (MQ=255)                                                                                                                                  aaaataCAGCTAACAACGGCAGGCTTAAGATGCACGTATCAGTAGGAGAAAGACCGCATCGCCGTCACCGACAGCCTGGACCGCGGGGAAGGCGGGCGCACGCAGGGCGAAGCCGGGGTGAAGGGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGGGGGCAGGCTCAGGGCACAG‑ACGGATGCAGCAGGCAGGACAACAGAGTACAGCCCGGAGGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGtttt                                                                 <  2:709050/293‑1 (MQ=17)                                                                                                                                       aCAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGTGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCATGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACATAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACtat                                                            <  1:314391/298‑1 (MQ=255)                                                                                                                                            tAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAacca                                                       >  1:68078/1‑298 (MQ=255)                                                                                                                                             aaaCCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAaccac                                                      >  1:440265/1‑298 (MQ=255)                                                                                                                                                            aaGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACg                                       >  1:513677/1‑298 (MQ=255)                                                                                                                                                             agctACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGt                                      >  1:138663/4‑298 (MQ=255)                                                                                                                                                                tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAg                                   >  1:741225/1‑298 (MQ=255)                                                                                                                                                                   acGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCa                                <  1:604841/298‑1 (MQ=255)                                                                                                                                                                     gTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCAcc                              >  1:391091/1‑298 (MQ=255)                                                                                                                                                                      tagcagaaagagaaaGACCGCATCACCATCACCGACAGCCGGGACCGCCGTGAAGTGCGGCGCACGGAGGGAGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGGACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCg                             <  2:660064/286‑1 (MQ=17)                                                                                                                                                                            tATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCCCCACAACCAGTTAACGTCAGCCACCGGGCCTg                       >  1:791044/1‑298 (MQ=255)                                                                                                                                                                             aTGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGa                      >  1:690842/1‑298 (MQ=255)                                                                                                                                                                             aTGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGa                      >  1:191353/1‑298 (MQ=255)                                                                                                                                                                               gagaaaGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACg                    <  1:742149/295‑1 (MQ=255)                                                                                                                                                                                  aaaGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGAGGGCAGGGCATCGGCGTTTTCCTATAACCCCCACCACCCGTTAACGTCATCcccccgccccgcccggg                 >  2:57167/1‑285 (MQ=17)                                                                                                                                                                                         gCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAtt          >  1:413920/1‑298 (MQ=255)                                                                                                                                                                                         gCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGACTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAtt          >  1:559326/1‑298 (MQ=255)                                                                                                                                                                                                catcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCgg   <  1:57167/298‑1 (MQ=255)                                                                                                                                                                                                 atcaccGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCCTAACCACGCCCGAAGGCAGGGCATCGGCGCTTTACTATACCCCCCACCACCAGCTCACGTCAGCCCCCGGGCCTGGCGGGCCGGAATTGCCCCgcg  >  2:388993/1‑296 (MQ=17)                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                            GTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAAAAAGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACACAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAG‑ACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGG  >  CP006584/793778‑794267

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