Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,152	A→G	I5T (ATC→ACC)	phnC ←	phosphonate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,152	0	A	G	100.0%	295.0 / NA	81	I5T (ATC→ACC)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/44);  total (37/44)

GACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGG  >  NC_000913/4325021‑4325276                                                                                                                                    |                                                                                                                             gACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                       <  2:957586/139‑1 (MQ=255) gACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                       <  2:885675/139‑1 (MQ=255)    gTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTg                                                                                                                    >  1:411946/1‑139 (MQ=255)      aaaaGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCa                                                                                                                  <  1:1072503/139‑1 (MQ=255)       aaaGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCat                                                                                                                 <  2:528346/139‑1 (MQ=255)             tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                           >  2:1090530/1‑139 (MQ=255)               gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCGGATCAACCGCATGCAGCGCATGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         <  1:40772/139‑1 (MQ=255)                   ttCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTcc                                                                                                     <  2:1561027/139‑1 (MQ=255)                         aaCCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                                   <  2:8667/119‑1 (MQ=255)                         aaCCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                                   >  1:8667/1‑119 (MQ=255)                          aCCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCACGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         <  2:1502159/128‑1 (MQ=255)                          aCCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCACGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         >  1:1502159/1‑128 (MQ=255)                                 ggCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCt                                                                                       <  1:718498/139‑1 (MQ=255)                                   cccAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCa                                                                                                                  >  2:580671/1‑110 (MQ=255)                                   cccAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCa                                                                                                                  <  1:580671/110‑1 (MQ=255)                                     cAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                                                                                   >  1:1377805/1‑139 (MQ=255)                                        gCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                                                                                >  1:1103569/1‑139 (MQ=255)                                         cAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGCTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAACAAAGTGGCTCCATCGTg                                                                               <  1:955502/139‑1 (MQ=255)                                            agCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                                   >  1:969242/1‑100 (MQ=255)                                            agCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                                   <  2:969242/100‑1 (MQ=255)                                            agCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg                                                                            <  1:1090530/139‑1 (MQ=255)                                              ccaccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                            >  1:285647/1‑89 (MQ=255)                                                          cGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                       >  1:178302/1‑82 (MQ=255)                                                          cGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                       <  2:178302/82‑1 (MQ=255)                                                            tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                            <  2:442423/75‑1 (MQ=255)                                                            tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                            >  1:442423/1‑75 (MQ=255)                                                              atgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTgg                                                                                         <  2:145207/108‑1 (MQ=255)                                                              atgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTgg                                                                                         >  1:145207/1‑108 (MQ=255)                                                               tgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                         <  2:1377805/139‑1 (MQ=255)                                                                      tCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         >  1:1014125/1‑84 (MQ=255)                                                                      tCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         <  2:1014125/84‑1 (MQ=255)                                                                         gATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGG‑GGATTATTGTGACGAt                                                             <  2:1632851/124‑1 (MQ=255)                                                                         gATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGG‑GGATTATTGTGACGAt                                                             >  1:1632851/1‑124 (MQ=255)                                                                          aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                                                                               >  1:354231/1‑106 (MQ=255)                                                                          aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                                                                               <  2:354231/106‑1 (MQ=255)                                                                           tCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGt                                             >  2:60493/1‑139 (MQ=255)                                                                           tCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGt                                             >  1:1486825/1‑139 (MQ=255)                                                                             aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                                       <  1:837669/79‑1 (MQ=255)                                                                             aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                                       >  2:837669/1‑79 (MQ=255)                                                                             aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                                                                                   >  2:590417/1‑99 (MQ=255)                                                                             aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                                                                                   <  1:590417/99‑1 (MQ=255)                                                                             aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGa                                           >  2:806095/1‑139 (MQ=255)                                                                             aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGa                                           >  2:1294168/1‑139 (MQ=255)                                                                               ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg                                                                            >  1:1262692/1‑104 (MQ=255)                                                                               ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg                                                                            <  2:1262692/104‑1 (MQ=255)                                                                                  cATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAAc                                      <  2:1486825/139‑1 (MQ=255)                                                                                    tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAAt                                                                                                                          <  1:797985/53‑1 (MQ=255)                                                                                    tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAAt                                                                                                                          >  2:797985/1‑53 (MQ=255)                                                                                    tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGtt                                                                                                                      <  2:318630/57‑1 (MQ=255)                                                                                    tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGtt                                                                                                                      >  1:318630/1‑57 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         <  2:117903/65‑1 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         >  1:117903/1‑65 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                         >  2:855168/1‑113 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                         <  1:855168/113‑1 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAAc                                      <  2:4582/132‑1 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAAc                                      >  1:4582/1‑132 (MQ=255)                                                                                         cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                               >  1:852831/1‑139 (MQ=255)                                                                                                     aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaattaa                                                 >  2:1151502/1‑109 (MQ=255)                                                                                                     aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaattaa                                                 <  1:1151502/109‑1 (MQ=255)                                                                                                     aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGa                   >  1:1524437/1‑139 (MQ=255)                                                                                                      ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         <  2:176185/52‑1 (MQ=255)                                                                                                      ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         >  1:176185/1‑52 (MQ=255)                                                                                                      ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGaa                  <  2:1524437/139‑1 (MQ=255)                                                                                                      ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGaa                  <  2:852831/139‑1 (MQ=255)                                                                                                         aaGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAt                                                             >  2:1639456/1‑93 (MQ=255)                                                                                                         aaGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAt                                                             <  1:1639456/93‑1 (MQ=255)                                                                                                              tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa          <  1:1294168/139‑1 (MQ=255)                                                                                                              tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa          <  2:411946/139‑1 (MQ=255)                                                                                                              tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa          <  2:1103569/139‑1 (MQ=255)                                                                                                              tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa          <  1:806095/139‑1 (MQ=255)                                                                                                              tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa          <  1:60493/139‑1 (MQ=255)                                                                                                               ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         <  1:1417305/43‑1 (MQ=255)                                                                                                               ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                         >  2:1417305/1‑43 (MQ=255)                                                                                                               ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg                                                                            <  2:184716/72‑1 (MQ=255)                                                                                                               ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg                                                                            >  1:184716/1‑72 (MQ=255)                                                                                                               ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa                                                                 <  2:1650186/83‑1 (MQ=255)                                                                                                               ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa                                                                 >  1:1650186/1‑83 (MQ=255)                                                                                                                      gCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                         >  2:1157165/1‑84 (MQ=255)                                                                                                                      gCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                         <  1:1157165/84‑1 (MQ=255)                                                                                                                      gCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCgg  >  2:159272/1‑139 (MQ=255)                                                                                                                       cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAg                                              >  2:522018/1‑94 (MQ=255)                                                                                                                       cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAg                                              <  1:522018/94‑1 (MQ=255)                                                                                                                                    |                                                                                                                             GACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGG  >  NC_000913/4325021‑4325276

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