Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,527	A→C	9.4%	intergenic (‑119/+37)	PP_1662 ← / ← PP_1663	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,859,527	0	A	C	9.4%	183.2 / 10.1	64	intergenic (‑119/+37)	PP_1662/PP_1663	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/32);  new base C (2/4);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

AAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCC  >  NC_002947/1859383‑1859669                                                                                                                                                   |                                                                                                                                               aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTgg                                                                                                                                            >  8:77117/1‑149 (MQ=255)   gATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGc                                                                                                                                          <  7:356522/149‑1 (MQ=255)       gCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCCCAGTGCgggg                                                                                                                                      <  3:176950/149‑1 (MQ=255)          gCGACCGTGTGTCGCACGTCGTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGACGGGCTGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCtt                                                                                                                                   <  6:52234/149‑1 (MQ=255)          gCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCCCAGTGCGGGGCtt                                                                                                                                   <  7:49531/149‑1 (MQ=255)          gCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCAT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCt                                                                                                                                    >  7:382950/1‑148 (MQ=255)            gACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCa                                                                                                                                 >  5:381863/1‑149 (MQ=255)             aCAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAg                                                                                                                                <  7:305599/149‑1 (MQ=255)              cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCCCAGTGCGGGGCttttct                                                                                                                               <  7:188070/149‑5 (MQ=255)              cAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAgg                                                                                                                               >  6:54038/1‑149 (MQ=255)               aGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAgtt                                                                                                                              <  6:381863/149‑3 (MQ=255)                 tgtgtCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTAt                                                                                                                            <  5:54038/149‑1 (MQ=255)                 tgtgtCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTAt                                                                                                                            <  6:83660/149‑1 (MQ=255)                   tgtCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTg                                                                                                                          <  2:112746/149‑1 (MQ=255)                       gCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTTGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTtcagg                                                                                                                      <  5:285430/149‑1 (MQ=255)                              ttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgcg                                                                                                               >  4:350116/1‑149 (MQ=255)                                cttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                                                                                                      <  8:168590/140‑1 (MQ=255)                                cttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTgaagg                                                                                                                      >  7:168590/1‑140 (MQ=255)                                 tttATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGTTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGa                                                                                                            <  1:126635/149‑1 (MQ=255)                                        agTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTg                                                                                                     <  7:77117/149‑1 (MQ=255)                                         gTTTGTTATCAGACTAAGTAAGTAAAACAGGCATGAAACAAGGTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                                                                                    <  6:330668/149‑1 (MQ=255)                                          tttGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGcc                                                                                                   >  8:5282/1‑149 (MQ=255)                                             gTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCtt                                                                                                >  6:165395/1‑149 (MQ=255)                                                aGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCGTCAc                                                                                             >  4:217396/1‑149 (MQ=255)                                                 gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCAcc                                                                                            >  8:217370/1‑149 (MQ=255)                                                 gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCAc                                                                                             <  6:244795/148‑1 (MQ=255)                                                 gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCAc                                                                                             >  5:244795/1‑148 (MQ=255)                                                  cAGACTAAGTAAGTAAAACAGGCATGAAACAATTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCg                                                                                           <  2:47127/149‑1 (MQ=255)                                                    gACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                                                                                         >  4:66667/1‑149 (MQ=255)                                                     aCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCt                                                                                        >  4:99132/1‑149 (MQ=255)                                                      cTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCtt                                                                                       <  7:88853/149‑1 (MQ=255)                                                         agtaagtaAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGcc                                                                                    <  3:99132/149‑1 (MQ=255)                                                             agtaaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATc                                                                                 <  7:5282/148‑1 (MQ=255)                                                                 aaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCg                                                                            >  5:133088/1‑149 (MQ=255)                                                                   aCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg                                                                          <  6:133088/149‑1 (MQ=255)                                                                     aGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTg                                                                                                                          <  5:374292/99‑1 (MQ=255)                                                                     aGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTg                                                                                                                          >  6:374292/1‑99 (MQ=255)                                                                     aGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGc                                                                         >  6:149609/1‑148 (MQ=255)                                                                       gCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCAt                                                                      <  5:165395/149‑1 (MQ=255)                                                                        cATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATc                                                                     <  7:62671/149‑1 (MQ=255)                                                                             aaCAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCTTTCACCGGCTTGCCATCGACCGTGCCATCCTgca                                                                <  3:350116/149‑1 (MQ=255)                                                                              aCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTg                                                                  >  1:60348/1‑147 (MQ=255)                                                                              aCAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcag                                                               <  8:382950/149‑1 (MQ=255)                                                                               cAAGTTTCATGAAGCTCTAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcagc                                                              <  5:149609/149‑1 (MQ=255)                                                                               cAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcag                                                               >  5:234111/1‑148 (MQ=255)                                                                                  gTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGTCGGGTCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACg                                                           <  7:234923/149‑1 (MQ=255)                                                                                   tttCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGa                                                          >  1:127833/1‑149 (MQ=255)                                                                                       aTGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATca                                                       <  6:234111/148‑1 (MQ=255)                                                                                         gAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATcaca                                                     >  7:268254/1‑148 (MQ=255)                                                                                           aGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATAc                                                  <  2:60348/149‑1 (MQ=255)                                                                                                gAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCtt                                             >  8:89929/1‑149 (MQ=255)                                                                                                gAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCtt                                             <  8:268254/149‑1 (MQ=255)                                                                                                                  gtgggCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCg                                      >  2:351791/1‑139 (MQ=255)                                                                                                                  gtgggCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCg                                      <  1:351791/139‑1 (MQ=255)                                                                                                                   tgggCTTTTTCGCAGCAAAGAAAAGCTCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCt                            >  1:4670/1‑149 (MQ=255)                                                                                                                     ggCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTgg                          <  7:217370/148‑1 (MQ=255)                                                                                                                       cTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGgcgc                       <  2:4670/149‑1 (MQ=255)                                                                                                                                  caAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCAc               <  5:111956/146‑1 (MQ=255)                                                                                                                                  caAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCAc               >  6:111956/1‑146 (MQ=255)                                                                                                                                   aaaGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGACATCGGTCTCTACCTGGCGCAGTTGCACcag            >  8:239495/1‑148 (MQ=255)                                                                                                                                     aGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACcagcag         >  4:178250/1‑149 (MQ=255)                                                                                                                                           gC‑CCCGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcagca                                                             <  4:77007/91‑1 (MQ=255)                                                                                                                                           gC‑CCCGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTgcagca                                                             >  3:77007/1‑91 (MQ=255)                                                                                                                                              cccGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTcc  >  2:225059/1‑148 (MQ=255)                                                                                                                                                   aCTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATc                                                                                 <  5:133475/64‑1 (MQ=255)                                                                                                                                                   aCTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATc                                                                                 >  6:133475/1‑64 (MQ=255)                                                                                                                                                   |                                                                                                                                               AAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCT‑GGGTGGGCTTTTTCGCAGCAAAGAAAAGC‑CCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCC  >  NC_002947/1859383‑1859669

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