Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,185,729	T→G	9.1%	S310A (TCC→GCC)	purL →	phosphoribosylformylglycinamidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,185,729	0	T	G	9.1%	172.5 / 9.8	66	S310A (TCC→GCC)	purL	phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base T (26/34);  new base G (2/4);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

TGATCGTCGGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTATCGTCGAC  >  NC_002947/1185580‑1185878                                                                                                                                                      |                                                                                                                                                       tGATCGTCGGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCt                                                                                                                                                        <  1:119516/150‑1 (MQ=255)    tcgtcgGCAACGTTGCCGGCCGTTTCTTCACACACCTTCTCTTCCGCCAGGACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCg                                                                                                                                                     <  1:1173159/150‑1 (MQ=255)     cgtcgGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCgg                                                                                                                                                    <  1:1101077/150‑1 (MQ=255)       tcgGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTg                                                                                                                                                  >  2:278778/1‑150 (MQ=255)         ggCAAGGTTGCCGGCCGTTTGTGCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTgga                                                                                                                                                <  2:538077/150‑2 (MQ=255)         ggCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGc                                                                                                                                                <  1:315824/150‑1 (MQ=255)           cAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTggagg                                                                                                                                              <  2:325640/150‑4 (MQ=255)           cAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCATGATGAAGGTGGAAACGCACAACCAACCGCCCGCAATCGCCCCGTGCTCCGGCGCCTCCACCGGCTCCGGTGGCGa                                                                                                                                              >  1:570026/1‑150 (MQ=255)                    cggccgTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCTGGTGGCGAAATCCGcga                                                                                                                                     >  2:789259/1‑150 (MQ=255)                              ttcCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCt                                                                                                                           >  2:574221/1‑150 (MQ=255)                              ttcCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCt                                                                                                                           >  1:362154/1‑150 (MQ=255)                                 ccGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAcc                                                                                                                        >  2:75860/1‑150 (MQ=255)                                 ccGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAcc                                                                                                                        >  2:538777/1‑150 (MQ=255)                                 ccGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAAGTGCTAcc                                                                                                                        >  1:51203/1‑150 (MQ=255)                                  agAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccg                                                                                                                       >  2:1203978/2‑150 (MQ=255)                                    aaCCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggc                                                                                                                     <  2:1148145/150‑1 (MQ=255)                                     aCCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggcc                                                                                                                    <  1:503033/150‑1 (MQ=255)                                      cccTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggccg                                                                                                                   <  2:51203/150‑1 (MQ=255)                                         tGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTgga                                                                                                                                                >  1:161548/1‑117 (MQ=255)                                         tGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTgga                                                                                                                                                <  2:161548/118‑2 (MQ=255)                                                       aCGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAgg                                                                                                  <  1:464037/150‑1 (MQ=255)                                                        cGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAggc                                                                                                 <  2:362154/150‑1 (MQ=255)                                                          gcgcGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAggccg                                                                                               <  1:75860/150‑1 (MQ=255)                                                              ggTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCt                                                                                           <  2:1109354/150‑1 (MQ=255)                                                                  cAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑Acc                                                                                      <  1:278778/150‑1 (MQ=255)                                                                   aGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCg                                                                                     <  2:1127532/150‑1 (MQ=255)                                                                         cGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCa                                                                               <  2:570026/150‑1 (MQ=255)                                                                           gTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCg                                                                                                                                               <  1:590400/85‑1 (MQ=255)                                                                           gTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCg                                                                                                                                               >  2:590400/1‑85 (MQ=255)                                                                           gTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCcgc                                                                                                                                       >  1:587719/1‑93 (MQ=255)                                                                           gTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCcgc                                                                                                                                       <  2:587719/93‑1 (MQ=255)                                                                                  tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCa                                                                      >  2:642011/1‑150 (MQ=255)                                                                                  tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCCGTGGCGAAATCCGCGACCCCCGTCATCCCGGCCGCGGACCCCAGCCCAAGGCTTGCCGG‑ACCGCCGTCACAgtgtaaa                                                                      >  2:317605/1‑147 (MQ=255)                                                                                         gaAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCaa                                                                                                    <  2:358769/114‑1 (MQ=255)                                                                                         gaAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCaa                                                                                                    >  1:358769/1‑114 (MQ=255)                                                                                          aaGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTgc                                                                <  2:406799/148‑1 (MQ=255)                                                                                          aaGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTgc                                                                >  1:406799/1‑148 (MQ=255)                                                                                                 aaaCGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCgg                                                       <  1:600397/150‑1 (MQ=255)                                                                                                  aaCGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGc                                                                                                        >  1:1072724/1‑101 (MQ=255)                                                                                                  aaCGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGc                                                                                                        <  2:1072724/101‑1 (MQ=255)                                                                                                   aCGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGc                                                     <  2:746399/150‑1 (MQ=255)                                                                                                     gCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCtt                                                   >  1:913062/1‑150 (MQ=255)                                                                                                     gCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCtt                                                   <  1:1203978/150‑1 (MQ=255)                                                                                                      cacaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCTGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑CCCGGATTCACCGTGTCCAACCTGCGCATTCCGGGCTTc                                                  <  1:789259/150‑1 (MQ=255)                                                                                                       acaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTgc                                                                <  2:1036010/135‑1 (MQ=255)                                                                                                       acaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTgc                                                                >  1:1036010/1‑135 (MQ=255)                                                                                                               ccGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCt                                         >  1:329170/1‑150 (MQ=255)                                                                                                                 gACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTg                                       <  2:329170/150‑1 (MQ=255)                                                                                                                  aCCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGTCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGTGCCAAGCCCAAGGTCGGCCTG‑ACCGGCTTCACCTTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTgg                                      <  1:370683/150‑1 (MQ=255)                                                                                                                          cGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGc                              <  1:642011/150‑1 (MQ=255)                                                                                                                             cccATTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGTCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTa                           >  2:866858/1‑150 (MQ=255)                                                                                                                                   ctcCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCtt                                        >  1:431339/1‑131 (MQ=255)                                                                                                                                   ctcCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCtt                                        <  2:431339/131‑1 (MQ=255)                                                                                                                                     ccGGCGCCTCCACCGGCTCAGGTGGCAAAATCCGCGACGGAGGTGCTACCGGCCGGGGCCCCACGCCCCAGTCCTGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGc                   <  1:812063/150‑1 (MQ=255)                                                                                                                                         tttcgccACCGGAGCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTgg                                      >  1:825979/6‑127 (MQ=255)                                                                                                                                         tttcgccACCGGAGCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTgg                                      <  2:825979/122‑1 (MQ=255)                                                                                                                                          gcCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAg              <  1:538777/150‑1 (MQ=255)                                                                                                                                             tCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAg                                                                                                         <  1:710145/57‑1 (MQ=255)                                                                                                                                             tCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAg                                                                                                         >  2:710145/1‑57 (MQ=255)                                                                                                                                             tCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTAt                          >  2:376625/1‑135 (MQ=255)                                                                                                                                             tCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTAt                          <  1:376625/135‑1 (MQ=255)                                                                                                                                              ccACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCAcc                                                                             >  2:998917/1‑83 (MQ=255)                                                                                                                                              ccACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCAcc                                                                             <  1:998917/83‑1 (MQ=255)                                                                                                                                               cACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAg                                <  2:583151/127‑1 (MQ=255)                                                                                                                                               cACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAg                                >  1:583151/1‑127 (MQ=255)                                                                                                                                               cACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTAt         >  1:45721/1‑150 (MQ=255)                                                                                                                                                aCCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGACAAG‑CCATACCTTAACTGTCCCGTCTGCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTAtc        <  1:574221/150‑1 (MQ=255)                                                                                                                                                      tCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTATCGTCGAc  >  2:387950/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       TGATCGTCGGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTG‑ACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTATCGTCGAC  >  NC_002947/1185580‑1185878

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