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%	341.0 / NA	94	I5T (ATC→ACC)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (46/48);  total (46/48)

TAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGAC  >  NC_000913/4325016‑4325269                                                                                                                                         |                                                                                                                      tAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCGAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                     <  1:1234442/139‑1 (MQ=255)      gACGTAAAAGGGTGGATTTTCCCGAACCCGACGGCCCAAGCAGTGCCACCATTTCACCGTGATTAAGGTTCAGATCAACCGCATTCAGCGCCTGATGCTGCTTGAAGGTTTTGGCGATCTTCTCTACACGGGTAACCGt                                                                                                                <  1:1572882/139‑1 (MQ=255)          tAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                            >  2:804487/1‑139 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACTGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                    >  2:1283592/1‑139 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                     <  1:2001395/122‑1 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTa                                                                                                                     >  2:2001395/1‑122 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                >  2:1205300/1‑127 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                    >  1:1053398/1‑139 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                    >  1:1510061/1‑139 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                    >  2:1123352/1‑139 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCa                                                                                                    >  2:723807/1‑139 (MQ=255)                  tGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCAAGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                <  1:1205300/127‑1 (MQ=255)                   ggATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                <  1:129813/126‑1 (MQ=255)                   ggATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                >  2:129813/1‑126 (MQ=255)                   ggATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAg                                                                                                   <  1:182973/139‑1 (MQ=255)                    gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  2:711344/139‑1 (MQ=255)                    gATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  1:836422/139‑1 (MQ=255)                          ccGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  2:1160226/1‑133 (MQ=255)                          ccGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  1:1160226/133‑1 (MQ=255)                                     cGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGc                                                                                 >  1:988955/1‑139 (MQ=255)                                        cccAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCg                                                                                                                 >  2:999484/1‑104 (MQ=255)                                        cccAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCg                                                                                                                 <  1:999484/104‑1 (MQ=255)                                        cccAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTcc                                                                              >  2:687216/1‑139 (MQ=255)                                           aaGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCT‑AAAAAGTGGCTCCATCg                                                                          <  2:1053398/139‑1 (MQ=255)                                              cAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                                                                        <  1:423931/139‑1 (MQ=255)                                                  gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                                >  2:1169713/1‑111 (MQ=255)                                                  gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                                <  1:1169713/111‑1 (MQ=255)                                                  gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCgg                                                                    >  1:1841486/1‑139 (MQ=255)                                                       caTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTTAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCa                                                                                                           <  2:1087185/95‑1 (MQ=255)                                                       caTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCa                                                                                                           >  1:1087185/1‑95 (MQ=255)                                                       caTTACACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGATTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGAGGCTCCATCGTGGCGGattat                                                               >  1:1197793/1‑139 (MQ=255)                                                           tCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg                                                           <  2:988955/139‑1 (MQ=255)                                                                 tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                            <  1:1245823/84‑1 (MQ=255)                                                                 tgatgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGc                                                                                                            >  2:1245823/1‑84 (MQ=255)                                                                    tgaATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                  <  1:804487/139‑1 (MQ=255)                                                                            cAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGttaattaa                                          <  1:354261/139‑1 (MQ=255)                                                                              gATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTaaaa                                        <  1:687216/139‑1 (MQ=255)                                                                               aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAg                                       >  1:306084/1‑139 (MQ=255)                                                                                tCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGt                                      <  2:1841486/139‑1 (MQ=255)                                                                                  aaCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGa                                    >  2:222411/1‑139 (MQ=255)                                                                                   aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGAc                                   >  2:499820/1‑139 (MQ=255)                                                                                    ccGCATGCAGCGCCTGATGGTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTa                                                                                           >  2:1370254/1‑82 (MQ=255)                                                                                    ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTa                                                                                           <  1:1370254/82‑1 (MQ=255)                                                                                    ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                  >  1:410882/1‑139 (MQ=255)                                                                                    ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                  >  1:1680254/1‑139 (MQ=255)                                                                                    ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                  <  1:1123352/139‑1 (MQ=255)                                                                                      gCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                <  2:1510061/139‑1 (MQ=255)                                                                                       cATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                   >  2:1426888/1‑103 (MQ=255)                                                                                       cATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                   <  1:1426888/103‑1 (MQ=255)                                                                                         tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  2:206696/70‑1 (MQ=255)                                                                                         tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  1:206696/1‑70 (MQ=255)                                                                                         tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAg                             >  1:1915739/1‑139 (MQ=255)                                                                                         tGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAg                             >  1:1448062/1‑139 (MQ=255)                                                                                            aGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  1:1006258/67‑1 (MQ=255)                                                                                            aGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  2:1006258/1‑67 (MQ=255)                                                                                              cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGGCGATTTGGTTAATTAAAAGTGACTAACAGATGaa                        >  2:685719/1‑139 (MQ=255)                                                                                              cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                >  2:934206/1‑131 (MQ=255)                                                                                              cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                <  1:934206/131‑1 (MQ=255)                                                                                              cgcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGaa                        >  2:802974/1‑139 (MQ=255)                                                                                               gcCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAag                       >  2:10406/1‑139 (MQ=255)                                                                                                    aTGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                                                                        >  2:1177334/1‑85 (MQ=255)                                                                                                    aTGCTGATTGAAGCTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAATAAAGTGGCTCCATCGTg                                                                        <  1:1177334/85‑1 (MQ=255)                                                                                                        tGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACgg              >  1:333702/1‑139 (MQ=255)                                                                                                         gATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGGCGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg             <  1:685719/139‑1 (MQ=255)                                                                                                         gATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg             <  1:802974/139‑1 (MQ=255)                                                                                                          aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  1:1983562/53‑1 (MQ=255)                                                                                                          aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  1:326682/1‑53 (MQ=255)                                                                                                          aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  2:326682/53‑1 (MQ=255)                                                                                                          aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  1:1635854/53‑1 (MQ=255)                                                                                                          aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  2:1983562/1‑53 (MQ=255)                                                                                                          aTTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  2:1635854/1‑53 (MQ=255)                                                                                                           ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                  <  1:68847/100‑1 (MQ=255)                                                                                                           ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTgg                                                  >  2:68847/1‑100 (MQ=255)                                                                                                           ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGaa           <  1:723807/139‑1 (MQ=255)                                                                                                           ttGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGaa           <  2:410882/139‑1 (MQ=255)                                                                                                             gAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGCTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAAtt         <  2:1252996/139‑1 (MQ=255)                                                                                                              aaGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTc        >  2:1670732/1‑139 (MQ=255)                                                                                                               aGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt                                                            <  2:410800/86‑1 (MQ=255)                                                                                                               aGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt                                                            >  1:410800/1‑86 (MQ=255)                                                                                                               aGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCg       >  2:1135093/1‑139 (MQ=255)                                                                                                               aGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCg       <  2:333702/139‑1 (MQ=255)                                                                                                                 gTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGAt     <  2:1915739/139‑1 (MQ=255)                                                                                                                 gTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGAt     <  1:499820/139‑1 (MQ=255)                                                                                                                 gTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGAt     <  1:1670732/139‑1 (MQ=255)                                                                                                                   tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa   <  2:306084/139‑1 (MQ=255)                                                                                                                   tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa   <  1:222411/139‑1 (MQ=255)                                                                                                                   tttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGa   <  2:1448062/139‑1 (MQ=255)                                                                                                                    ttGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGAc  >  1:1377093/1‑139 (MQ=255)                                                                                                                          aGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  <  1:284310/37‑1 (MQ=255)                                                                                                                          aGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGc                                                                                                  >  2:284310/1‑37 (MQ=255)                                                                                                                            cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                   <  1:268635/66‑1 (MQ=255)                                                                                                                            cTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                   >  2:268635/1‑66 (MQ=255)                                                                                                                                  gACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                   <  1:337407/60‑1 (MQ=255)                                                                                                                                  gACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa                                                                   >  2:337407/1‑60 (MQ=255)                                                                                                                                         |                                                                                                                      TAAGTGACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGAC  >  NC_000913/4325016‑4325269

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