Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,844,871	1	.	G	38.2%	95.2 / 136.8	76	intergenic (+48/‑249)	ccoP‑II/PP_4259	cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (25/22);  new base G (15/14);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAT  >  NC_002947/4844723‑4845020                                                                                                                                                      |                                                                                                                                                      tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                      <  1:1213690/150‑1 (MQ=255)  gCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCT‑Tg                                                                                                                                                     <  2:330880/150‑1 (MQ=255)  gCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑Tg                                                                                                                                                     <  1:696577/150‑1 (MQ=255)  gCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGCCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCtgtg                                                                                                                                                     <  1:1529962/150‑4 (MQ=37)    cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCCCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGTCTGGTCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTa                                                                                                                                                   <  1:413220/150‑1 (MQ=255)     ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAg                                                                                                                                                  <  2:1382141/150‑1 (MQ=255)       cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGa                                                                                                                                                <  2:1345838/150‑1 (MQ=255)           ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                            <  1:2005457/150‑1 (MQ=255)           ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                            <  2:1148907/150‑1 (MQ=255)           ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            <  2:1561851/150‑1 (MQ=255)            gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGGCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑GGTAGGAGCggg                                                                                                                                           <  1:2134379/150‑1 (MQ=255)            gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg                                                                                                                                           <  2:134968/150‑1 (MQ=255)                   tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                      <  1:811596/132‑1 (MQ=255)                   tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                      >  2:811596/1‑132 (MQ=255)                   tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAccc                                                                                                                                    <  1:1744652/150‑1 (MQ=255)                       gggCAACGACAAGGTGCACTTGCGGGCGGCTTACGTCTGCAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTCAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTAAACCcgctc                                                                                                                                <  1:27747/150‑3 (MQ=38)                        ggCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                      <  2:489731/127‑1 (MQ=255)                        ggCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                                      >  1:489731/1‑127 (MQ=255)                              gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAg                                                                                                                                                  >  2:1923453/1‑125 (MQ=255)                              gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAg                                                                                                                                                  <  1:1923453/125‑1 (MQ=255)                                cAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            >  1:842606/1‑129 (MQ=255)                                cAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            <  2:842606/129‑1 (MQ=255)                                    gTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            <  1:959272/125‑1 (MQ=255)                                    gTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            >  2:959272/1‑125 (MQ=255)                                          ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAg                                                                                                                                               >  1:1167640/1‑116 (MQ=255)                                          ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAg                                                                                                                                               <  2:1167640/116‑1 (MQ=255)                                                gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg                                                                                                                                           >  1:2241296/1‑114 (MQ=255)                                                gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCggg                                                                                                                                           <  2:2241296/114‑1 (MQ=255)                                                   gcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            <  2:2872506/110‑1 (MQ=25)                                                   gcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            >  1:2872506/1‑110 (MQ=25)                                                    cTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCg                                                                                                   <  1:3074939/150‑1 (MQ=255)                                                                         aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            >  2:2122245/1‑88 (MQ=21)                                                                         aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                            <  1:2122245/88‑1 (MQ=21)                                                                                tGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGc                                                                                                                                              <  1:1870752/79‑1 (MQ=255)                                                                                tGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGc                                                                                                                                              >  2:1870752/1‑79 (MQ=255)                                                                                               tGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                             >  1:2100663/1‑81 (MQ=255)                                                                                               tGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                             <  2:2100663/81‑1 (MQ=255)                                                                                                   tAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCg                                                                                                                                             >  2:1705042/1‑61 (MQ=11)                                                                                                   tAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCg                                                                                                                                             <  1:1705042/61‑1 (MQ=11)                                                                                                          ggcctggcctCTTCGCGGGTAAACCCGCCCCTACACAGGATCCT‑TGTAg                                                                                                                                                  <  1:1292794/49‑1 (MQ=255)                                                                                                          ggcctggcctCTTCGCGGGTAAACCCGCCCCTACACAGGATCCT‑TGTAg                                                                                                                                                  >  2:1292794/1‑49 (MQ=255)                                                                                                                            gTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTg                                                                                                                  <  1:3009118/63‑1 (MQ=255)                                                                                                                            gTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTg                                                                                                                  >  2:3009118/1‑63 (MQ=255)                                                                                                                            gTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg                           >  2:3064350/1‑150 (MQ=255)                                                                                                                                aCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGa                                                                                  >  2:2004306/1‑91 (MQ=255)                                                                                                                                aCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGa                                                                                  <  1:2004306/91‑1 (MQ=255)                                                                                                                                 cccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACGGCATGTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCAAATTCGAGCAAAGGTCGcc                      >  1:1838594/1‑149 (MQ=38)                                                                                                                                 cccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa                      >  2:1510923/1‑150 (MQ=255)                                                                                                                                 cccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAAAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACAACCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa                      >  1:981266/1‑150 (MQ=255)                                                                                                                                  ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  1:382733/1‑150 (MQ=255)                                                                                                                                  ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  1:1972399/1‑150 (MQ=255)                                                                                                                                  ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  1:827908/95‑1 (MQ=21)                                                                                                                                  ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            >  2:827908/1‑95 (MQ=21)                                                                                                                                  ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  2:1747114/1‑150 (MQ=255)                                                                                                                                   cGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTtct                                                                                                     <  1:257713/69‑1 (MQ=255)                                                                                                                                   cGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTtct                                                                                                     >  2:257713/1‑69 (MQ=255)                                                                                                                                    gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCt                                                                                     <  2:1242498/84‑1 (MQ=21)                                                                                                                                    gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCt                                                                                     >  1:1242498/1‑84 (MQ=21)                                                                                                                                    gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                                             <  2:2585932/92‑1 (MQ=21)                                                                                                                                    gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                                             >  1:2585932/1‑92 (MQ=21)                                                                                                                                    gCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  1:1220310/1‑150 (MQ=255)                                                                                                                                      tCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg                 >  1:1620126/1‑150 (MQ=255)                                                                                                                                      tCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg                 >  1:1473901/1‑150 (MQ=255)                                                                                                                                        cTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa               >  2:789712/1‑150 (MQ=255)                                                                                                                                         ctacaaGGATCCTGTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGCCACCCGCCATAACACGTTTAATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa               >  2:138943/6‑150 (MQ=38)                                                                                                                                          acacAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGAAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGAGTAAAGGTCGAAGTGCGAAcc             >  1:1406871/1‑150 (MQ=255)                                                                                                                                           cacaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACACGACTAAAGGTCGCAGTGCGAAccc            >  1:2731270/1‑150 (MQ=255)                                                                                                                                           acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            >  1:1138186/4‑87 (MQ=25)                                                                                                                                           acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  2:1138186/84‑1 (MQ=25)                                                                                                                                            caaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            >  1:2911067/3‑86 (MQ=255)                                                                                                                                            caaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  2:2911067/84‑1 (MQ=255)                                                                                                                                            acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACATCCGCCATAACACGTTTCATAGCAAGAACCCAACACGACTAAAGGTCGCAGTGCAAAcaca           >  2:1191359/1‑150 (MQ=39)                                                                                                                                              aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATa         >  1:2252565/1‑150 (MQ=255)                                                                                                                                              aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCAGGTACTGGCACAGCATGTCTCGCAGGCAAGTGACTTGGATCCATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATa         >  1:2525409/1‑150 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCAc                                               >  2:1538589/1‑110 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCAc                                               <  1:1538589/110‑1 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACt       >  2:116872/1‑150 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACt       <  2:382733/150‑1 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACt       <  1:116872/150‑1 (MQ=255)                                                                                                                                                    cT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                            >  2:1194641/1‑109 (MQ=255)                                                                                                                                                    cT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                            <  1:1194641/109‑1 (MQ=255)                                                                                                                                                    cT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCa   >  1:2834017/1‑150 (MQ=255)                                                                                                                                                     t‑tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAt  >  1:1834235/1‑150 (MQ=255)                                                                                                                                                     t‑tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAt  >  1:585355/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAT  >  NC_002947/4844723‑4845020

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