Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	360,147	G→C	29.5%	intergenic (‑81/+514)	lacZ ← / ← lacI	pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	360,147	0	G	C	29.5%	191.8 / 32.4	86	intergenic (‑81/+514)	lacZ/lacI	pseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (33/27);  new base C (22/4);  total (55/31)
Fisher's exact test for biased strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAA  >  CP009273/359998‑360295                                                                                                                                                      |                                                                                                                                                        ggcggatttgtcctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                         >  2:98792‑M1/107‑150 (MQ=255)  aCATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCgg                                                                                                                                                        <  1:724115/150‑1 (MQ=255)       cAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                                                                                                                                   <  1:101323/150‑1 (MQ=255)       cAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGcccg                                                                                                                                                   >  2:912662/1‑148 (MQ=255)            gCACTTCACCGCTTGCCAGCGGCTTTCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                              <  1:906831/150‑1 (MQ=255)            gCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCCGCCCTTGcct                                                                                                                                              >  2:423787/1‑148 (MQ=255)            gCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGCCGCTGCtt                                                                                                                                              >  2:665365/1‑150 (MQ=255)               cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTGGCCCCCCCCCCCCCGGCCCGTGCTTccc                                                                                                                                           >  2:505395/1‑148 (MQ=255)                 tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGGCCGCCACCCCCCCGGCCGCTGCTTCGC‑aa                                                                                                                                        >  1:676633/1‑150 (MQ=255)                    ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGCTGCTTCGC‑AAC‑G‑t                                                                                                                                   >  1:126647/1‑150 (MQ=38)                     cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑tt                                                                                                                                  <  1:581883/150‑1 (MQ=255)                        ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCCCTGCTTCGC‑CAC‑G‑GTCa                                                                                                                                >  2:498679/1‑149 (MQ=38)                            ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATcc                                                                                                                           >  2:379871‑M1/80‑150 (MQ=255)                            ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATcc                                                                                                                           >  1:134017‑M1/80‑150 (MQ=255)                            ttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATcc                                                                                                                           >  1:948069‑M1/80‑150 (MQ=255)                             tcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGGGTCTTGGCCCCCCCCCCCCCGGCCGTTGCTTCCCCAAC‑g‑ttcaatacc                                                                                                                           <  2:950180‑M1/72‑5 (MQ=255)                                ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTc                                                                                                                       <  1:980988/150‑1 (MQ=255)                                  cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCTC‑CAC‑G‑TTCAAAACCGCcccc                                                                                                                     >  2:5488/1‑150 (MQ=255)                                   acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCcg                                                                                                                    >  1:264962‑M1/73‑150 (MQ=255)                                    caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCcgg                                                                                                                   >  1:657522‑M1/72‑150 (MQ=255)                                     aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCCC‑AAC‑GCTTCAAACCCCCCCCcgc                                                                                                                   >  1:3175/1‑149 (MQ=37)                                      ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCcggc                                                                                                                  <  1:112397/149‑1 (MQ=255)                                      ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGCCGGTGCTTCGC‑AAC‑G‑CTCAAAACCGCTCCccgcg                                                                                                                 >  1:968804/1‑150 (MQ=255)                                       cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGGCGCTGCTTCCC‑CAC‑G‑GTTCAATCCCCTccccccgg                                                                                                                >  2:264962/1‑145 (MQ=255)                                       acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCcggcgg                                                                                                                >  2:76578‑M1/69‑150 (MQ=255)                                       acaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCcggcgg                                                                                                                >  2:400182‑M1/69‑150 (MQ=255)                                          ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGCCCCCTCCCCCCCCCCCCCGGCCGTTGTTTCGC‑AAC‑G‑GTCAAAACCCCCCCCCCCGGact                                                                                                             >  1:837883/1‑148 (MQ=255)                                           cagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGAttt                                                                                                            >  1:950180‑M1/65‑150 (MQ=255)                                           cagagaaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGAGGACTGCCCGCCACCCTCAGGGAAGTTGATTAGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGAttt                                                                                                            >  1:3966‑M1/65‑150 (MQ=255)                                           cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGAtt                                                                                                             <  1:776989/149‑1 (MQ=255)                                            aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTg                                                                                                           <  1:835517/150‑1 (MQ=255)                                                aaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGCTGCTTCCC‑AAC‑A‑TTTAAATCCCCTCCCCGCGGGTTTGTCCt                                                                                                       <  1:614721‑M1/91‑1 (MQ=255)                                                 aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                       >  1:458259‑M1/59‑149 (MQ=255)                                                 aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                      >  1:274018‑M1/59‑150 (MQ=255)                                                  aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGGCGGTGCTTCGC‑AAA‑G‑TTCAAAACCGCCCCcggcgggttttttcttc                                                                                                     <  1:644563‑M1/93‑8 (MQ=255)                                                  aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                     >  2:482312‑M1/58‑150 (MQ=255)                                                   acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                    >  2:687100‑M1/57‑150 (MQ=255)                                                    caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                   <  2:996790/150‑1 (MQ=255)                                                     gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                  >  1:452142‑M1/55‑150 (MQ=255)                                                      tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCAACTCAg                                                                                                 <  2:978393/150‑1 (MQ=255)                                                       aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAgg                                                                                                >  2:499014‑M1/53‑150 (MQ=255)                                                          gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGCGGGCCTCCCCCCCCCCCTTCCGGCCCTTGGTTTTC‑AAC‑C‑ACTAAATACGCACCCGGCGGGTTAGTGCTCCTCTggg                                                                                               >  1:23570/1‑147 (MQ=255)                                                           tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                            <  2:158313/150‑1 (MQ=255)                                                           tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                            <  1:936530/150‑1 (MQ=255)                                                           cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                            >  1:246373‑M1/49‑150 (MQ=255)                                                              agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTCGC‑AAC‑C‑TTCAAATCCGCTCCCCGCGGGTTTTTTCTACTCAGGAGAGcgg                                                                                         <  1:318099‑M1/105‑2 (MQ=255)                                                              aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCCGGCGTTGGTTTGC‑GACGG‑TTCAAAACC‑CTCCCCGCGGGTTTTTTCTTCCCCGGAGAGCgg                                                                                         >  2:658904/1‑149 (MQ=255)                                                                tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                       >  2:487660‑M1/44‑150 (MQ=255)                                                                tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                       >  2:126009‑M1/44‑150 (MQ=255)                                                                 ggCTCCGATAATAAGCTGTCAAACATGAGATTTCCGAAGCCCAATACCGGGCAGTTCATGCCGGGCCTCCTCCCCGCCGCCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGTTCa                                                                                      <  1:324663/149‑1 (MQ=17)                                                                  gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                     <  1:863818/150‑1 (MQ=255)                                                                  gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                     <  1:505395/150‑1 (MQ=255)                                                                   cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                    <  1:385816/150‑1 (MQ=255)                                                                     cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                   >  1:62996‑M1/39‑149 (MQ=255)                                                                      cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                  <  2:968760/149‑1 (MQ=255)                                                                       actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                >  1:344418‑M1/37‑150 (MQ=255)                                                                            tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                           <  1:377948/150‑1 (MQ=255)                                                                                   tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGCTTTGC‑AAA‑G‑TTCAAATCCCCCCCCCGCGGGGTTGTCCTACTCAGGAGAGCGGTCACCCccaaaaaaaaaaaa                                                                     >  1:98792/1‑143 (MQ=17)                                                                                   gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                    >  2:610020‑M1/25‑150 (MQ=255)                                                                                       aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                <  1:924767/150‑1 (MQ=255)                                                                                            gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                            >  2:968804‑M1/16‑149 (MQ=255)                                                                                             atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                           >  1:1010426‑M1/15‑149 (MQ=255)                                                                                              aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                         <  1:154389/150‑1 (MQ=255)                                                                                                cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                        >  2:325786‑M1/12‑149 (MQ=255)                                                                                                   cGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTGCTTCGC‑AAC‑G‑TTTAAATCCCCTCCCCGCGGGTTTGTCCTTCTCAGGGGAGCGGTCACCCACAAAAAACAAAAAAAACGAAAAGCCCAggc                                                    >  1:400182/1‑148 (MQ=14)                                                                                                    gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                   >  2:186426‑M1/8‑150 (MQ=255)                                                                                                    gaaaattAGTTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                   >  2:114582‑M1/8‑150 (MQ=255)                                                                                                    gATCCCCAATTCCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                   <  2:837883/150‑1 (MQ=38)                                                                                                    gATCCCCAATTCCTGGCAGTTTATGGCGGGGGGCCTGCCCGCCCCCCCCCCGGCCGTTGCTTCGC‑CAC‑G‑GTCAAAAACGCCCCCCGCGGGTTTGTTCCACTCCGGAGAGCGGTCACCCACAAAAAACAGAAAAAAAGAAAAAGCCAGTCt                                                   >  2:368620/1‑150 (MQ=255)                                                                                                     aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGC‑AAA‑G‑GTCAAATCCCCCCCCCGCGGGTTTGTCCTACTCAGGAGAGCGGTCACCGCCAAACAACAAATAAAAAGAAAAGccccatcct                                                  >  1:76578/1‑144 (MQ=11)                                                                                                      gaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                 >  2:596963‑M1/6‑150 (MQ=255)                                                                                                       ccccAATTCCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCCTTGCTTTGC‑AAC‑G‑CTCCAATCCCCCCCCCGCGGGTTTGTTCTTCCCCGGAGGGGGGTCACCGACAAACAAAAGATAAAAAAAAAAGCCCACTCTTTc                                                >  2:246373/1‑150 (MQ=2)                                                                                                       aattaaTTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCTC‑CAA‑G‑TTCAAAATCCCCCCCCGCGGGTTTTTTCTACTCCGGGGAGCGGTCACCGCCAAACAACAAAAAAAAAAAAAAGACCAGACTTTc                                                >  2:164473/5‑150 (MQ=2)                                                                                                       aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                >  1:1011971‑M1/5‑150 (MQ=255)                                                                                                        cccAATTCCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGC‑AAC‑G‑TTCAAAACCCCTCCCCGCGGGGTTTTTCTACTCAGGGGAGCGTTCACCCCCAAAAAACAAATAAAAAAAAAAGCCCAGGCTTTc                                                >  1:308148/1‑149 (MQ=2)                                                                                                          taaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTTCTTCGC‑AAC‑G‑TTCAAAACCCCCCCCCGCGGGTTTTTCCTTCTCCAGAGAGCGTTCCCCCACAAACAACAGAAAAAAAGAAAAGCCCACGCTTTCCAc                                             >  1:968443/2‑150 (MQ=2)                                                                                                           aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                             <  1:423787/149‑1 (MQ=32)                                                                                                             ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                          <  1:5488/150‑1 (MQ=31)                                                                                                                                         cccGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                <  1:841205/148‑1 (MQ=11)                                                                                                                                          ccGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa             <  1:912662/150‑1 (MQ=32)                                                                                                                                           cGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat            >  2:3966/1‑150 (MQ=33)                                                                                                                                           cGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACAGAGCCTTTCGTTTTATTTGATGCCTGgggaat            >  1:418180/1‑150 (MQ=18)                                                                                                                                             ccACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt           >  2:274018/1‑149 (MQ=34)                                                                                                                                              cACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgg         >  1:287132/1‑150 (MQ=35)                                                                                                                                                  cTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat     <  2:340298/150‑1 (MQ=39)                                                                                                                                                     cGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaa  >  1:530933/1‑150 (MQ=255)                                                                                                                                                     cGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaa  <  1:405069/150‑1 (MQ=255)                                                                                                                                                     cGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaatta   <  1:908124/149‑1 (MQ=40)                                                                                                                                                      |                                                                                                                                                        GACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGC‑AAC‑G‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAA  >  CP009273/359998‑360295

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