Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	503,813	A→G	11.4%	intergenic (+16/‑135)	PP_0414 → / → rpe	polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	503,813	0	A	G	11.4%	161.4 / 9.5	62	intergenic (+16/‑135)	PP_0414/rpe	polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (4/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACG  >  NC_002947/503665‑503960                                                                                                                                                     |                                                                                                                                                     tATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTa                                                                                                                                                      <  2:271781/149‑1 (MQ=255)      gcgAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGc                                                                                                                                                 >  6:328909/1‑149 (MQ=255)        gagaACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCct                                                                                                                                               >  1:243872/1‑149 (MQ=255)        gagaACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCct                                                                                                                                               <  6:117159/149‑1 (MQ=255)        gacaaCCCGAGCCCGACCAATGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCct                                                                                                                                               <  8:68835/146‑1 (MQ=255)                    ccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCttt                                                                                                                                   <  4:291066/149‑3 (MQ=255)                              ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAaga                                                                                                                        <  2:162503/149‑1 (MQ=39)                              ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTC‑CGGCTTTAGCCGCGAaga                                                                                                                        <  5:156295/149‑1 (MQ=255)                               cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAagag                                                                                                                       <  5:248098/149‑1 (MQ=255)                               cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTC‑CGGCTTTAGCCGCGAagag                                                                                                                       <  6:209403/149‑1 (MQ=255)                               cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTC‑CGGCTTTAGCCGCGAagag                                                                                                                       <  5:139302/149‑1 (MQ=255)                                 cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGc                                                                                                                     <  4:108233/149‑1 (MQ=255)                                 cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGc                                                                                                                     <  2:243872/149‑1 (MQ=255)                                 cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGc                                                                                                                     <  6:102790/149‑1 (MQ=255)                                           cTTGCTGATTGGGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGt                                                                                                           >  6:346239/1‑149 (MQ=255)                                              gCTGATAGCGTTCGCGGTCTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTCGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAggcaa                                                                                                         <  8:332908/148‑4 (MQ=25)                                               cTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGc                                                                                                        <  2:291745/148‑1 (MQ=255)                                                      cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                 >  2:311441/1‑148 (MQ=255)                                                      cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                 <  2:229997/148‑1 (MQ=255)                                                              gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTC‑CGGCTTTAGcc                                                                                                                               <  5:88727/110‑1 (MQ=38)                                                              gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTACCTGTGCCGGCCTCTTC‑CGGCTTTAGcc                                                                                                                               >  6:88727/1‑110 (MQ=38)                                                                   gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                    <  5:328909/148‑1 (MQ=255)                                                                   gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                   >  3:296236/1‑149 (MQ=255)                                                                       ctgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCaa                                                                               <  6:206400/149‑1 (MQ=255)                                                                               ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTACAGCCG‑GACGAGGCCGGCACAGGTGGCGTTGGTCGCCACCCATTGCTACCAATCCGCCtt                                                                       >  6:96909/1‑149 (MQ=255)                                                                               ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGcc                                                                                                                               >  8:275198/1‑94 (MQ=255)                                                                               ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGcc                                                                                                                               <  7:275198/94‑1 (MQ=255)                                                                               ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCtt                                                                       <  3:253501/149‑1 (MQ=255)                                                                               ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCtt                                                                       >  8:59058/1‑149 (MQ=255)                                                                               ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCt                                                                        >  5:3244/1‑148 (MQ=255)                                                                                   gAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCa                                                                     <  7:59058/147‑1 (MQ=255)                                                                                        gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                       >  1:362764/1‑124 (MQ=255)                                                                                        gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                       <  2:362764/124‑1 (MQ=255)                                                                                        gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCttt                                                              >  1:150052/1‑149 (MQ=255)                                                                                           gCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTccc                                                           >  4:358860/1‑149 (MQ=255)                                                                                            cGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCt                                                          <  5:346239/149‑1 (MQ=255)                                                                                              tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                        <  4:296236/149‑1 (MQ=255)                                                                                                      gCGAGATGCTGAGGACTCATCAGGCGTGTGGATGATGGGTGGCGTGTACGGGCGTCTTCGCGGCTTAAGCCG‑GAAGTGGCCGGCACAGGTTGCGTTGGGCGACAAGAATTGATACCAATCAGCCTTCATCGCTTTCGCTGAGTTACTAt                                                <  4:28110/149‑1 (MQ=255)                                                                                                       cGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATg                                               >  5:201794/1‑149 (MQ=255)                                                                                                               tGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAc                                       >  6:184712/1‑149 (MQ=255)                                                                                                                gAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAcc                                      >  7:136879/1‑149 (MQ=255)                                                                                                                gAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAcc                                      >  4:170194/1‑149 (MQ=255)                                                                                                                 aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCg                                     >  4:96681/1‑149 (MQ=255)                                                                                                                      tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                >  1:184988/1‑149 (MQ=255)                                                                                                                      tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                >  1:95631/1‑149 (MQ=255)                                                                                                                        aagaGCCGGCTTGATGATGGCTTGCGTGTTCCGGGCTCTTCGCGGCTAATGCCG‑GAAGAGGCCGGCAGAGGTTGCGt                                                                                                      <  1:96109/74‑1 (MQ=255)                                                                                                                        aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                            <  5:206990/135‑1 (MQ=255)                                                                                                                        aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                              <  5:184712/149‑1 (MQ=255)                                                                                                                        aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                              <  3:358860/149‑1 (MQ=255)                                                                                                                        aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCGTCATCGCTTTCCCTGAGTGACAATtcgc                                            >  6:206990/1‑135 (MQ=255)                                                                                                                        aaCAGCCGGCTTGAAGATGGCTTGCCTGTACCGGCCTCTTACCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGGTGCGt                                                                                                      >  2:96109/1‑77 (MQ=255)                                                                                                                           aGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGcc                           <  8:335568/149‑1 (MQ=255)                                                                                                                                    gatgatGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGt                    <  2:95631/147‑1 (MQ=255)                                                                                                                                    gatgatGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgc                  >  4:115776/1‑149 (MQ=255)                                                                                                                                          ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                             >  4:177228/1‑68 (MQ=255)                                                                                                                                          ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                             <  3:177228/68‑1 (MQ=255)                                                                                                                                          ggCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATg            >  5:120376/1‑149 (MQ=255)                                                                                                                                           cgcaaccTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGc           >  3:180383/6‑149 (MQ=37)                                                                                                                                             ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAg         >  2:174233/1‑149 (MQ=255)                                                                                                                                             caaccTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAg         >  6:286387/4‑149 (MQ=37)                                                                                                                                                 cTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGccc      >  6:140879/1‑148 (MQ=38)                                                                                                                                                 cTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCt     <  6:120376/149‑1 (MQ=255)                                                                                                                                                 cTGTACCGGCCTCTTACCGGCTAAAGCCG‑GAAGAGGCCCGCACAGGGTGCGTTGGTCGGCACCAAATGATACCAAACACCCTTCAACGCTTTCCCAGATTTACAAAGCGCGCCCCCGTGAAATTGCCCACAGGTGCGCCATGCAGccc      >  2:356258/1‑148 (MQ=255)                                                                                                                                                  tgtgCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTa    >  4:285355/1‑149 (MQ=38)                                                                                                                                                     aCCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACg  <  2:184988/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACG  >  NC_002947/503665‑503960

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