Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	24.5%	S27S (AGT→AGC)	tpn ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,471	0	A	G	24.5%	102.7 / 35.5	68	S27S (AGT→AGC)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (23/28);  new base G (6/11);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCG  >  CP000730/36331‑36608                                                                                                                                             |                                                                                                                                          cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGGATCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                           <  1:308992/141‑1 (MQ=255) cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAATCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                           <  1:308990/141‑1 (MQ=21)  gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                                          <  1:148196/141‑1 (MQ=21)   ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                         >  2:315957/1‑141 (MQ=255)   ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAAcc                                                                                                                                         >  1:201574/1‑140 (MQ=19)   ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTGAATATTTCAGATATATCACGATAAcc                                                                                                                                         >  2:112488/1‑140 (MQ=19)    cATTTGTAATAAGCTTTTTTATGTTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                        <  1:102302/141‑1 (MQ=21)    cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATCACTc                                                                                                                                        <  2:247852/141‑1 (MQ=255)    cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGGTAGCta                                                                                                                                        >  2:30517/1‑140 (MQ=11)    cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                        <  2:304128/141‑1 (MQ=11)    cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                                        <  2:161364/141‑1 (MQ=11)     aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTTTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAACtaa                                                                                                                                       >  1:112061/1‑139 (MQ=255)     aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                       <  2:81945/141‑1 (MQ=255)      tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  2:210928/141‑1 (MQ=11)      tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  2:295846/141‑1 (MQ=11)      tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  1:146054/141‑1 (MQ=39)      tttGTAATAAGCTTTTTTATGCTTTTTCGTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                      <  2:341641/141‑1 (MQ=21)      tttGTAATAAGCTCTTTTATGCTTCTTCTTCCAAATTTGATATAAAATTCGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                      <  1:70822/141‑1 (MQ=11)       ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                     >  2:333132/1‑141 (MQ=11)       ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                     >  1:18726/1‑141 (MQ=11)       ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGTTATATCACGATAACTCAAt                                                                                                                                     <  1:197048/141‑1 (MQ=21)        tggaataaGCTTTGTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTGTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                    <  2:59215/138‑1 (MQ=11)        tGTAATAAGCTTTTTTATGCTTTTTCTTTCAAAGTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGGAGAGCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCGCAATg                                                                                                                                    <  2:70160/141‑1 (MQ=11)        tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGATCGTTTACACCACGTTCCCTTAATTTTTCAGATATATCACGATAACTCAATg                                                                                                                                    <  2:253898/141‑1 (MQ=255)        tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                    >  1:153226/1‑141 (MQ=21)         ttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                                                   <  2:227894/140‑2 (MQ=11)         gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                                                   <  1:119745/141‑2 (MQ=18)          taataaGCTTTTTTTTTTTTTTTTTTCCAAAATTGAGATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGCCTGCCGAACTCTTCCACCACGTTTTCTTATTTTTTCTGATTTTTCTCGATAGCTCAATgca                                                                                                                                  >  2:38327/1‑139 (MQ=2)          taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                  <  1:71094/141‑3 (MQ=11)          taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                  <  1:94592/141‑3 (MQ=17)          taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                  <  2:18623/141‑1 (MQ=40)           aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                 >  1:14379/1‑138 (MQ=11)           aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgcat                                                                                                                                 >  2:301604/1‑138 (MQ=11)           aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtat                                                                                                                                 >  1:198179/1‑141 (MQ=255)                 cTTTTGTGTTCTTTTTCTTACAAATTTGATACAAAATTGGGGTATATTCTTGAACCCAACGGTAGAACGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTa                                                                                                                           <  2:157256/141‑1 (MQ=11)                    ttgtATGCTTGTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAt                                                                                                                        <  1:165218/138‑1 (MQ=25)                       tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATCTTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:161027/1‑141 (MQ=25)                       tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                     >  1:238464/1‑141 (MQ=35)                         tGCTTTTTCTTCCAAATTTGATACAAAATTGGGGTATATTCTTGAACCCAACGGTAGACCGTTAAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAGATATCTCACGATAACTCCATGTATTTCTTACAtagtag                                                                                                                   >  2:63843/1‑141 (MQ=11)                          gCTTTTTCTTCCAAATTTTATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                  <  1:13687/141‑1 (MQ=21)                                 cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                           <  2:63515/141‑1 (MQ=35)                                   tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                         <  1:102523/141‑1 (MQ=35)                                   aacAAAGTTTATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTCTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                         <  2:266484/139‑1 (MQ=255)                                          ttGTTACAAAATTTGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGCTGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                  <  1:160250/141‑1 (MQ=255)                                          ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                  <  2:7160/141‑1 (MQ=35)                                           tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCAAACGGCTACAGTGATa                                                                                                 >  2:267110/1‑141 (MQ=25)                                                      tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                      <  2:105824/141‑1 (MQ=32)                                                      tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                      <  2:65203/141‑1 (MQ=32)                                                      tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                      >  2:285829/1‑141 (MQ=32)                                                       ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                     >  1:320257/1‑141 (MQ=32)                                                       ggggCATATTCTTGAACCCAACGGTAGACCGGTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAAATCCTTGTTa                                                                                     >  1:183682/1‑141 (MQ=17)                                                                tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                            <  1:359210/141‑1 (MQ=32)                                                                tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                            <  1:256047/141‑1 (MQ=32)                                                                     aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                       <  2:156032/141‑1 (MQ=32)                                                                                aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                            >  1:205800/1‑141 (MQ=32)                                                                                    cGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACa                                                        <  1:208113/141‑1 (MQ=32)                                                                                           tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAcc                                                 >  1:341572/1‑140 (MQ=17)                                                                                                   gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGACACAACAGTTTACAGCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                         >  1:170544/1‑141 (MQ=11)                                                                                                     ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACCTCCTTGTTAAATTGTTTCTATCTGAAATAGTTCATACAGAAGACTCCTTTTAGtt                                       >  2:207319/1‑141 (MQ=12)                                                                                                             cGTTTCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                               <  2:302329/141‑1 (MQ=17)                                                                                                             cGTTCCCTGAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                               <  2:302431/141‑1 (MQ=17)                                                                                                                   cTTAATATTCCAGATATATCACGATAACTCAATGTATATCTTAGATAGCAGCCAACGGCTACAGTGATAACATCCTTGTCAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtata                         <  1:261740/141‑1 (MQ=12)                                                                                                                        tttgtCAGATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTc                    <  2:256962/137‑1 (MQ=255)                                                                                                                            tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCaaat                >  1:311228/1‑139 (MQ=21)                                                                                                                            tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                >  1:143815/1‑141 (MQ=35)                                                                                                                            tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTCCATACAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCAACt                >  1:164453/1‑141 (MQ=21)                                                                                                                                atatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAATGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTGAAAATTATACTATAAATTCAACTTTGc            >  1:49910/1‑141 (MQ=25)                                                                                                                                   tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGCATACAGTGATAACATCCTTGTTAAATTGCTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTCCaaa         >  2:172205/1‑140 (MQ=255)                                                                                                                                      cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa      >  1:92812/1‑141 (MQ=35)                                                                                                                                          aTAACTCAATGTATTTCTTAGATAGTAGCCAACGGCGACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCg  <  2:98122/141‑1 (MQ=25)                                                                                                                                             |                                                                                                                                          CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCG  >  CP000730/36331‑36608

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: ‑