Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	733,647	A→T	5.4%	intergenic (+21/‑26)	PP_t09 → / → PP_t10	tRNA‑Lys/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	733,647	0	A	T	5.4%	223.4 / 2.6	74	intergenic (+21/‑26)	PP_t09/PP_t10	tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base T (2/2);  total (38/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.06e-01

GGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAC  >  NC_002947/733504‑733793                                                                                                                                                 |                                                                                                                                                   ggAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGGAGAAAAGGGCCTCTCAAAGAGAt                                                                                                                                               >  4:242888/1‑149 (MQ=255)  gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATg                                                                                                                                              >  7:217813/1‑149 (MQ=255)  gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATg                                                                                                                                              <  7:131692/149‑1 (MQ=255)    gCAGGAGGTCAGGAGTTCGAATCCCCTCCCCCCGACCATTTTCCGGGTCGGTAGCTCAGGTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑GCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGcc                                                                                                                                            <  3:242888/149‑1 (MQ=17)       ggATGTCAGGAGTTCGAATCCCCTCACCACGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCATACACGACCCACCATGTAGAAAATGGCATCTCAAAGAGATGCCAtt                                                                                                                                         >  8:347633/1‑149 (MQ=255)        gATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttt                                                                                                                                        <  8:217813/149‑1 (MQ=255)              aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCttttctaca                                                                                                                                  <  6:202440/149‑6 (MQ=255)              aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttttt                                                                                                                                   >  5:336135/1‑148 (MQ=255)              aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGAATCTCAAAGAGATGCCCtttttttt                                                                                                                                   >  5:341911/1‑148 (MQ=255)               ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTg                                                                                                                                  >  6:275814/1‑148 (MQ=255)               ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTg                                                                                                                                  >  3:181376/1‑148 (MQ=255)                gAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaa                                                                                                                                <  2:79706/149‑1 (MQ=255)                 aGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTCAGGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttttttga                                                                                                                                <  1:111773/149‑4 (MQ=17)                 aGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATTCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaaa                                                                                                                               <  1:103022/149‑1 (MQ=255)                   ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                            >  5:358741/1‑134 (MQ=255)                   ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                            <  6:358741/134‑1 (MQ=255)                       aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                         <  7:48958/149‑1 (MQ=255)                       aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                         <  5:148391/149‑1 (MQ=255)                         tCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                       <  2:177911/149‑1 (MQ=255)                           cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                     <  7:289799/149‑1 (MQ=255)                             cTCACCCCGACCATTTTCCGGGTGGTTAGCTCAGTTGGTAGAGCAGTGGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtgt                                                                                                                   <  6:299183/149‑1 (MQ=255)                              tCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtg                                                                                                                    >  8:90113/1‑147 (MQ=255)                                  cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAgcg                                                                                                               >  8:316090/1‑148 (MQ=255)                                   ccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCATTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                             <  7:347633/149‑1 (MQ=255)                                    cGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAg                                                                                                            >  5:366082/1‑149 (MQ=255)                                               cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgc                                                                                                 <  8:200533/149‑1 (MQ=255)                                               cGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgc                                                                                                 >  7:200533/1‑149 (MQ=255)                                                      ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCt                                                                                          >  6:72874/1‑149 (MQ=255)                                                         gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCttt                                                                                        <  5:275814/148‑1 (MQ=255)                                                          cTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                      >  7:321737/1‑149 (MQ=255)                                                             aGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGC‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAg                                                                                   >  4:164140/1‑149 (MQ=255)                                                                 ggTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                <  7:90113/148‑1 (MQ=255)                                                                 ggTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                <  7:316090/148‑1 (MQ=255)                                                                  gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGATACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGa                                                                              >  7:31602/1‑149 (MQ=255)                                                                     gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGt                                                                           <  4:181376/149‑1 (MQ=255)                                                                     gagCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGt                                                                           <  6:336135/149‑1 (MQ=255)                                                                     gagCAGGTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGAAGGATGt                                                                           <  6:341911/149‑1 (MQ=255)                                                                       gCAGTTGGCTTTTAACCAATTGGTCGTAGGC‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCa                                                                         <  3:164140/149‑1 (MQ=255)                                                                           ttGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAg                                                                     >  4:349049/1‑149 (MQ=255)                                                                              gCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTc                                                                  <  5:72874/149‑1 (MQ=255)                                                                                ttttAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCg                                                                 >  6:329164/1‑148 (MQ=255)                                                                                ttttAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCg                                                                 >  2:118136/1‑148 (MQ=255)                                                                                   tAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAt                                                              <  8:321737/148‑1 (MQ=255)                                                                                    aaCCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATc                                                             <  6:366082/148‑1 (MQ=255)                                                                                        aaTTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCt                                                         >  1:131535/1‑148 (MQ=255)                                                                                                 tAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAc                                                >  3:345502/1‑148 (MQ=255)                                                                                                     t‑tCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAttt                                           >  7:193167/1‑149 (MQ=255)                                                                                                        cGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAtttt                                          <  2:131535/148‑1 (MQ=255)                                                                                                         gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCCGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTg                                        <  8:31602/149‑1 (MQ=255)                                                                                                          aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCTAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                       >  6:45305/1‑149 (MQ=255)                                                                                                          aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGAt                                                                             <  8:296454/111‑1 (MQ=255)                                                                                                          aaTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGAt                                                                             >  7:296454/1‑111 (MQ=255)                                                                                                                    gACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCt                             >  6:183830/1‑148 (MQ=255)                                                                                                                       ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAg                          >  8:241598/1‑149 (MQ=255)                                                                                                                          ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTg                       >  7:271178/1‑149 (MQ=255)                                                                                                                          ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCATTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTg                       >  1:234582/1‑149 (MQ=255)                                                                                                                            aTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                      >  8:22009/1‑148 (MQ=255)                                                                                                                            aTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                     <  5:329164/149‑1 (MQ=255)                                                                                                                            aTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                     <  1:58353/149‑1 (MQ=255)                                                                                                                              aaaaaaaaGGGCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATGCCCTCACCCAGACCATTTTTGGGTCGTTAGCTAAGTTgt                      >  8:73738/5‑145 (MQ=18)                                                                                                                               tAGAAAAGGGCATCTCTAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTag                   <  5:45305/148‑1 (MQ=255)                                                                                                                               tAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTaga                  <  7:241598/149‑1 (MQ=255)                                                                                                                                 gAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGTCGTTTTTTGGGTGGTTTGGTGTGTTGGttgtgg                >  2:146619/1‑144 (MQ=18)                                                                                                                                   aaaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAg              >  6:172563/1‑149 (MQ=255)                                                                                                                                       ggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg          <  5:172563/149‑1 (MQ=255)                                                                                                                                       ggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATACCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg          >  3:327150/1‑149 (MQ=255)                                                                                                                                        gCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGc         <  8:193167/149‑1 (MQ=255)                                                                                                                                         cATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAt                                                              >  3:219060/1‑95 (MQ=255)                                                                                                                                         cATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAt                                                              <  4:219060/95‑1 (MQ=255)                                                                                                                                          aTCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATc                                                             >  6:369273/1‑95 (MQ=255)                                                                                                                                          aTCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATc                                                             <  5:369273/95‑1 (MQ=255)                                                                                                                                             tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGCCCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt     >  1:328942/1‑148 (MQ=255)                                                                                                                                             tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt     >  1:323648/1‑148 (MQ=255)                                                                                                                                             tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt     >  5:288541/1‑148 (MQ=255)                                                                                                                                               aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCTGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCATTAGCTCAGTTGGTAGAGCAGTTGGCTTTTa    <  1:118136/147‑1 (MQ=255)                                                                                                                                               aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCTGTTGGTAGAGCAGTTGGCTTTTAAc  >  4:230697/1‑149 (MQ=255)                                                                                                                                               aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg           <  1:311064/140‑1 (MQ=255)                                                                                                                                               aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTg           >  2:311064/1‑140 (MQ=255)                                                                                                                                               aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAc  >  2:314096/1‑149 (MQ=255)                                                                                                                                                aaGCGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAc  >  5:259828/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGT‑TCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAC  >  NC_002947/733504‑733793

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