Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	453,467	G→A	5.0%	P27S (CCG→TCG)	ampG ←	muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	453,467	0	G	A	5.0%	252.3 / 2.5	79	P27S (CCG→TCG)	ampG	muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base A (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTTGATGAATT  >  NC_000913/453321‑453614                                                                                                                                                   |                                                                                                                                                    tGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTa                                                                                                                                                   <  1:299567/149‑1 (MQ=255)         ttttgCGGTCGATCAGCGGTGACCTGAGGAATTTAATAACGTAAGCCTGGGCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGTTACCGGATGTCAGGGCGAGCGGTAGCCCGGaa                                                                                                                                           <  4:138992/145‑1 (MQ=255)           tAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAg                                                                                                                                          >  7:447280/1‑148 (MQ=255)            aGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGc                                                                                                                                         >  1:61545/1‑148 (MQ=255)            aGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGc                                                                                                                                         >  8:164113/1‑148 (MQ=255)                  ccATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAAcc                                                                                                                                   >  6:168904/1‑148 (MQ=255)                   cATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAAccc                                                                                                                                  <  1:20394/148‑1 (MQ=255)                   cATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAAcc                                                                                                                                   >  5:175074/1‑147 (MQ=255)                   cATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCa                                                                                                                                 <  6:175074/149‑1 (MQ=255)                   cATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCa                                                                                                                                 <  3:374233/149‑1 (MQ=255)                     tCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAg                                                                                                                                >  6:10609/1‑148 (MQ=255)                          ggTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAg                                                                                                                          <  5:344274/149‑1 (MQ=255)                          ggTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAg                                                                                                                          <  8:84418/149‑1 (MQ=255)                           gTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGc                                                                                                                         <  5:148979/149‑1 (MQ=255)                                     gAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATc                                                                                                                            >  5:48878/1‑136 (MQ=255)                                     gAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATc                                                                                                                            <  6:48878/136‑1 (MQ=255)                                         ttAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCa                                                                                                                                 >  6:184985/1‑127 (MQ=255)                                         ttAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCa                                                                                                                                 <  5:184985/127‑1 (MQ=255)                                            aaaaCGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCCATATTGCTGatcctg                                                                                                         <  2:43017/148‑6 (MQ=255)                                            aaaaCGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAAcgcg                                                                                                         <  5:10609/148‑1 (MQ=255)                                             aaaCGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGc                                                                                                       >  1:297255/1‑149 (MQ=255)                                                 gTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGtt                                                                                                   <  6:152373/149‑1 (MQ=255)                                                 gTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGtt                                                                                                   >  7:158838/1‑149 (MQ=255)                                                          gCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAAt                                                                                            <  8:224234/147‑1 (MQ=255)                                                            cTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAACAATATGTCTGAACGCGGCTGTTGAAAAATACg                                                                                         >  8:289969/1‑148 (MQ=255)                                                            cTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGGCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGt                                                                                        <  4:351348/149‑1 (MQ=255)                                                               ccAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTaaa                                                                                     >  8:282815/1‑149 (MQ=255)                                                               ccAGAGAGAAGAAAACAATGGTTTTGAGATCGATATTCTAGACCGTCATCCAGGCCTGTAATGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTa                                                                                       <  7:349264/147‑1 (MQ=255)                                                                   agagaAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATAtt                                                                                 <  3:153237/149‑1 (MQ=255)                                                                       aagaaACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACt                                                                             <  7:164113/149‑1 (MQ=255)                                                                       aagaaACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACt                                                                             >  3:92749/1‑149 (MQ=255)                                                                        agaaACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCACGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTg                                                                            <  3:274932/149‑1 (MQ=255)                                                                            aCCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACa                                                                        >  3:344689/1‑149 (MQ=255)                                                                             ccAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACAt                                                                       <  8:447280/149‑1 (MQ=255)                                                                              cAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGAc                                                                         <  4:144308/146‑1 (MQ=255)                                                                              cAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGAc                                                                         >  3:144308/1‑146 (MQ=255)                                                                              cAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATg                                                                      >  6:221437/1‑149 (MQ=255)                                                                               aaTGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATgg                                                                     >  7:119521/1‑149 (MQ=255)                                                                               aaTGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATgg                                                                     <  8:158838/149‑1 (MQ=255)                                                                               aaTGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCGGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATgg                                                                     <  2:22839/149‑1 (MQ=255)                                                                                aTGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCACAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATggg                                                                    >  7:211572/1‑149 (MQ=255)                                                                                        gagaTCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTAc                                                             <  2:61545/148‑1 (MQ=255)                                                                                           aTCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGc                                                         >  1:210682/1‑149 (MQ=255)                                                                                            tCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGcc                                                        >  7:89131/1‑149 (MQ=255)                                                                                                 atTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGc                                                                   >  5:353864/1‑133 (MQ=255)                                                                                                 atTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGc                                                                   <  6:353864/133‑1 (MQ=255)                                                                                                  ttCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCATGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATTGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCtt                                                  >  7:323614/1‑149 (MQ=255)                                                                                                    ctcGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTg                                                 <  7:282815/148‑1 (MQ=255)                                                                                                         ccGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCgg                                           >  2:222398/1‑149 (MQ=255)                                                                                                          cGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCggg                                          >  3:373716/1‑149 (MQ=255)                                                                                                                    gCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCAGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGccc                                                       >  1:397609/1‑126 (MQ=255)                                                                                                                    gCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCAGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGccc                                                       <  2:397609/126‑1 (MQ=255)                                                                                                                      cTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAg                               >  7:114312/1‑148 (MQ=255)                                                                                                                       tGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGcctgta                                    <  3:241720/142‑5 (MQ=255)                                                                                                                       tGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGcctgta                                    >  4:241720/1‑138 (MQ=255)                                                                                                                       tGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGa                              >  6:179285/1‑148 (MQ=255)                                                                                                                       tGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGa                              >  5:342831/1‑148 (MQ=255)                                                                                                                       tGTAAGGTACCGGATGTCAGGGCGAGCGATAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGa                              >  7:248800/1‑148 (MQ=255)                                                                                                                       tGTAAGGTACCGGATGTCAGGGCGAGCGATAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGa                              >  7:255482/1‑148 (MQ=255)                                                                                                                        gTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGaaa                            >  7:67797/1‑149 (MQ=255)                                                                                                                         tAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAAt                           >  7:448055/1‑149 (MQ=255)                                                                                                                             ggcctggATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAAc                       <  5:21149/144‑1 (MQ=255)                                                                                                                             gTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACt                                                                             <  1:224985/95‑1 (MQ=255)                                                                                                                             gTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACt                                                                             >  2:224985/1‑95 (MQ=255)                                                                                                                              tACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAAcg                      >  2:36693/1‑149 (MQ=255)                                                                                                                                ccGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAAcgc                     <  5:179285/148‑1 (MQ=255)                                                                                                                                 cGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTAGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTCCACCGGGCCAGAGGGCAGAAATTAAcgcgc                   >  4:70816/1‑149 (MQ=255)                                                                                                                                  ggATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCa                                >  4:395205/1‑135 (MQ=255)                                                                                                                                  ggATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCa                                <  3:395205/135‑1 (MQ=255)                                                                                                                                     tGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTc               <  5:168904/149‑1 (MQ=255)                                                                                                                                       tCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTg             >  6:413230/1‑149 (MQ=255)                                                                                                                                         aGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTaa                        >  4:89998/1‑136 (MQ=255)                                                                                                                                         aGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTaa                        <  3:89998/136‑1 (MQ=255)                                                                                                                                         aGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCt           >  3:206355/1‑149 (MQ=255)                                                                                                                                            gCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTtga        >  1:267741/1‑149 (MQ=255)                                                                                                                                             cGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACTTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTtgat       <  4:344689/149‑1 (MQ=255)                                                                                                                                              gAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTtgat       >  7:26519/1‑148 (MQ=255)                                                                                                                                               aGCGATAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTtgatga     <  8:248800/149‑1 (MQ=255)                                                                                                                                               aGCGATAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTtgatga     <  8:255482/149‑1 (MQ=255)                                                                                                                                                  ggTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTTGATGAAtt  >  5:233241/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TGGAGGCGTGTAGCGGTCCATCAGCGGTGACCAGAGGAATTTAAAAACGTAAGCCTGGCCTACCAGAGAGAAGAAACCAATGGTTTTGAGATCGATATTCTCGACCGTCATCCAGGCCTGTAAGGTACCGGATGTCAGGGCGAGCGGTAGCCCGGAAGCAAAACCCAGGATCAGCAATATGGCTGAACGCGGCTGTTGAAAAATACGTAAATATTGACTGGACATGGGCGTCTACAGGCCCGGCTTGCACCGGGCCAGAGGGCAGAAATTAACGCGCGTTCTGCTTGATGAATT  >  NC_000913/453321‑453614

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