Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	8,261	T→G	29.8%	L8L (CTT→CTG)	talB →	transaldolase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	8,261	0	T	G	29.8%	132.5 / 98.6	94	L8L (CTT→CTG)	talB	transaldolase B
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base G (20/8);  total (56/38)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

AAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCG  >  CP009273/8113‑8406                                                                                                                                                     |                                                                                                                                                  aaaTCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTc                                                                                                                                                  <  1:327794/150‑1 (MQ=255)   aTCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGt                                                                                                                                                <  2:609776/150‑1 (MQ=255)      cGAATCAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCa                                                                                                                                              <  2:916686/149‑1 (MQ=255)        aaaCAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGGTATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTa                                                                                                                                           >  2:270791/1‑150 (MQ=255)        aaaCAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTa                                                                                                                                           >  1:181677/1‑150 (MQ=255)        aaaCAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGAGCTGACGCTTCGTCGAGCTATCATCCGGGCAGACCCGCTATATCCCCCTAACAAGCTGTTTAAAGAGAAAAACTATCATGACGGAAAAAATGAACTCCCTTACTCAGTa                                                                                                                                           >  1:223357/1‑150 (MQ=255)            aGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGTCGTTGCGGCGTGATATCATCAGGGCAGACCGGGTACATCCCCCTAACAAGCTGGTTTAAGAGAAAAACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTacac                                                                                                                                        >  2:912578/1‑149 (MQ=255)              tctcGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTAcaccac                                                                                                                                     <  1:583714/150‑1 (MQ=255)                  gCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACaccacc                                                                                                                                    <  2:301793/147‑1 (MQ=255)                    cTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGt                                                                                                                               >  1:899347/1‑150 (MQ=255)                    cTGGCGATAACCGTCTTGTCGGCGGTTGCCCTGACCGTGCCGCGTGATATCATCAGGGTAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGCAATACTATCAAGACGGACAAATTGCCCTCCCCTCGTCAGAACACCACCGTAGt                                                                                                                               >  1:1576/1‑150 (MQ=255)                        cGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATTCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGAAAAATTGACCTCCCCTCGTCAGTAAACCACCGTAGTgg                                                                                                                             >  1:484546/1‑148 (MQ=255)                          aTAACCGTCTTGTCGGCGGTTGGGCCGACGGTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTAACATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGa                                                                                                                         >  2:807579/1‑150 (MQ=255)                           tAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGac                                                                                                                        <  2:925338/150‑1 (MQ=255)                               cGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTg                                                                                                                    <  2:727104/150‑1 (MQ=255)                                gTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTgg                                                                                                                   <  2:134287/150‑1 (MQ=255)                                        ggcggTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATcg                                                                                                           <  2:695452/150‑1 (MQ=255)                                              tGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCaa                                                                                                     >  1:101295/1‑150 (MQ=255)                                                     aCGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCt                                                                                              <  1:270791/150‑1 (MQ=255)                                                     aCGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTCTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCt                                                                                              >  2:259322/1‑150 (MQ=255)                                                      cGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTg                                                                                             <  2:560431/150‑1 (MQ=255)                                                        ttGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGATATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTa                                                                                           <  1:383395/150‑1 (MQ=255)                                                        ttGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTa                                                                                           >  1:408920/1‑150 (MQ=255)                                                              cgtGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAAcc                                                                                     >  2:436348/1‑150 (MQ=255)                                                                  atatCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAg                                                                                 >  2:297230/1‑150 (MQ=255)                                                                       atcaGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATg                                                                             <  2:194512/149‑1 (MQ=255)                                                                       atcaGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATg                                                                             <  2:194508/149‑1 (MQ=255)                                                                       atcaGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGc                                                                            >  1:860577/1‑150 (MQ=255)                                                                          aGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCca                                                                          >  2:604296/1‑149 (MQ=255)                                                                               aGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATg                                                                                                   >  1:496633/1‑119 (MQ=255)                                                                               aGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaacc                                                                     >  1:768289/1‑149 (MQ=255)                                                                                gACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTTAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCAcaaccaa                                                                   <  1:513427/150‑1 (MQ=255)                                                                                gACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGa                                                                                                  <  2:496633/119‑1 (MQ=255)                                                                                     gTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCtt                                                              <  2:157439/150‑1 (MQ=255)                                                                                      ttACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCtcc                                                             <  2:853760/150‑3 (MQ=255)                                                                                      ttACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTtc                                                             >  1:260707/1‑150 (MQ=255)                                                                                        aCATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGGCAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTtctc                                                           >  1:568320/1‑150 (MQ=255)                                                                                        aCATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTtct                                                            <  1:170824/149‑1 (MQ=255)                                                                                        aCATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTtct                                                            <  2:971781/149‑1 (MQ=255)                                                                                                cTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCt                                                    >  2:51569/1‑149 (MQ=255)                                                                                                cTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCt                                                    >  1:466837/1‑149 (MQ=255)                                                                                                 tAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCtt                                                   >  2:875442/1‑149 (MQ=255)                                                                                                   aCAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTaa                                                 <  2:1576/149‑1 (MQ=255)                                                                                                    cAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcg                                               >  2:6535/1‑150 (MQ=255)                                                                                                    cAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcg                                               <  1:451398/150‑1 (MQ=255)                                                                                                           tttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcgc                                         >  2:765271/1‑149 (MQ=255)                                                                                                           tttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcg                                          >  2:190836/1‑148 (MQ=255)                                                                                                           tttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcg                                          >  2:249534/1‑148 (MQ=255)                                                                                                           tttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcg                                          >  1:152745/1‑148 (MQ=255)                                                                                                            ttAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAAcgcagcgca                                        <  1:765271/149‑1 (MQ=255)                                                                                                                 gagaAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTcc                                  >  2:526182/1‑150 (MQ=255)                                                                                                                  agaAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCg                                 <  2:181677/150‑1 (MQ=255)                                                                                                                       tACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATa                            >  2:273270/1‑150 (MQ=255)                                                                                                                         cTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATAcc                          <  1:912578/150‑1 (MQ=255)                                                                                                                         acaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTAcc                          >  1:418829/6‑150 (MQ=255)                                                                                                                         acaaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTAcc                          >  2:134071/6‑150 (MQ=255)                                                                                                                          caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                         >  2:186198/5‑150 (MQ=255)                                                                                                                          caaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCg                         >  2:712036/5‑150 (MQ=255)                                                                                                                           aaaaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                        <  1:1054275/147‑1 (MQ=255)                                                                                                                           aTCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCg                         >  2:711548/1‑149 (MQ=255)                                                                                                                             cATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGt                        >  1:590199/1‑148 (MQ=255)                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                        >  2:707357/2‑148 (MQ=255)                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                        >  1:639322/2‑148 (MQ=255)                                                                                                                             aaTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGt                        >  1:789879/2‑148 (MQ=255)                                                                                                                              aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                     >  2:943630/1‑150 (MQ=255)                                                                                                                              aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                     >  2:836999/1‑150 (MQ=255)                                                                                                                              aTGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAAACCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAg                     >  1:206375/1‑150 (MQ=255)                                                                                                                              aTGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTaa                      <  2:101295/149‑1 (MQ=255)                                                                                                                              aTGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAg                     <  2:484546/150‑1 (MQ=255)                                                                                                                              aTGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAg                     <  1:997582/150‑1 (MQ=255)                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCGCCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                    <  1:55690/150‑2 (MQ=255)                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                    <  1:752383/150‑2 (MQ=255)                                                                                                                               tGACGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgc                    <  1:89109/150‑2 (MQ=255)                                                                                                                                 aCGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAgct                   >  1:192695/1‑147 (MQ=255)                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga                 <  2:731227/150‑1 (MQ=255)                                                                                                                                  cGGATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctga                 >  2:893542/1‑150 (MQ=255)                                                                                                                                   ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                <  2:507651/150‑1 (MQ=255)                                                                                                                                   ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                <  1:503680/150‑1 (MQ=255)                                                                                                                                   ggATAAATTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgat                <  1:311959/150‑1 (MQ=255)                                                                                                                                   ggACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgat                >  1:311964/1‑150 (MQ=255)                                                                                                                                     ataaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgatt               >  2:664426/3‑149 (MQ=255)                                                                                                                                     ataaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgatt               >  1:169481/3‑149 (MQ=255)                                                                                                                                     ataaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgatt               >  1:659944/3‑149 (MQ=255)                                                                                                                                     ataaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgatt               >  1:879118/3‑149 (MQ=255)                                                                                                                                     aCAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGttgattg              <  1:807579/150‑1 (MQ=255)                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGTTATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga             >  2:508066/2‑150 (MQ=255)                                                                                                                                      taaaTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGctgattga             >  1:458598/2‑150 (MQ=255)                                                                                                                                       aaaTTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTtgattgat            <  1:1042211/150‑1 (MQ=255)                                                                                                                                       aaaTTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTtgattgat            >  1:929464/1‑150 (MQ=255)                                                                                                                                       aaaTTGACCTCCCTTCGTCAGTACACCACC‑‑‑GTGGCCGCCACTAGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGCATGCCACAACCACCCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGATGATTgttgat         >  2:259358/1‑150 (MQ=255)                                                                                                                                        aaTTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTtgattgat            >  2:985862/1‑149 (MQ=255)                                                                                                                                         aTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATtgacga          >  2:228577/1‑150 (MQ=255)                                                                                                                                         aTTAACTTCTCTGCGTCAGTACACCACTGTCGTCGCTGACACCGGAGATATTGCGGCAATGAAGCTGTATCAGCCTCAGGATGCCACCACTAACCCTTCCCTGATTCTCGGCGCGGCTCAGATCCCTGAGTACCGTAAGCTGATtgacga          >  1:876880/1‑150 (MQ=255)                                                                                                                                                 cccTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCg  >  1:475640/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AAATCCGAAACAGTCTCGCCTGGCGATAACCGTCTTGTCGGCGGTTGCGCTGACGTTGCGTCGTGATATCATCAGGGCAGACCGGTTACATCCCCCTAACAAGCTGTTTAAAGAGAAATACTATCATGACGGACAAATTGACCTCCCTTCGTCAGTACACCACCGTAGTGGCCGACACTGGGGACATCGCGGCAATGAAGCTGTATCAACCGCAGGATGCCACAACCAACCCTTCTCTCATTCTTAACGCAGCGCAGATTCCGGAATACCGTAAGTTGATTGATGATGCTGTCG  >  CP009273/8113‑8406

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