Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,241,203:1	+G	6.1%	coding (72/663 nt)	yqjA →	general envelope maintenance protein; DedA family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,241,202	1	.	G	6.1%	278.4 / 12.5	66	coding (71/663 nt)	yqjA	general envelope maintenance protein; DedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (32/30);  new base G (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

TACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTC  >  CP009273/3241054‑3241348                                                                                                                                                      |                                                                                                                                                   tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat                                                                                                                                                     <  1:28226/149‑1 (MQ=255)     aaaaCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GAtt                                                                                                                                                <  1:219643/149‑1 (MQ=255)           gCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAAT‑CATCGATGGATTGGccag                                                                                                                                          <  1:1104103/149‑4 (MQ=255)            cAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGCGTTTGAAACCCTCGCCAAGCCATCGAT‑GATTg                                                                                                                                               >  2:126128/1‑143 (MQ=255)              aTAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCAtgttg                                                                                                                                       >  1:415980/1‑149 (MQ=255)               tAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAACCCATCGCT‑GGCTGGCa                                                                                                                                            >  2:266549/1‑143 (MQ=255)               tAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCAAATTAATGCAATTACTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGcg                                                                                                                                            >  2:680098/1‑142 (MQ=255)               tAACACTGCGTAAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGCTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGAGTTTGAAACCCTGGCCAAT‑CATCGATGGATTGGccagggtt                                                                                                                                      <  1:262577/149‑8 (MQ=255)                aaCATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTa                                                                                                                                     <  1:817703/149‑1 (MQ=255)                               cATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCtt                                                                                                                            >  2:937338/1‑143 (MQ=255)                                aTCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCACAGGGTTTTGAAACCCTGGCCAATCCCTCGAT‑GATTGGCAGGTTGTATTTTGTCttg                                                                                                                           >  2:892725/1‑143 (MQ=255)                                aTCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTGGACCCAATTGCTGCAAGCCCTGTGTGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GAGTGGCATGTTGGATTTTGTCTTGTTTGTa                                                                                                                     >  1:276320/1‑149 (MQ=255)                                           ccAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAAtttt                                                                                                                <  2:791684/143‑1 (MQ=255)                                               cGTAGCGTTAACCGATTCACCAGGAATAATGAATGTAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaa                                                                                                      <  1:124939/149‑1 (MQ=255)                                                  ggCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCtt                                                                                                         <  2:948675/143‑1 (MQ=255)                                                    cGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGa                                                                                                       >  2:876163/1‑143 (MQ=255)                                                    cGTTAACCAATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGa                                                                                                       <  2:45955/143‑1 (MQ=255)                                                        aaCCTATCCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCACTCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttg                                                                                             <  1:1240454/149‑1 (MQ=255)                                                        aaCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAAc                                                                                                   >  2:254076/1‑143 (MQ=255)                                                             aTTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGCCTTGTTTGTAATTTTGTTCCTTGAAAACGgctt                                                                                              >  2:441684/1‑143 (MQ=255)                                                               tCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTCTTCCTTGAAAACGGCTTGCTTCcgg                                                                                      <  1:652752/149‑1 (MQ=255)                                                                 aCCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcg                                                                                    >  1:543298/1‑149 (MQ=255)                                                                   cAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc                                                                                  <  1:299306/149‑1 (MQ=255)                                                                        ataatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAGAACGGCTTGCTTCCGGCGGCCtttt                                                                             <  1:775516/149‑1 (MQ=255)                                                                        ataatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCtttt                                                                             >  1:1243831/1‑149 (MQ=255)                                                                        ataatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGATTTcggcgt                                                                                   >  2:704511/1‑142 (MQ=255)                                                                         taatGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTCGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTTTTTGTAATTTTGTTCCTTCAAAACGGCTTGCTTCCGGCGGCCttttt                                                                            <  1:732018/149‑1 (MQ=255)                                                                           atgactgGAACTTTTGACCCAATCGGTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTTATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                                                                          <  1:256748/149‑1 (MQ=255)                                                                           atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAATGGCTTGCTTCCGGCGGCCTTTTTAc                                                                          <  1:727214/149‑1 (MQ=255)                                                                           atgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                                                                          >  1:1038004/1‑149 (MQ=255)                                                                               atgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCAGCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCggg                                                                      <  1:3005/149‑1 (MQ=255)                                                                                tgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGc                                                                     <  1:252103/149‑1 (MQ=255)                                                                                       tttGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTCCCGGGCg                                                                    >  2:523705/1‑143 (MQ=255)                                                                                             ccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg                                                        >  1:870137/1‑149 (MQ=255)                                                                                                aTTGCTGCAAGCCCTATGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTAc                                                           >  2:623461/1‑143 (MQ=255)                                                                                                  tgctgcCCGCCCTGTGGGCGCAGTATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTg                                                   <  1:1266797/149‑1 (MQ=255)                                                                                                   gctgcAAGCCCTGTGGGCGCATGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATTTTGTATTTTGTCTTGTTTGTAATTTTGTTCATTGTATTCTGCGGGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTgg                                                  <  1:255771/149‑1 (MQ=255)                                                                                                    ctgcACGTCCTGTGGGCGCAGGATTTTGAAGCCCTGCCCAAGCCATCGAC‑TATTGGCAGGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGt                                                 <  1:93830/149‑1 (MQ=255)                                                                                                    ctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGTCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGt                                                 <  1:842151/149‑1 (MQ=255)                                                                                                       cAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGCCGACAGTTTACTGGTCtt                                                    >  2:569114/1‑143 (MQ=255)                                                                                                        aaGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTg                                                   <  2:998518/143‑1 (MQ=255)                                                                                                         aGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCg                                            <  1:1206355/149‑1 (MQ=255)                                                                                                           cgcTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTc                                                >  2:300588/3‑143 (MQ=255)                                                                                                           cccTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGCCAGTTTACTGGTATTGGTCGGCgtg                                          >  1:277137/1‑149 (MQ=255)                                                                                                                  ggCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTg                                   <  1:999317/149‑1 (MQ=255)                                                                                                                    agccgcaTTTTGAAACCCTGGTCAATCCATCGGT‑GATTGGCATGTTGTTTTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGGTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTg                                       <  2:1085899/137‑1 (MQ=255)                                                                                                                     gcAGGATTTTGAAACCCTGGCCAATGCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGa                                <  1:1073559/149‑1 (MQ=255)                                                                                                                       aGGATTTTGAAACCCTGGCCAATCCATTGAT‑GAGTGGACTGTTGTATTTTGTTTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCCGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGaaa                              <  1:428134/149‑1 (MQ=255)                                                                                                                         gATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAgg                            >  1:1024270/1‑149 (MQ=255)                                                                                                                         gATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAgg                            >  1:107930/1‑149 (MQ=255)                                                                                                                           ttttGAACCCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGa                                <  2:638651/143‑1 (MQ=255)                                                                                                                           ttttGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGa                                >  2:29098/1‑143 (MQ=255)                                                                                                                              tGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGACAGGCGCGa                       >  1:1214831/1‑149 (MQ=255)                                                                                                                              tGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGa                       <  1:1162319/149‑1 (MQ=255)                                                                                                                                 aaCCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAACACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATgg                    >  1:789076/1‑149 (MQ=255)                                                                                                                                  caacatgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATggg                   >  1:168682/7‑149 (MQ=255)                                                                                                                                  caacatgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATggg                   >  1:979955/7‑149 (MQ=255)                                                                                                                                     cTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTa                <  1:674898/149‑1 (MQ=255)                                                                                                                                       ggCCAATCCATCGAT‑GATTGGCAGGTTGTATTTTGTCTTGGTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGCTTGCGAAAGGCGCGATgg                    >  2:1100406/1‑143 (MQ=255)                                                                                                                                        gtcaaaacatgcat‑gatTGGCATGTTGCCTTTGGTCTTGTTTCTATTTTTTTTCCGTGAGAACGGCTGGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATggg                   <  2:223519/131‑1 (MQ=255)                                                                                                                                        gCCAATCCATCGGT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATggg                   <  2:41342/143‑1 (MQ=255)                                                                                                                                          cAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTTAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGc           >  1:575177/1‑149 (MQ=255)                                                                                                                                           aaTCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACCGTTTATTGGTTTGGGTCGGCGTGTTGATTGCTGAAGGCGCGATGGGCTATCCGCa          >  1:298229/1‑149 (MQ=255)                                                                                                                                           aaTCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCa          >  1:258387/1‑149 (MQ=255)                                                                                                                                             tCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCaaa        <  1:550914/149‑1 (MQ=255)                                                                                                                                              ccATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATcc             >  2:293229/1‑143 (MQ=255)                                                                                                                                                 tCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCaaaaggt    >  1:614044/1‑145 (MQ=255)                                                                                                                                                 tCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGAt    >  1:606799/1‑149 (MQ=255)                                                                                                                                                   gat‑gatTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAACCGATTc  >  1:1103349/1‑149 (MQ=255)                                                                                                                                                   gat‑gatTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAAAGATTc  >  1:1097359/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTC  >  CP009273/3241054‑3241348

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