Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,662,017	G→A	6.4%	A49T (GCA→ACA)	purT →	phosphoribosylglycinamide formyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,662,017	0	G	A	6.4%	186.1 / 2.2	63	A49T (GCA→ACA)	purT	phosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base A (2/2);  total (30/33)
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.65e-01

GAAATGACCCGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCTGGAAAAC  >  NC_002947/1661870‑1662160                                                                                                                                                    |                                                                                                                                                gAAATGACCCGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCgg                                                                                                                                                 >  6:78229/1‑148 (MQ=255)         ccGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAgg                                                                                                                                        <  1:282863/149‑1 (MQ=255)         ccGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAgg                                                                                                                                        <  5:318044/149‑1 (MQ=255)         ccGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGAAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGATCCGGCAATGCAgg                                                                                                                                        >  3:782/1‑149 (MQ=255)          cGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGt                                                                                                                                       <  3:235311/149‑1 (MQ=255)           gTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGTGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTg                                                                                                                                      <  7:79495/149‑1 (MQ=255)            tATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTgg                                                                                                                                     <  8:188184/149‑1 (MQ=255)              tCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTgg                                                                                                                                     >  6:99787/1‑147 (MQ=255)              tCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTgg                                                                                                                                     <  5:99787/147‑1 (MQ=255)                 gAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCAc                                                                                                                                <  4:170894/149‑1 (MQ=255)                        aTTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGcc                                                                                                                         <  2:148135/149‑1 (MQ=255)                            tCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCAtgt                                                                                                                     <  4:782/149‑1 (MQ=255)                              gCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTgg                                                                                                                   >  2:40577/1‑149 (MQ=255)                               cccACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGt                                                                                                                  <  5:78229/149‑1 (MQ=255)                               cccACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGt                                                                                                                  >  5:240792/1‑149 (MQ=255)                                  aCCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCaa                                                                                                                                              <  6:160559/118‑1 (MQ=255)                                  aCCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCaa                                                                                                                                              <  6:264103/118‑1 (MQ=255)                                  aCCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCaa                                                                                                                                              >  5:264103/1‑118 (MQ=255)                                  aCCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCaa                                                                                                                                              >  5:160559/1‑118 (MQ=255)                                          ccGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTc                                                                                                        >  3:18239/1‑148 (MQ=255)                                            gtgtACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGAt                                                                                                     >  1:91616/1‑149 (MQ=255)                                            gtgtACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGAt                                                                                                     <  1:40577/149‑1 (MQ=255)                                              gtACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATgg                                                                                                   >  1:259680/1‑149 (MQ=255)                                                aCTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATg                                                                                                    >  8:23608/1‑146 (MQ=255)                                                aCTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATg                                                                                                    <  7:23608/146‑1 (MQ=255)                                                  tGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGtt                                                                                               >  8:294074/1‑149 (MQ=255)                                                  tGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGtt                                                                                               >  7:130453/1‑149 (MQ=255)                                                         cGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGTCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCt                                                                                          <  5:175612/147‑1 (MQ=255)                                                         cGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCAt                                                                                                                       >  8:195476/1‑118 (MQ=255)                                                         cGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCAt                                                                                                                       <  7:195476/118‑1 (MQ=255)                                                            cTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGACAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTgcgcg                                                                                     >  5:1654/1‑149 (MQ=255)                                                               tgGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGcgcgcg                                                                                   <  2:91616/148‑1 (MQ=255)                                                                         ggCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAg                                                                         <  6:299894/148‑1 (MQ=255)                                                                         ggCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGc                                                                        <  3:266642/149‑1 (MQ=255)                                                                          gCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGcc                                                                       <  4:18239/149‑1 (MQ=255)                                                                              ggAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCgaga                                                                   <  2:59227/149‑1 (MQ=255)                                                                               gAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGAACGATGGGGTTGCCATGCGCGCGGTGAACGAAGCCGAGaa                                                                  >  8:49778/1‑149 (MQ=255)                                                                                              cTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGAGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGt                                                   >  8:267381/1‑149 (MQ=255)                                                                                                   gcgcCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCg                                               >  8:140082/1‑148 (MQ=255)                                                                                                   gcgcCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGACAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCg                                               >  6:247017/1‑148 (MQ=255)                                                                                                       cTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAAt                                          <  2:259680/149‑1 (MQ=255)                                                                                                           gCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGACAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGaa                                      <  6:1654/149‑1 (MQ=255)                                                                                                            cGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAg                                     <  5:56363/149‑1 (MQ=255)                                                                                                            cGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGGTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAg                                     <  4:214208/149‑1 (MQ=255)                                                                                                             gTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCAGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGc                                    >  8:155622/1‑149 (MQ=255)                                                                                                                gAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCgg                                              >  6:206966/1‑136 (MQ=255)                                                                                                                gAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCgg                                              <  5:206966/136‑1 (MQ=255)                                                                                                                  aGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATcg                               >  2:182944/1‑149 (MQ=255)                                                                                                                   gTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCAc                                                           <  1:79572/120‑1 (MQ=255)                                                                                                                   gTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCAc                                                           >  2:79572/1‑120 (MQ=255)                                                                                                                     aaTCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATcgcca                            >  1:142095/1‑149 (MQ=255)                                                                                                                                aCCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCAc                                                           <  6:249472/107‑1 (MQ=255)                                                                                                                                aCCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCAc                                                           >  5:249472/1‑107 (MQ=255)                                                                                                                                  cGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGaa                                      <  6:7950/126‑1 (MQ=255)                                                                                                                                  cGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGaa                                      >  5:7950/1‑126 (MQ=255)                                                                                                                                       cGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGGGTTGCCCTGCGCGCGGTGATCGAACCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGc          <  7:267381/149‑1 (MQ=255)                                                                                                                                        gCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCTCCACCGCCACCCTGGTCGAGCt         <  7:49778/149‑1 (MQ=255)                                                                                                                                        gCCAATGCTCCGACAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCt         <  5:247017/149‑1 (MQ=255)                                                                                                                                         ccAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCTg        >  4:233723/1‑149 (MQ=255)                                                                                                                                         ccAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCTg        >  3:7315/1‑149 (MQ=255)                                                                                                                                          cAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCTg        >  1:51116/1‑148 (MQ=255)                                                                                                                                           aaTGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCaa                                                                                                               <  5:20955/44‑1 (MQ=255)                                                                                                                                           aaTGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCaa                                                                                                               >  6:20955/1‑44 (MQ=255)                                                                                                                                                tCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCTGGAAAAc  >  5:33078/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                GAAATGACCCGTATCGGAACCCCATTGTCGCCCACCGCGACCCGTGTACTGCTCTGCGGCTGTGGCGAGCTGGGCAAGGAAGTGGTAATCGAGCTGCAGCGCCTGGGCGTCGAAGTAATCGCCGTCGACCGCTACGCCAATGCTCCGGCAATGCAGGTGGCGCACCGCAGCCATGTGGTCAACATGCTCGATGGCGTTGCCCTGCGCGCGGTGATCGAAGCCGAGAAGCCGCACTACATCGTGCCGGAAATCGAAGCCATCGCCACCGCCACCCTGGTCGAGCTGGAAAAC  >  NC_002947/1661870‑1662160

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