Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,373,558	A→T	M73K (ATG→AAG)	pbuG ←	NCS2 family nucleobase:cation symporter‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,373,558	0	A	T	100.0%	244.9 / NA	68	M73K (ATG→AAG)	pbuG	NCS2 family nucleobase:cation symporter‑2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (28/40);  total (28/40)

AGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCAT  >  CP000730/2373421‑2373694                                                                                                                                          |                                                                                                                                         aGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGa                                                                                                                                       >  2:783809/1‑141 (MQ=255)      ataCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGc                                                                                                                                  >  1:269621/1‑141 (MQ=255)       taCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCg                                                                                                                                 >  2:351381/1‑141 (MQ=255)          cTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGc                                                                                                                              <  1:533966/141‑1 (MQ=255)                  aaaaCGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAAt                                                                                                                      <  2:561757/141‑1 (MQ=255)                   aaaCGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAgatacattcct                                                                                                                     >  2:108680/1‑131 (MQ=255)                       gCCGGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAAtgctgc                                                                                                                 <  2:759592/141‑1 (MQ=255)                          tGTTAAACCAACTTGCCAAGTAATACCCATGGTTAACACAACAGTAAATACAAATTATGCATTTAACCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTaa                                                                                                              >  1:383250/1‑141 (MQ=255)                          tGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTaa                                                                                                              <  1:789078/141‑1 (MQ=255)                              aaaCCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCa                                                                                                          >  1:524063/1‑141 (MQ=255)                               aaCCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAg                                                                                                         <  1:370156/141‑1 (MQ=255)                                  cAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCg                                                                                                      >  1:227069/1‑141 (MQ=255)                                      ttGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTAc                                                                                                  <  1:315867/141‑1 (MQ=255)                                         ccAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaa                                                                                               >  2:619627/1‑141 (MQ=255)                                          cAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaaa                                                                                              <  1:330756/141‑1 (MQ=255)                                          cAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaaa                                                                                              <  1:137633/141‑1 (MQ=255)                                          cAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaaaa                                                                                              <  1:687998/141‑1 (MQ=255)                                            aGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAAt                                                                                            <  1:115265/141‑1 (MQ=255)                                            aGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAAt                                                                                            >  1:115602/1‑141 (MQ=255)                                                 tACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCAc                                                                                       >  1:770520/1‑141 (MQ=255)                                                   cccATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCt                                                                                     >  1:742510/1‑141 (MQ=255)                                                         gTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATggaccaaggtgcc                                                                               <  1:648424/141‑12 (MQ=255)                                                          ttAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATggaccaaggtgcca                                                                              <  2:712817/141‑13 (MQ=255)                                                          ttAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCa                                                                              <  1:669151/141‑1 (MQ=255)                                                           tAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATTCCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                                                                             <  1:467037/141‑1 (MQ=255)                                                             acacAACAGTAAATGCAAAGAATGCATTTAATCCCATGCCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCTAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                                                                           >  1:572262/1‑141 (MQ=255)                                                             acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                                                                           >  1:115105/1‑141 (MQ=255)                                                             acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                                                                           >  1:272089/1‑141 (MQ=255)                                                              cacaACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAtttt                                                                          <  2:173100/141‑1 (MQ=255)                                                               acaacaGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTc                                                                         <  1:731059/141‑1 (MQ=255)                                                                  acaGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCata                                                                      >  1:661243/1‑141 (MQ=255)                                                                        aaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttct                                                                <  2:783366/141‑1 (MQ=255)                                                                         aTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttct                                                                <  2:647695/140‑1 (MQ=255)                                                                         aTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttct                                                                >  1:647695/1‑140 (MQ=255)                                                                               aGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc                                                         <  1:150799/141‑1 (MQ=255)                                                                                  aTGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTc                                                                         <  2:468885/122‑1 (MQ=255)                                                                                  aTGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTc                                                                         >  1:468885/1‑122 (MQ=255)                                                                                       tttAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAaca                                                 <  2:115105/141‑1 (MQ=255)                                                                                        ttAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAacac                                                <  2:131196/141‑1 (MQ=255)                                                                                        ttAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAacac                                                <  2:524063/141‑1 (MQ=255)                                                                                         tAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACAcc                                               >  2:21964/1‑141 (MQ=255)                                                                                            tgcgatgccTGGTGCTAAGGCAATGGGATAGTGAGCGATGAGGCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGc                                            <  2:758859/134‑1 (MQ=255)                                                                                             cccATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCa                                                                              <  2:568061/106‑1 (MQ=255)                                                                                             cccATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCa                                                                              >  1:568061/1‑106 (MQ=255)                                                                                               cATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTaa                                         >  1:680014/1‑141 (MQ=255)                                                                                               cATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGAAACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTaa                                         <  2:416887/141‑1 (MQ=255)                                                                                                     tGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaa                                   <  1:344391/141‑1 (MQ=255)                                                                                                     tGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaa                                   >  1:495277/1‑141 (MQ=255)                                                                                                     tGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaa                                   <  2:770520/141‑1 (MQ=255)                                                                                                     tGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaa                                   <  1:230583/141‑1 (MQ=255)                                                                                                       gTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaaaa                                 >  2:413721/1‑141 (MQ=255)                                                                                                       gTGCTAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTaaaa                                 >  2:770448/1‑141 (MQ=255)                                                                                                           tAACGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTg                             <  2:383250/141‑1 (MQ=255)                                                                                                             aCGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCg                           <  1:546327/141‑1 (MQ=255)                                                                                                              cGCAATTGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCgg                          >  2:624325/1‑141 (MQ=255)                                                                                                                    tGGATATTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCAACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAAc                    <  1:32854/141‑1 (MQ=255)                                                                                                                       atatTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGc                 <  2:742510/141‑1 (MQ=255)                                                                                                                       atatTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGc                 <  2:719140/141‑1 (MQ=255)                                                                                                                         atTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTCCATTGATTGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTa               <  1:14468/141‑1 (MQ=255)                                                                                                                         atTTAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCAGAGCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTa               <  1:14544/141‑1 (MQ=255)                                                                                                                          tttAGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTaa              >  1:648973/1‑141 (MQ=255)                                                                                                                             aGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAat           <  2:306042/141‑1 (MQ=255)                                                                                                                             aGCTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAat           <  1:588379/141‑1 (MQ=255)                                                                                                                               cTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAatat         >  2:759726/1‑141 (MQ=255)                                                                                                                               cTATTAGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAatat         >  2:169109/1‑141 (MQ=255)                                                                                                                                    aGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGcc    <  1:503451/141‑1 (MQ=255)                                                                                                                                    aGCCCCTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGcc    >  2:503451/1‑141 (MQ=255)                                                                                                                                      ccccTTGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCAt  <  1:181399/141‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         AGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTAACACCTGCTAAACTTAAAACTTGCGGGTTAACTGCTAAAATATAGGCCAT  >  CP000730/2373421‑2373694

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑