Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,455,998	A→C	T72P (ACC→CCC)	mpl →	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,455,998	0	A	C	100.0%	241.3 / NA	65	T72P (ACC→CCC)	mpl	UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (37/28);  total (37/28)

TAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCAC  >  NC_000913/4455864‑4456134                                                                                                                                       |                                                                                                                                         tAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGt                                                                                                                                     >  1:40641/1‑140 (MQ=255)          cGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCg                                                                                                                            >  2:452778/1‑140 (MQ=255)           ggACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCgt                                                                                                                           >  2:1198759/1‑140 (MQ=255)               gCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCgtgtgt                                                                                                                       <  1:619208/140‑1 (MQ=255)                         aTCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTCCGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTa                                                                                                             <  2:37946/140‑1 (MQ=255)                          tCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTAc                                                                                                            <  2:687180/140‑1 (MQ=255)                            cgccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTg                                                                                                          <  1:1198759/140‑1 (MQ=255)                             gccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTgg                                                                                                                     <  1:753594/128‑1 (MQ=255)                             gccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTgg                                                                                                                     >  2:753594/1‑128 (MQ=255)                             gccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTgg                                                                                                         >  1:1053029/1‑140 (MQ=255)                              ccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGa                                                                                                        <  2:978890/140‑1 (MQ=255)                               cgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaa                                                                                                       <  2:355781/140‑1 (MQ=255)                               cgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaa                                                                                                       >  2:747902/1‑140 (MQ=255)                                  tGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaaaaa                                                                                                    >  2:746179/1‑140 (MQ=255)                                    aGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaaaaaaa                                                                                                  <  2:226271/140‑1 (MQ=255)                                      ctccTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAata                                                                                                <  2:617874/138‑1 (MQ=255)                                      cACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAata                                                                                                <  2:40641/140‑1 (MQ=255)                                      cACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAata                                                                                                >  2:15676/1‑140 (MQ=255)                                        ccTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATc                                                                                              <  1:1213300/140‑1 (MQ=255)                                        ccTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATc                                                                                              <  1:15676/140‑1 (MQ=255)                                         cTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATcc                                                                                             >  1:720289/1‑140 (MQ=255)                                           tACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCt                                                                                           >  2:48326/1‑140 (MQ=255)                                           tACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCt                                                                                           >  2:24893/1‑140 (MQ=255)                                                gagaAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTatat                                                                                      <  2:720289/140‑1 (MQ=255)                                                  gaAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    <  1:747902/140‑1 (MQ=255)                                                   aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    <  1:258171/139‑1 (MQ=255)                                                   aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    >  1:378729/1‑139 (MQ=255)                                                   aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    <  2:378729/139‑1 (MQ=255)                                                   aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    >  2:258171/1‑139 (MQ=255)                                                   aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTc                                                                                   >  2:744130/1‑140 (MQ=255)                                                   aaGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTc                                                                                   <  1:746179/140‑1 (MQ=255)                                                      cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTc                                                                                   >  1:1083768/1‑137 (MQ=255)                                                      cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTc                                                                                   <  2:1083768/137‑1 (MQ=255)                                                      cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAgg                                                                                >  2:370811/1‑140 (MQ=255)                                                      cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATACCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    <  2:665562/136‑1 (MQ=255)                                                      cAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATACCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    >  1:665562/1‑136 (MQ=255)                                                       aaGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGt                                                                               >  2:979886/1‑140 (MQ=255)                                                           cATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCcac                                                                           <  2:1040628/140‑1 (MQ=255)                                                                aGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTgg                                                                      <  1:452778/140‑1 (MQ=255)                                                                    gATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    >  2:13596/1‑122 (MQ=255)                                                                    gATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGt                                                                                    <  1:13596/122‑1 (MQ=255)                                                                     aTTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCa                                                                 >  2:44885/1‑140 (MQ=255)                                                                               aCGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTtgtg                                                       >  1:517625/1‑140 (MQ=255)                                                                                cGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGtgc                                                      >  2:17279/1‑140 (MQ=255)                                                                                   tGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGtgctgc                                                   >  2:718325/1‑140 (MQ=255)                                                                                   tGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGtgctgc                                                   >  1:998941/1‑140 (MQ=255)                                                                                    gccagccaGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTgcg                                                  <  2:517625/140‑1 (MQ=255)                                                                                        gccagcTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGAc                                              >  2:377262/1‑140 (MQ=255)                                                                                        gccagcTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGAc                                              >  1:940197/1‑140 (MQ=255)                                                                                            gcTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCt                                          <  2:1008138/140‑1 (MQ=255)                                                                                             cTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTg                                         >  2:680274/1‑140 (MQ=255)                                                                                             cTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTg                                         >  2:942216/1‑140 (MQ=255)                                                                                                       aGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGcc                               >  1:403765/1‑140 (MQ=255)                                                                                                       aGCCGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGcc                               <  1:17279/140‑1 (MQ=255)                                                                                                         ccGGATCTGGTGATTATTGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGt                             <  1:718325/140‑1 (MQ=255)                                                                                                                         ttGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGc             <  1:744130/140‑1 (MQ=255)                                                                                                                          tGGCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCa            >  1:954581/1‑140 (MQ=255)                                                                                                                           ggCAACGCATTGCCCCGTGGAAAGCCGTGTGTGGAAGCGGTACTGTACAAAAATATCCCTTAGATTTCAGGTCCACAGGGGCTCCACGATTTTCTGCTGCTCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCaa           <  2:724926/140‑1 (MQ=255)                                                                                                                            gCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCaaa          >  1:724926/1‑140 (MQ=255)                                                                                                                            gCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCaaa          >  1:833444/1‑140 (MQ=255)                                                                                                                            gCAACGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCaaa          >  1:1070955/1‑140 (MQ=255)                                                                                                                              aaCGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAac        >  1:1126458/1‑140 (MQ=255)                                                                                                                               aCGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAacc       >  2:680228/1‑140 (MQ=255)                                                                                                                                cGCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAacca      >  1:271896/1‑140 (MQ=255)                                                                                                                                 gCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAaccac     <  1:212113/140‑1 (MQ=255)                                                                                                                                 gCCATGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAaccac     >  2:257436/1‑140 (MQ=255)                                                                                                                                    aTGCCCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAaccaccac  <  1:388726/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         TAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTTATATGTCAGGTCCACAGTGGCTGCACGATTTTGTGCTGCGCGACCGCTGGGTGCTGGCCGTTGCCGGTACACACGGCAAAACCACCAC  >  NC_000913/4455864‑4456134

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