Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,234,944	G→C	6.0%	R41G (CGC→GGC)	PP_3711 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,234,944	0	G	C	6.0%	210.1 / 2.2	66	R41G (CGC→GGC)	PP_3711	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (26/36);  new base C (2/2);  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 = 9.69e-01

AGGCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCAT  >  NC_002947/4234799‑4235092                                                                                                                                                  |                                                                                                                                                     aGGCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTg                                                                                                                                                   <  2:1787/149‑1 (MQ=255)   gCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTggg                                                                                                                                                 >  3:183784/1‑149 (MQ=255)         gcCCGTTTCTTCACAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCGCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaa                                                                                                                                           <  4:83710/149‑2 (MQ=255)           ccGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAggg                                                                                                                                         >  8:153074/1‑149 (MQ=255)              tttCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTc                                                                                                                                      <  7:21829/149‑1 (MQ=255)                 cttcAACGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCa                                                                                                                                    <  1:325609/148‑1 (MQ=255)                   tcAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGtt                                                                                                                                 >  7:104684/1‑149 (MQ=255)                     aaaGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTgg                                                                                                                               >  5:260321/1‑149 (MQ=255)                     aaaGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTgg                                                                                                                               <  4:183784/149‑1 (MQ=255)                         cGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGcgac                                                                                                                           >  2:194880/1‑149 (MQ=255)                             cGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGAt                                                                                                                        <  5:354281/148‑1 (MQ=255)                             cGGTGCGGCTCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATa                                                                                                                       <  7:240243/149‑1 (MQ=255)                               gTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATAcc                                                                                                                     <  2:14328/149‑1 (MQ=255)                                  ctgcGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaag                                                                                                                                          <  8:334196/125‑3 (MQ=255)                                  ctgcGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCCCAGCGCGCCaag                                                                                                                                          >  7:334196/1‑123 (MQ=255)                                     cgcgcCTGCTCCAGGTTGAATGCCGGCAACGCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAg                                                                                                                <  7:340701/148‑1 (MQ=255)                                      gcgcCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAggg                                                                                                              <  6:260321/149‑1 (MQ=255)                                      gcgcCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAggg                                                                                                              <  2:138990/149‑1 (MQ=255)                                          cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGTGTTCCAGTTGGCGACGATACCTTCAGGGGTGa                                                                                                          <  1:27469/149‑1 (MQ=255)                                          cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGt                                                                                                            >  4:349750/1‑147 (MQ=255)                                          cTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGa                                                                                                          <  1:21668/149‑1 (MQ=255)                                             cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAg                                                                                                         <  3:28669/147‑1 (MQ=255)                                             cTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGGGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCa                                                                                                       <  6:242162/149‑1 (MQ=255)                                                    ttGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGAt                                                                                                <  6:199720/149‑1 (MQ=255)                                                            cGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAAt                                                                                        >  2:142364/1‑149 (MQ=255)                                                            cGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAAt                                                                                        >  4:18236/1‑149 (MQ=255)                                                                aaCCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccac                                                                                    >  1:107636/1‑149 (MQ=255)                                                                  cccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccacca                                                                                  >  3:193979/1‑149 (MQ=255)                                                                  cccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATccacca                                                                                  >  5:62556/1‑149 (MQ=255)                                                                   ccGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCcaccac                                                                                 >  4:182963/1‑149 (MQ=255)                                                                         gcgcgcTCGGCATCGGAGTAGAACAACGAGTAGTGCTGGCCGACGATATCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCACGGGTCAGCACGTAGATGGAATAATCCAACACGctctc                                                                           >  2:185777/1‑148 (MQ=255)                                                                          cgcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGTGGTGAGCAGGTAGATGGCATAATCCACCACGctct                                                                            <  8:104684/147‑1 (MQ=255)                                                                          cgcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATTATCCACCACGctct                                                                            <  6:8928/147‑1 (MQ=255)                                                                           gcgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGa                                                                         <  1:142364/149‑1 (MQ=255)                                                                            cgcTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGAt                                                                        <  6:62556/149‑1 (MQ=255)                                                                              cTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATACTCCACCACGCTCTGGATc                                                                       >  7:321343/1‑148 (MQ=255)                                                                              cTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCa                                                                      <  7:91710/149‑1 (MQ=255)                                                                                     tCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCg                                                               <  3:182963/149‑1 (MQ=255)                                                                                       ggTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCgg                                                              <  3:303219/148‑1 (MQ=255)                                                                                         tgtAGAACAACGAGTAGTGCTGGCCGACGATTTCATCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGTGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGAAGGTAGATGGCATAATCCACAACGCTCTGGATCAACAGACGGTa                                                            >  6:241888/1‑148 (MQ=255)                                                                                          gtAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACa                                                          <  3:349750/149‑1 (MQ=255)                                                                                           tAGAACTACGAGTAGTGCTGGCCGTCGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGTTACa                                                          <  5:47486/148‑1 (MQ=255)                                                                                             gAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACatga                                                       <  7:75855/149‑1 (MQ=255)                                                                                              aacaacGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACatgat                                                      >  1:269505/1‑149 (MQ=255)                                                                                                 aacGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTc                                                   >  1:359862/1‑149 (MQ=255)                                                                                                    gAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCggg                                                <  3:152852/149‑1 (MQ=255)                                                                                                    gAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATAAACAGCCGGTACATGATGTCggg                                                >  7:161681/1‑149 (MQ=255)                                                                                                      gtagtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGctct                                                                            <  4:346706/119‑1 (MQ=255)                                                                                                      gtagtGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGctct                                                                            >  3:346706/1‑119 (MQ=255)                                                                                                            cTGGCGGACGATTTCCTCGGCAGGGTACCCCTTGGCGCGCTGGGCGCGAGGGTTCCAGTTGGCGACGATACCTTCTGGGGTGAGCAGGTAGTTGGCATAATCCACCACGCTGTGGATCATCAGCCGGTACATGATGGCGGGGTTtgct                                         <  1:185777/148‑1 (MQ=255)                                                                                                             tGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCgg                                                              >  2:42234/1‑126 (MQ=255)                                                                                                             tGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCgg                                                              <  1:42234/126‑1 (MQ=255)                                                                                                              tgCCGAGGAAATCGTCGGCCCGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                                       >  6:277347/2‑148 (MQ=255)                                                                                                               gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTtgctgc                                       <  2:107636/147‑1 (MQ=255)                                                                                                               gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcga                                     <  4:193979/149‑1 (MQ=255)                                                                                                               gCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCGGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcga                                     <  4:108734/149‑1 (MQ=255)                                                                                                                 cgacgaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTgcgagc                                   >  5:9892/1‑149 (MQ=255)                                                                                                                     gaTTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACa                               >  3:3312/1‑149 (MQ=255)                                                                                                                        ttCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTACAGTTGGCGACGATACCTTGAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCAGGATCAACAGCCGGTAAATGATGTCGGGGTTTGCTGGGAGCGACATg                             >  5:258151/1‑148 (MQ=255)                                                                                                                           ctcggctcggTACCCCTTGGCGCGTTGGGCGCCAGTGTTCCAGTTGGCGACGATTCCTTCAGGGGTGAGCAGGTTGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAAc                          <  5:241888/148‑1 (MQ=255)                                                                                                                              ggCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGc                                                                 <  4:288856/106‑1 (MQ=255)                                                                                                                              ggCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGc                                                                 >  3:288856/1‑106 (MQ=255)                                                                                                                                         ccctGGCGCCCAGCGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGCTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc           >  5:367196/4‑149 (MQ=255)                                                                                                                                         ccTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc           <  6:9892/149‑1 (MQ=255)                                                                                                                                         ccTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc           <  2:359862/149‑1 (MQ=255)                                                                                                                                          cTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTc           >  3:220111/1‑148 (MQ=255)                                                                                                                                              tcgtgcTGGGCGCCATGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCATGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTAGATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGTAgc      <  6:186634/145‑1 (MQ=255)                                                                                                                                                gcgcTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAgcg     >  2:366299/1‑148 (MQ=255)                                                                                                                                                  gcTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCAt  >  1:99784/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     AGGCGTGCGCCCGTTTCTTCAAAGCGTCCGGTGCTGCGCGCCTGCTCCAGGTTGAATGCCGGCAACCCGGCTGCGCGCTCGGCATCGGTGTAGAACAACGAGTAGTGCTGGCCGACGATTTCCTCGGCACGGTACCCCTTGGCGCGCTGGGCGCCAGGGTTCCAGTTGGCGACGATACCTTCAGGGGTGAGCAGGTAGATGGCATAATCCACCACGCTCTGGATCAACAGCCGGTACATGATGTCGGGGTTTGCTGCGAGCGACATGAACTTTACCTGAACCTTCGAAGCGCAT  >  NC_002947/4234799‑4235092

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