Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,425,668	G→A	5.2%	intergenic (+63/‑442)	PP_1249 → / → PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,425,668	0	G	A	5.2%	238.0 / 2.5	76	intergenic (+63/‑442)	PP_1249/PP_1250	DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (38/34);  new base A (2/2);  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 = 7.20e-01

ACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG  >  NC_002947/1425523‑1425808                                                                                                                                                  |                                                                                                                                              aCTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgg                                                                                                                                           <  2:995990/150‑1 (MQ=255)  cTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATTGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCggg                                                                                                                                          <  1:830058/150‑1 (MQ=255)   tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGTAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                         <  2:174616/150‑1 (MQ=255)     ccTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                       >  1:737140/1‑150 (MQ=255)     ccTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                       >  1:1359992/1‑150 (MQ=255)     ccTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                       >  1:511482/1‑150 (MQ=255)       tGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttt                                                                                                                                     >  1:1050530/1‑150 (MQ=255)             aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                               >  2:252757/1‑150 (MQ=255)               ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgc                                                                                                                             >  1:101678/1‑150 (MQ=255)                           aTCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgc                                                                                                                 <  2:101678/150‑1 (MQ=255)                           aTCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgc                                                                                                                 >  1:786860/1‑150 (MQ=255)                               ccAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgg                                                                                                                                           >  2:1097171/1‑120 (MQ=255)                               ccAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgg                                                                                                                                           <  1:1097171/120‑1 (MQ=255)                                cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGTCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGtt                                                                                                            >  2:1192134/1‑150 (MQ=255)                                  aGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgc                                                                                                          <  2:511482/150‑1 (MQ=255)                                          ggTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgcgcatacgc                                                                                                  <  2:1050530/150‑7 (MQ=255)                                               cTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCa                                                                                                                         <  2:478779/122‑1 (MQ=255)                                               cTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCa                                                                                                                         >  1:478779/1‑122 (MQ=255)                                                     cccGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCATGCCGCtt                                                                                       <  2:21195/150‑1 (MQ=255)                                                       cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                     >  2:593574/1‑150 (MQ=255)                                                            gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑Ttgtg                                                                               <  2:1117119/150‑1 (MQ=255)                                                             atAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTgtgg                                                                              <  1:671571/150‑1 (MQ=255)                                                              tAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTgtggt                                                                             <  1:193257/150‑1 (MQ=255)                                                                     ccTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGTCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAg                                                                      >  1:1540128/1‑150 (MQ=255)                                                                       tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTg                                                                    >  1:484477/1‑150 (MQ=255)                                                                       tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCggg                                                                                                                                          >  2:641129/1‑81 (MQ=37)                                                                       tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCggg                                                                                                                                          <  1:641129/81‑1 (MQ=37)                                                                        cagcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTgg                                                                   <  1:40802/150‑1 (MQ=255)                                                                               agTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGAGTCGCCAGCAGACTCAACCTGCTTGG‑GTGTGGTGGCCAAGTGGGTGACtt                                                            >  1:1599026/1‑150 (MQ=37)                                                                                 tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                       <  1:1064143/122‑1 (MQ=255)                                                                                 tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                       >  2:1064143/1‑122 (MQ=255)                                                                                  aaTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCACCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGa                                                         <  2:72979/150‑1 (MQ=255)                                                                                   aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAg                                                        <  1:252757/150‑1 (MQ=255)                                                                                     cTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGACCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                               >  2:1229236/1‑94 (MQ=255)                                                                                     cTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGACCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                               <  1:1229236/94‑1 (MQ=255)                                                                                      ttGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                     >  2:1133328/1‑150 (MQ=255)                                                                                         gCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  <  2:571418/150‑1 (MQ=255)                                                                                               cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                            <  1:1133328/150‑1 (MQ=255)                                                                                                ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGtt                                                                                                            <  1:664886/86‑1 (MQ=255)                                                                                                ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGtt                                                                                                            >  2:664886/1‑86 (MQ=255)                                                                                                  tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                       <  1:53627/105‑1 (MQ=255)                                                                                                  tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                       >  2:53627/1‑105 (MQ=255)                                                                                                           cGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                    <  2:1017701/146‑1 (MQ=255)                                                                                                           cGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                    >  1:1017701/1‑146 (MQ=255)                                                                                                            gCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCa                               >  1:261554/1‑150 (MQ=255)                                                                                                                  ggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGa                         >  2:1027766/1‑150 (MQ=255)                                                                                                                    cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg                       >  1:865740/1‑150 (MQ=255)                                                                                                                          ccAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                 >  1:1039739/1‑150 (MQ=255)                                                                                                                          ccAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                 >  2:1440487/1‑150 (MQ=255)                                                                                                                           cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTgg                                                                   <  2:647479/99‑1 (MQ=255)                                                                                                                           cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTgg                                                                   >  1:647479/1‑99 (MQ=255)                                                                                                                           cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAg                                                        >  2:1106054/1‑110 (MQ=255)                                                                                                                           cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAg                                                        <  1:1106054/110‑1 (MQ=255)                                                                                                                            aaGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                            >  1:540251/1‑121 (MQ=255)                                                                                                                            aaGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                            <  2:540251/121‑1 (MQ=255)                                                                                                                            aaGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCCAGGGGAAGTGttcgac               >  2:451288/1‑146 (MQ=255)                                                                                                                             aaacccccGGCCATGTGCATGACCGGGGGCTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  >  1:385532/4‑114 (MQ=11)                                                                                                                             aaacccccGGCCATGTGCATGACCGGGGGCTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  <  2:385532/111‑1 (MQ=11)                                                                                                                              ggCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCgg             >  1:126881/1‑150 (MQ=255)                                                                                                                               ccccccGGTCATTCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGAGTAGCAAG‑CCGCTGAATTTGTGTGGCTTGTGGTGGCCAATCGCGTGAATTGAAAGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            <  1:1192134/149‑1 (MQ=9)                                                                                                                                cccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg           <  2:737140/150‑1 (MQ=255)                                                                                                                                 ccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑CTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg          >  2:1037066/1‑150 (MQ=255)                                                                                                                                      gTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg     <  2:484477/150‑1 (MQ=255)                                                                                                                                      gTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg     <  2:1359992/150‑1 (MQ=255)                                                                                                                                       tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAggg    >  2:1627645/1‑150 (MQ=255)                                                                                                                                       tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGCGGCCAAGTGGGt                                                                 <  1:843625/89‑1 (MQ=255)                                                                                                                                       tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGCGGCCAAGTGGGt                                                                 >  2:843625/1‑89 (MQ=255)                                                                                                                                         aTGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtg  >  2:266253/1‑150 (MQ=255)                                                                                                                                           gCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                            >  2:1531507/1‑122 (MQ=255)                                                                                                                                           gCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                            <  1:1531507/122‑1 (MQ=255)                                                                                                                                            cacaTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑Ttgt                                                                                <  1:835544/69‑1 (MQ=255)                                                                                                                                            cacaTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑Ttgt                                                                                >  2:835544/1‑69 (MQ=255)                                                                                                                                              caTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  <  1:246417/97‑1 (MQ=255)                                                                                                                                              caTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  >  2:246417/1‑97 (MQ=255)                                                                                                                                                tGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGt                                                                     <  1:325145/76‑1 (MQ=255)                                                                                                                                                tGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGt                                                                     >  2:325145/1‑76 (MQ=255)                                                                                                                                                  |                                                                                                                                              ACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGG‑TTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG  >  NC_002947/1425523‑1425808

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