Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,529,302	G→A	6.0%	R272C (CGT→TGT)	yddE ←	PhzC‑PhzF family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,529,302	0	G	A	6.0%	192.7 / 2.2	67	R272C (CGT→TGT)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base G (20/42);  new base A (2/2);  total (23/44)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGCATTGCCAGTCACCG  >  CP009273/1529154‑1529448                                                                                                                                                     |                                                                                                                                                   aaaCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACg                                                                                                                                                    >  2:419393/1‑149 (MQ=255)   aCGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAACTGGTGCATTTGTCCGGTTGGTTATAACGaa                                                                                                                                                  >  1:835710/1‑149 (MQ=255)    cGCCTTATCCGGCATACATTCGCCCGGTTGTAGACCTGAGAGGCGAAGGACACCGGCCAACGGTTTCAGAGCTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAAt                                                                                                                                                 >  1:294166/1‑149 (MQ=255)    cGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATGGCCCATGCGGCATGAAACAAAGTCACCGCTGTGCCTAGAATGGTGACTTTGTCCGGTTGGTTATCACTAAt                                                                                                                                                 >  1:258448/1‑149 (MQ=255)     gCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCCCCGCTGTGGCTAAAATGGTGAATTTTTCCGGTTGGTGATCGCGAAtt                                                                                                                                                >  1:31266/1‑149 (MQ=255)      ccTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTg                                                                                                                                               >  2:1201693/1‑149 (MQ=255)       cTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGt                                                                                                                                              <  1:330654/149‑1 (MQ=255)         tATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATTGTGACTTTTTCCGGTTGGTTATCACGAATTGTCa                                                                                                                                            <  2:1191019/149‑1 (MQ=255)           tCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCAcc                                                                                                                                          <  1:132196/149‑1 (MQ=255)             cGGCATACATTAGCCCGGTTGTCGTCCTGACCGGCGAAGCGCATCAGACAACGGTTTCAGCGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTGTCACGAATTGTCACCGt                                                                                                                                        <  2:1017302/149‑1 (MQ=255)              ggCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTc                                                                                                                                       <  2:1020632/149‑1 (MQ=255)               gCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCa                                                                                                                                      <  1:102083/149‑1 (MQ=255)                 aTACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATGACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACAAATTGTCACCGTCACt                                                                                                                                    >  1:800984/1‑149 (MQ=255)                 aTACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACAAATTGTCACCGTCACt                                                                                                                                    >  1:1022481/1‑149 (MQ=255)                  tACATTAGCCCGGTTGTAGACCTGACAGTCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACtt                                                                                                                                   <  1:946033/149‑1 (MQ=255)                  tACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACtt                                                                                                                                   <  1:996320/149‑1 (MQ=255)                  tACATTAGCCCGGTTGGATACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACtt                                                                                                                                   <  2:587023/149‑1 (MQ=255)                            cGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATAcc                                                                                                                         <  1:1107999/149‑1 (MQ=255)                                tGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCa                                                                                                                     <  1:1045247/149‑1 (MQ=255)                                tGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCa                                                                                                                     <  1:440298/149‑1 (MQ=255)                                         gACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAAt                                                                                                            >  2:613381/1‑149 (MQ=255)                                           cAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATgc                                                                                                          <  1:813525/149‑1 (MQ=255)                                           cAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATgc                                                                                                          <  1:735642/149‑1 (MQ=255)                                                          aGACAACGGTTTCTGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCGGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGcc                                                                                           <  2:1018089/149‑1 (MQ=255)                                                          aGACAACGGTTTCAGAGTTCAATTGGCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGcc                                                                                           <  2:219547/149‑1 (MQ=255)                                                          aGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTTGTGATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGcc                                                                                           >  2:444234/1‑149 (MQ=255)                                                          aGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGcc                                                                                           <  2:956528/149‑1 (MQ=255)                                                            aCAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGGCtc                                                                                         <  2:1198858/149‑2 (MQ=255)                                                            aCAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCtt                                                                                         <  2:252822/149‑1 (MQ=255)                                                               aCGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTaa                                                                                      <  1:1069426/149‑1 (MQ=255)                                                                 ggTTTCAGGGGTCAGTTGCCCATTCGGCATGAAACAAACTCACCGCTGGGCCTGAACTGGTGCCTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAaca                                                                                    <  2:814454/149‑1 (MQ=255)                                                                 ggTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTTACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAaca                                                                                    <  1:182338/149‑1 (MQ=255)                                                                 ggTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTCCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAaca                                                                                    <  2:814281/149‑1 (MQ=255)                                                                   tttCAGAGTTCAATTGTCCATTCGGTATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTCATCACGAATTGTCACCGTCACTTCAATCATACCGCCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACg                                                                                  <  1:656479/149‑1 (MQ=255)                                                                              aaTTGGCGAGTCGGCATGGAACAAAAGCGCCGCTGGGCCTGAAATGGGGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCGATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGc                                                                       <  2:815367/149‑1 (MQ=255)                                                                               aTTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTAcc                                                                      >  1:785343/1‑149 (MQ=255)                                                                                 tGCCCATTCGGCATGAAACAAAATCACCCCAGTGACAAAAATTTGCGCTTTGGCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCATGATGGCCTTTAACACGCAACACGTTGCCAt                                                                    <  2:258448/149‑1 (MQ=255)                                                                                   cccATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTTCCATCg                                                                  >  1:1188083/1‑149 (MQ=255)                                                                                       ttCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTggg                                                              <  2:423398/149‑1 (MQ=255)                                                                                       ttCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTggg                                                              <  1:522832/149‑1 (MQ=255)                                                                                       ttCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTAACCGTCACTTCAATCATCCCGTCCCGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAAAACGTTGCCATCATggg                                                              >  2:36458/1‑149 (MQ=255)                                                                                        tCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACAAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTgggg                                                             <  2:800984/149‑1 (MQ=255)                                                                                        tCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACAAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTgggg                                                             <  2:1022481/149‑1 (MQ=255)                                                                                          ggCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCaaacaggggca                                                           >  1:974428/1‑139 (MQ=255)                                                                                          ggCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCa                                                           >  1:1118797/1‑149 (MQ=255)                                                                                           gCATGAAACAAAATCACCAGTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCac                                                          >  1:911291/1‑148 (MQ=255)                                                                                            cATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGATAAGTTGTTTTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAAt                                                         >  1:658948/1‑149 (MQ=255)                                                                                              tGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATAc                                                       >  1:313972/1‑149 (MQ=255)                                                                                              tGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATAc                                                       >  1:994349/1‑149 (MQ=255)                                                                                                 aaCAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGtt                                                    >  1:773865/1‑149 (MQ=255)                                                                                                 aaCAAAACCCCCGCTGTGCCTGAACTGGTTACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCATTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGtt                                                    <  2:413496/149‑1 (MQ=255)                                                                                                      caTCACCGCCGTGCCTGAAATGGTGACTTTTACCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGAt                                               <  1:421242/148‑1 (MQ=255)                                                                                                      aaTCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGAt                                               <  1:1191156/149‑1 (MQ=255)                                                                                                      aaTCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGAt                                               <  1:1201693/149‑1 (MQ=255)                                                                                                       aTCCCCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCATCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATg                                              <  2:407722/149‑1 (MQ=255)                                                                                                                  gCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGACACTTCATTCATACCGTCACGCCCCGATGCGCGCCCCTGATGGCCTTTAACGCGCAACACTTTGCCATCGTGGTGCTATACGTTGTGATGTACCAACCATg                                   <  1:914855/149‑1 (MQ=255)                                                                                                                  gCCTGAAATGGTGACTTTTGCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATg                                   <  1:1131558/149‑1 (MQ=255)                                                                                                                    cTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCa                                 >  1:1195483/1‑149 (MQ=255)                                                                                                                           ggctacctgttCCGGGTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCg                          <  1:486963/140‑1 (MQ=255)                                                                                                                                 tttttCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCAt                    <  1:691665/149‑1 (MQ=255)                                                                                                                                     tCCGGTTCGTTAGCAAGAATTGTCACCGTCACTTCAATCGTACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGAGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGc                <  2:31266/149‑1 (MQ=255)                                                                                                                                       cGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGGCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCAGCGGGCCATTGGCAt              <  2:456493/149‑1 (MQ=255)                                                                                                                                         gttggttATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGCATTg            <  1:529060/149‑1 (MQ=255)                                                                                                                                            ggttATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGCATTGCCa         <  1:303308/149‑1 (MQ=255)                                                                                                                                              ttATCACGACTTGTCACCGTCACTTCAATCATACCGCCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGCATTGCCAGt       <  1:419393/149‑1 (MQ=255)                                                                                                                                                 tCACGAACTGTCACCGTCATTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCCCCACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCGACCATGCCCCCATCGGGCCATTGGCATTGccctttac    >  2:146920/1‑143 (MQ=255)                                                                                                                                                   aCGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCCCCCATCGGGCCATTGGCATTGCCAGTCACCg  >  1:603103/1‑149 (MQ=255)                                                                                                                                                   aCGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGCATTGCCAGTCACCg  >  1:470666/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AAACGCCTTATCCGGCATACATTAGCCCGGTTGTAGACCTGACAGGCGAAGCGCATCAGACAACGGTTTCAGAGTTCAATTGCCCATTCGGCATGAAACAAAATCACCGCTGTGCCTGAAATGGTGACTTTTTCCGGTTGGTTATCACGAATTGTCACCGTCACTTCAATCATACCGTCACGCCCCAATGCGCGCCCCTGATGGCCTTTAACACGCAACACGTTGCCATCGTGGGGCAATACGTTGTGATGTACCAACCATGCACCCATCGGGCCATTGGCATTGCCAGTCACCG  >  CP009273/1529154‑1529448

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