Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,599,749	T→C	G4G (GGA→GGG)	dapA ←	4‑hydroxy‑tetrahydrodipicolinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,599,749	0	T	C	100.0%	258.2 / NA	69	G4G (GGA→GGG)	dapA	4‑hydroxy‑tetrahydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/36);  total (33/36)

GCGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAC  >  NC_000913/2599615‑2599885                                                                                                                                       |                                                                                                                                         gCGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGt                                                                                                                                     >  2:505685/1‑140 (MQ=255) gCGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGt                                                                                                                                     >  1:730575/1‑140 (MQ=255) gCGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGt                                                                                                                                     >  2:596554/1‑140 (MQ=255)  cGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTg                                                                                                                                    >  1:415549/1‑140 (MQ=255)  cGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTg                                                                                                                                    >  1:877284/1‑140 (MQ=255)     aCTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAAc                                                                                                                                 <  1:163149/140‑1 (MQ=255)       tCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACAt                                                                                                                               <  2:322851/140‑1 (MQ=255)            aGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTGCCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGcc                                                                                                                          <  2:877284/140‑1 (MQ=255)            aGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGcc                                                                                                                          <  2:866299/140‑1 (MQ=255)                gtgCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATcc                                                                                                                      <  1:829601/140‑1 (MQ=255)                 tgCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCct                                                                                                                     <  2:229290/140‑1 (MQ=255)                          aaaCGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCaaa                                                                                                            >  1:939600/1‑140 (MQ=255)                          aaaCGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCaaa                                                                                                            >  1:164475/1‑140 (MQ=255)                           aaCGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCct                                                                                                                     <  1:852643/130‑1 (MQ=255)                           aaCGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCct                                                                                                                     >  2:852643/1‑130 (MQ=255)                            aCGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACa                                                                                                          <  2:358085/140‑1 (MQ=255)                                 cGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAgtgt                                                                                                     <  2:164475/140‑1 (MQ=255)                                  gCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGtc                                                                                                    >  2:895526/1‑140 (MQ=255)                                        gTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATg                                                                                              >  1:788480/1‑140 (MQ=255)                                          aCCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGt                                                                                            <  2:415549/140‑1 (MQ=255)                                          aCCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGt                                                                                            <  1:333675/140‑1 (MQ=255)                                            cGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGt                                                                                            >  2:761150/1‑138 (MQ=255)                                            cGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGt                                                                                            <  1:761150/138‑1 (MQ=255)                                              cTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGt                                                                                            >  2:925143/1‑136 (MQ=255)                                              cTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGt                                                                                            <  1:925143/136‑1 (MQ=255)                                               tGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGtt                                                                                       <  2:788480/140‑1 (MQ=255)                                               tGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGtt                                                                                       <  2:59864/140‑1 (MQ=255)                                                  cGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTgg                                                                                    >  2:588032/1‑140 (MQ=255)                                                   gACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACTTTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGt                                                                                   <  2:48393/140‑1 (MQ=255)                                                     cATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTAt                                                                                 <  2:898927/140‑1 (MQ=255)                                                      aTGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATg                                                                                <  2:730575/140‑1 (MQ=255)                                                           aatcaatcAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCAtt                                                                           >  2:700995/1‑140 (MQ=255)                                                                atcaGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaag                                                                      >  2:40073/1‑140 (MQ=255)                                                                   aGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaa                                                                   <  1:504303/140‑1 (MQ=255)                                                                   aGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaa                                                                   <  2:939600/140‑1 (MQ=255)                                                                     ttttttCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaagc                                                                 <  1:895526/140‑1 (MQ=255)                                                                     ttttttCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAaagcaagc                                                                 >  2:282952/1‑140 (MQ=255)                                                                             aGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGaa                                                         <  1:40073/140‑1 (MQ=255)                                                                                  gCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGtt                                                    >  2:683225/1‑140 (MQ=255)                                                                                   cccGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTc                                                   <  1:406619/140‑1 (MQ=255)                                                                                      gacagacaTTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGa                                                <  1:904892/140‑1 (MQ=255)                                                                                        cagacaTTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGt                                                     <  2:293066/133‑1 (MQ=255)                                                                                        cagacaTTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGt                                                     >  1:293066/1‑133 (MQ=255)                                                                                           acaTTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgtt                                           <  1:700995/140‑1 (MQ=255)                                                                                                aCCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATg                                      <  2:602823/140‑1 (MQ=255)                                                                                                      ttCATCCATCGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                >  1:711830/1‑140 (MQ=255)                                                                                                         atccatcGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttt                             >  2:871039/1‑140 (MQ=255)                                                                                                           ccatcGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttttt                           >  1:1094091/1‑140 (MQ=255)                                                                                                           ccatcGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttttt                           >  1:231297/1‑140 (MQ=255)                                                                                                               cGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                       >  2:135832/1‑140 (MQ=255)                                                                                                               cGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                       <  1:47271/140‑1 (MQ=255)                                                                                                               cGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                       <  1:683225/140‑1 (MQ=255)                                                                                                               cGGAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                       >  2:152772/1‑140 (MQ=255)                                                                                                                ggAGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGc                      >  2:495744/1‑140 (MQ=255)                                                                                                                  aGTAACAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                    <  1:871039/140‑1 (MQ=255)                                                                                                                       cAATCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCtt               >  2:104182/1‑140 (MQ=255)                                                                                                                          tCGCGACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAg            <  1:135832/140‑1 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCGATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCGTGTTTTTTTTATGCTTTCCTTAAGaa          >  1:772182/1‑140 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                    >  2:773838/1‑130 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                    <  1:773838/130‑1 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaa          >  1:730356/1‑140 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaa          >  1:64834/1‑140 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaa          <  1:505685/140‑1 (MQ=255)                                                                                                                            gcgACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaa          >  1:333730/1‑140 (MQ=255)                                                                                                                              gACAATACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaaca        >  1:230007/1‑140 (MQ=255)                                                                                                                                 aaTACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACt     <  2:1094091/140‑1 (MQ=255)                                                                                                                                   tACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa   >  2:521637/1‑140 (MQ=255)                                                                                                                                   tACTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa   <  2:198781/140‑1 (MQ=255)                                                                                                                                    aCTCCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAc  <  1:495744/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         GCGGACTCGCCAGTGGTGCCAACAGAAACGATCGCCGAAGTACCGCTGGCGACATGATAATCAATCAGTTTTTTCAAGCTAGCCCGACAGACATTACCTTTTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAC  >  NC_000913/2599615‑2599885

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: ‑