Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	101,284	A→G	6.3%	N35S (AAC→AGC)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	101,284	0	A	G	6.3%	182.5 / 2.6	64	N35S (AAC→AGC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

GATGGGGCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGTGTGCAACAGTGCCGA  >  NC_002947/101144‑101431                                                                                                                                                 |                                                                                                                                                    gATGGGGCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTaaa                                                                                                                                              <  6:43167/148‑1 (MQ=255) gATGGGGCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTa                                                                                                                                                >  5:43167/1‑146 (MQ=255)  aTGGGGCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGc                                                                                                                                            <  3:151814/149‑1 (MQ=255)     gggCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCAt                                                                                                                                         <  7:180618/149‑1 (MQ=255)       gCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTg                                                                                                                                       <  8:51090/149‑1 (MQ=255)         aGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGTTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCAt                                                                                                                                         <  8:15251/145‑1 (MQ=255)         aGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGTTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCAt                                                                                                                                         >  7:15251/1‑145 (MQ=255)           tgtTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTCCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACACACGTAAAGCCATTGcc                                                                                                                                     >  7:255068/1‑147 (MQ=255)                   ctGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTCCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCt                                                                                                                             <  8:255068/147‑1 (MQ=255)                    tGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCg                                                                                                                           >  5:253739/1‑148 (MQ=255)                     gCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCga                                                                                                                                   >  6:263422/1‑139 (MQ=255)                     gCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCga                                                                                                                                   <  5:263422/139‑1 (MQ=255)                     gCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGa                                                                                                                          >  6:56369/1‑148 (MQ=255)                      ccAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGc                                                                                                                                      >  8:196087/1‑135 (MQ=255)                      ccAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGc                                                                                                                                      <  7:196087/135‑1 (MQ=255)                          ggTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAg                                                                                                                      >  8:256262/1‑147 (MQ=255)                                cAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCa                                                                                                               <  6:253739/148‑1 (MQ=255)                                 aaaCCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAACTAAGTTTTTGCCGGCGTCGTCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAggg                                                                                                                 >  1:256416/1‑149 (MQ=255)                                 aaaCCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACaa                                                                                                                                                     <  5:187738/109‑1 (MQ=255)                                 aaaCCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACaa                                                                                                                                                     >  6:187738/1‑109 (MQ=255)                                 aaaCCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGc                                                                                                                >  1:72936/1‑146 (MQ=255)                                      gTGACTGCGAACAACCCGCTTCTGCAATCGCACGATCTGCCGCCCTTCTCGCAAATTCTTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAAccc                                                                                                        <  2:86997/149‑1 (MQ=255)                                           tGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTaaa                                                                                                                                              <  7:193218/106‑1 (MQ=255)                                           tGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTaaa                                                                                                                                              >  8:193218/1‑106 (MQ=255)                                             cGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAGCCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATgg                                                                                                 >  4:122691/1‑149 (MQ=255)                                             cGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATgg                                                                                                 >  3:21575/1‑149 (MQ=255)                                                  aacCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGggccgg                                                                                            >  3:94296/1‑149 (MQ=255)                                                   acCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGggccggc                                                                                           <  4:63060/149‑1 (MQ=255)                                                    cccGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGggccggcc                                                                                          >  8:160639/1‑149 (MQ=255)                                                       gCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTgg                                                                                       <  7:256262/149‑1 (MQ=255)                                                       gCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTgg                                                                                       <  7:160639/149‑1 (MQ=255)                                                            tGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTg                                                                                  >  5:315018/1‑149 (MQ=255)                                                            tGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTg                                                                                  >  6:302634/1‑149 (MQ=255)                                                              cAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTgg                                                                                 <  5:302634/148‑1 (MQ=255)                                                               aaTCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAGCCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGc                                                                                >  1:239319/1‑148 (MQ=255)                                                                  cccACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATg                                                                            >  5:201374/1‑149 (MQ=255)                                                                    cACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGaa                                                                                                           >  8:195288/1‑116 (MQ=255)                                                                    cACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGaa                                                                                                           >  8:190689/1‑116 (MQ=255)                                                                    cACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGaa                                                                                                           <  7:195288/116‑1 (MQ=255)                                                                    cACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGaa                                                                                                           <  7:190689/116‑1 (MQ=255)                                                                             ccgccCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAGCCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAAc                                                                     <  3:122691/145‑1 (MQ=255)                                                                             ccgccCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAAc                                                                     <  5:56369/145‑1 (MQ=255)                                                                             ccgccCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGaa                                                                 <  5:170501/149‑1 (MQ=255)                                                                             ccgccCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGa                                                                  <  6:201374/148‑1 (MQ=255)                                                                             ccgccCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGa                                                                  <  2:256416/148‑1 (MQ=255)                                                                                cccTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGa                                                              >  2:181206/1‑149 (MQ=255)                                                                                   ttCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCg                                                           >  7:78266/1‑149 (MQ=255)                                                                                      tcGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACAt                                                                                                   >  2:172528/1‑106 (MQ=255)                                                                                      tcGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACAt                                                                                                   <  1:172528/106‑1 (MQ=255)                                                                                      tcGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCt                                                        >  8:108876/1‑149 (MQ=255)                                                                                                 tGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAGCCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCt                                             <  2:239319/149‑1 (MQ=255)                                                                                                   ccGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTTGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGctccaggcgg                                           >  4:28516/1‑143 (MQ=255)                                                                                                   ccGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGctccaggcgg                                           >  6:77450/1‑143 (MQ=255)                                                                                                    cGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTg                                                                                        <  4:104074/103‑1 (MQ=255)                                                                                                    cGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTg                                                                                        >  3:104074/1‑103 (MQ=255)                                                                                                      aaC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGc                                        >  8:45404/1‑149 (MQ=255)                                                                                                          acG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGccag                                     >  8:270761/1‑147 (MQ=255)                                                                                                                                 gCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGgt                 <  2:272507/149‑1 (MQ=255)                                                                                                                                   gATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGgtgt               <  4:21575/149‑1 (MQ=255)                                                                                                                                    aTCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGtgtg              >  6:238528/1‑149 (MQ=255)                                                                                                                                     gccTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGtgtg              <  7:193351/147‑1 (MQ=255)                                                                                                                                        tGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGTGTGCaa           >  8:60778/1‑148 (MQ=255)                                                                                                                                         ggCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGgt                 >  4:179733/1‑141 (MQ=255)                                                                                                                                         ggCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGgt                 <  3:179733/141‑1 (MQ=255)                                                                                                                                          gCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGTGTGCAACa         <  6:315018/148‑1 (MQ=255)                                                                                                                                             gACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCTCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGTTTGCAAAAGTGc     <  8:78266/149‑1 (MQ=255)                                                                                                                                                aaCCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGTGTGCAACAGTGCCGa  >  8:101355/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GATGGGGCAGTGTTTCTTCTGCCAAGGTATCCAAACCGTGACTGCGAACAACCCGCTTCTGCAATCCCACGATCTGCCGCCCTTCTCGCAAATTCGTGCCGAAC‑ACG‑TGTTGCCGGCG‑‑ATCGAAGCGATCCTGGCCGACAACCGTAAAGCCATTGCCGAGATCCTCGAAAAGCAGGGCAAGAACCCTACATGGGCCGGCCTGGTGCTGGCCATGGACGAACTGAACGACCGCCTGGGCGCCGCCTGGAGCCCGGTCAGCCACCTCAATGCGGTGTGCAACAGTGCCGA  >  NC_002947/101144‑101431

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑