Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	503,672	T→C	E342G (GAG→GGG)	fsr ←	fosmidomycin efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	503,672	0	T	C	98.4%	208.5 / NA	63	E342G (GAG→GGG)	fsr	fosmidomycin efflux pump
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/32);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCC  >  NC_000913/503537‑503807                                                                                                                                        |                                                                                                                                        gCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTg                                                                                                                                     >  1:216899/1‑140 (MQ=255)     gacagttATAGACTACCTCGATTCTTGTGCGATCGGCAATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCg                                                                                                                                 <  1:548477/135‑1 (MQ=255)     aGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCg                                                                                                                                 <  1:39444/140‑1 (MQ=255)     aGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCg                                                                                                                                 <  1:292537/140‑1 (MQ=255)     aGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCg                                                                                                                                 <  2:339707/140‑1 (MQ=255)       aTTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTa                                                                                                                               >  2:1036350/1‑140 (MQ=255)       aTTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTa                                                                                                                               >  2:1086018/1‑140 (MQ=255)              aGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGa                                                                                                                        <  1:857368/140‑1 (MQ=255)                 cTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATg                                                                                                                     >  2:1230098/1‑140 (MQ=255)                 cTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATg                                                                                                                     >  1:1387508/1‑140 (MQ=255)                   aaCTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGc                                                                                                                   >  2:267044/1‑140 (MQ=255)                   aaCTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGc                                                                                                                   >  2:379550/1‑140 (MQ=255)                    aCTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCa                                                                                                                  >  2:131045/1‑140 (MQ=255)                    aCTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCa                                                                                                                  <  1:883982/140‑1 (MQ=255)                    aCTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACAGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCa                                                                                                                  <  1:533273/140‑1 (MQ=255)                     cTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCGTACCGATATGTCGTGGAAGCAGCCCCTGAGGGTAGACCAGAAtctgtc                                                                                                                 >  2:548477/1‑135 (MQ=255)                          tGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGGAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAAt                                                                                                            >  1:1370508/1‑140 (MQ=255)                          tGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCACCTAAGCGTAGACAAGAATGGCAGAGAAt                                                                                                            >  1:639420/1‑140 (MQ=255)                          tGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGAAAGCAGCTCCTAAGCGTAGACAAGAATGGCAGAGAAt                                                                                                            >  1:635146/1‑140 (MQ=255)                             tGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGcc                                                                                                         <  2:349616/140‑1 (MQ=255)                                  tgATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGc                                                                                                    <  1:302735/140‑1 (MQ=255)                                       ggCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGa                                                                                               >  2:204278/1‑140 (MQ=255)                                       ggCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGa                                                                                               <  1:267044/140‑1 (MQ=255)                                          gATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATaa                                                                                            <  2:1174907/140‑1 (MQ=255)                                            tAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAAt                                                                                          >  1:60300/1‑140 (MQ=255)                                            tAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAAt                                                                                          <  1:1230098/140‑1 (MQ=255)                                              aGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATcc                                                                                        >  2:307190/1‑140 (MQ=255)                                                     gAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCATAACGGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATc                                                                                 <  2:558238/140‑1 (MQ=255)                                                     gAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATc                                                                                 >  2:1015742/1‑140 (MQ=255)                                                        cTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCAcc                                                                              <  1:1129468/140‑1 (MQ=255)                                                          gCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGt                                                                            >  1:974755/1‑140 (MQ=255)                                                          gCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGt                                                                            >  1:44139/1‑140 (MQ=255)                                                             gCTCCCAGACCTCCCATGCCAATAGCAAAACGGAACAAGGGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTaa                                                                         <  2:635146/140‑1 (MQ=255)                                                                  cAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAAccc                                                                    <  1:1086018/140‑1 (MQ=255)                                                                  cAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAAccc                                                                    <  1:131045/140‑1 (MQ=255)                                                                      ccTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGt                                                                <  1:204278/140‑1 (MQ=255)                                                                             tGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGAATATATCGAATAATCAGCGTTGAAACCCCCGTCCAGTGc                                                         >  2:506832/1‑140 (MQ=255)                                                                               ccAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCa                                                        <  1:1535232/139‑1 (MQ=255)                                                                               ccAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCa                                                        >  2:1535232/1‑139 (MQ=255)                                                                                cAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAgg                                                      <  1:1036350/140‑1 (MQ=255)                                                                                 aaaaGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATCAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGc                                                     <  2:974755/140‑1 (MQ=255)                                                                                       aaaaCCGAAAAAGAGTCCAGAATCCATACCGTTACCTCCTGGAAGCAGTCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGt                                               <  2:639420/140‑1 (MQ=255)                                                                                       aaaaCCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGt                                               <  1:933347/140‑1 (MQ=255)                                                                                         aaCCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAg                                             >  1:1297430/1‑140 (MQ=255)                                                                                         aaCCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAAACACCGTTAAAAGCCCCGTCCATTGCAGGCTGGCGTAg                                             >  1:618064/1‑140 (MQ=255)                                                                                                aaaGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTaaaa                                      <  1:902569/140‑1 (MQ=255)                                                                                                      tCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACGCCCGTCCAGTGTAGTCTGGCGGAGGGTAAAATCAGCg                                >  1:489386/1‑140 (MQ=255)                                                                                                      tCCAGAAACCATACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCg                                >  2:266862/1‑140 (MQ=255)                                                                                                               cATACCGATACGTCCAGGAAGCGGCGCCTGAGCGTAGTCCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGgcg                       <  1:506832/140‑1 (MQ=255)                                                                                                                 tACCGATACGTCTTGGAAGCAGCCCCCGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCa                     <  1:924046/140‑1 (MQ=255)                                                                                                                 tACCGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCa                     <  2:347866/140‑1 (MQ=255)                                                                                                                   ccGATACGTCCTGGACGCAGCCTCTCAGCCTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAAc                   <  1:476549/140‑1 (MQ=255)                                                                                                                   ccGATACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAAc                   <  2:60300/140‑1 (MQ=255)                                                                                                                      aTACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGcc                <  1:379550/140‑1 (MQ=255)                                                                                                                       tACGTCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCGGGCGTAGTGTAAAATCAGCGTAAACGGCGCAACGCCg               >  1:469310/1‑140 (MQ=255)                                                                                                                           tCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGa           <  2:317118/140‑1 (MQ=255)                                                                                                                           tCCTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGa           >  1:317118/1‑140 (MQ=255)                                                                                                                            ccTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGAt          >  1:291140/1‑140 (MQ=255)                                                                                                                            ccTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGAt          >  2:918405/1‑140 (MQ=255)                                                                                                                            ccTGGAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGGAAGATCTGCATAAACGACGCAACGCCGAGGAt          >  1:722930/1‑140 (MQ=255)                                                                                                                             cTGGAAGCAGCCCCTGAACGTAAACCAGAATGGCAGAGAGTGCCGAAGTGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATa         <  2:618064/140‑1 (MQ=255)                                                                                                                                gAAGCAGCCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATagag      >  1:1414754/1‑140 (MQ=255)                                                                                                                                    cagcCCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGcccc  <  2:1420751/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCC  >  NC_000913/503537‑503807

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