Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,917	G→T	7.0%	K395N (AAG→AAT)	PP_4432 →	Xaa‑Pro aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,028,917	0	G	T	7.0%	185.5 / 3.7	66	K395N (AAG→AAT)	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base T (2/3);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

CAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATT  >  NC_002947/5028770‑5029065                                                                                                                                                     |                                                                                                                                                                caacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAg                                                                                                                                                                 >  2:75642/1‑148 (MQ=255)  aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGa                                                                                                                                                               <  3:19159/149‑1 (MQ=255)  aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGCGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGa                                                                                                                                                               >  4:194863/1‑149 (MQ=255)       cGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGAc                                                                                                                                                          <  3:128685/149‑1 (MQ=255)            tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGATACGTTATGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGaa                                                                                                                                                       >  6:81812/1‑147 (MQ=255)              aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGt                                                                                                                                                  >  2:268680/1‑149 (MQ=255)              aaaCATGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGTGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGaaa‑a                                                                                                                                                    <  7:214920/147‑1 (MQ=255)               aaCCTGGTATGACCTTCCTCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGtt                                                                                                                                                 <  2:210443/149‑1 (MQ=255)               aaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGtt                                                                                                                                                 <  1:183333/149‑1 (MQ=255)                aCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCGCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGttt                                                                                                                                                <  2:253217/149‑1 (MQ=255)                aCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGAGTGGCGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGttt                                                                                                                                                <  8:217274/149‑1 (MQ=255)                aCCTGGTATGAACTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑ACTTTTCGTCAGTCCTTGACGAAA‑AGctt                                                                                                                                                <  8:122534/149‑1 (MQ=255)                 ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTg                                                                                                                                               <  6:254897/149‑1 (MQ=255)                     gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCaa                                                                                                                                            <  6:263567/148‑1 (MQ=255)                     gTATGACCTTCCACTTCATCCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑ACTTTTCGTCAGTCCTTGACGAAA‑AGCTTGCGAg                                                                                                                                           <  7:115737/149‑1 (MQ=38)                          aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACTTTATGCTCGACTCCTCGCAA‑GATTTTCGTAAAGGACTGACGAAA‑AGTTTGCAATTCGCt                                                                                                                                      >  3:169492/1‑149 (MQ=255)                           ccTTCCACTTCATGCCTGGGCTATGGATGGATGATTGGCGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGAAGAAA‑AGTTTGCAAGTCGCTg                                                                                                                                     <  4:54774/149‑1 (MQ=255)                            cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTg                                                                                                                                     >  8:199306/1‑148 (MQ=255)                                  cTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTtat                                                                                                                     <  1:308586/149‑1 (MQ=255)                                  cTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTta                                                                                                                      >  8:84814/1‑148 (MQ=255)                                       tGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGc                                                                                                                <  2:125203/149‑1 (MQ=255)                                          cTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCAc                                                                                                             <  5:116306/149‑1 (MQ=255)                                             gCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACgg                                                                                                           <  5:269622/148‑1 (MQ=255)                                               cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAATTCGCTCG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCt                                                                                                         >  4:72085/1‑148 (MQ=255)                                               cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTa                                                                                                        >  4:76482/1‑149 (MQ=255)                                                tATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTAt                                                                                                       >  2:100318/1‑149 (MQ=255)                                                  tggatggaTGATTGGTGGCTCTAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAACGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTAt                                                                                                       <  5:81812/147‑1 (MQ=255)                                                  tggatggaTGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCt                                                                                                     >  1:290172/1‑149 (MQ=255)                                                      tggatgATTGGGGGCTCGAGATTACCGAGCGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAg                                                                                                 <  6:290272/149‑1 (MQ=255)                                                      tggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAg                                                                                                 >  3:257160/1‑149 (MQ=255)                                                          tgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGAc                                                                                             <  7:84814/149‑1 (MQ=255)                                                                 tgtcTCGAGATTACCGAGATCATCCGGATGACTGATACTGGTGTCGAGTCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTGAATGACTGACGAAA‑AGTTTGAAAGTCTCTTG‑‑‑‑‑‑‑‑‑ATTTATAGTTCGACGGCTATCTTCAGCGACAATGTGa                                                                                      <  6:294180/146‑1 (MQ=255)                                                                       gagaTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGc                                                                                <  3:313909/149‑1 (MQ=255)                                                                       gagaTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGc                                                                                >  4:313909/1‑149 (MQ=255)                                                                         gaTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCt                                                                              >  2:37459/1‑149 (MQ=255)                                                                              ccGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGACGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTaaa                                                                          <  3:194863/148‑1 (MQ=255)                                                                                gagaGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCtt                                                                                                    >  8:87949/1‑120 (MQ=255)                                                                                gagaGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCtt                                                                                                    <  7:87949/120‑1 (MQ=255)                                                                                     cATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGGCTTTTCTTCAAGGACTGACGAAA‑AGTTTGAAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCAACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGAACCAt                                                                   >  1:248774/1‑149 (MQ=255)                                                                                       tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACCATCAc                                                                <  4:189807/149‑1 (MQ=255)                                                                                       tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAc                                                                >  4:215439/1‑149 (MQ=255)                                                                                           gATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGa                                                                                              >  6:110887/1‑115 (MQ=255)                                                                                                 tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACt                                                     >  4:94445/1‑149 (MQ=255)                                                                                                 tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACt                                                     >  6:102999/1‑149 (MQ=255)                                                                                                       tGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTg                                               >  5:266793/1‑149 (MQ=255)                                                                                                         gggTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGtt                                             <  3:94445/149‑1 (MQ=255)                                                                                                         gggTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGtt                                             <  1:37459/149‑1 (MQ=255)                                                                                                              gagaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTAc                                         >  5:174132/1‑148 (MQ=255)                                                                                                              gagaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTAc                                         >  4:116902/1‑148 (MQ=255)                                                                                                                gaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCCTCTCCCTTGCCCCCGCCTAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCttt                                                <  1:268680/148‑1 (MQ=255)                                                                                                                gaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCa                                                           <  7:132806/128‑1 (MQ=255)                                                                                                                gaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCa                                                           >  8:132806/1‑128 (MQ=255)                                                                                                                     tCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACaaaaaaa                                  <  3:184051/148‑1 (MQ=255)                                                                                                                      cTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATTGTGCCTTAAAACACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAg                                 <  3:76482/148‑1 (MQ=255)                                                                                                                      cTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAACACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAg                                 <  3:72085/148‑1 (MQ=255)                                                                                                                            gACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAAGAGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCCCGGCTATGTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACg                             >  1:53113/1‑147 (MQ=255)                                                                                                                                  ttgCAA‑ACTTTTCGTCAGTCCTTGACGAAA‑AGCTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                     >  2:149420/3‑148 (MQ=38)                                                                                                                                  ttgCAA‑ACTTTTCGTCAGTCCTTGACGAAA‑AGCTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                     >  1:159111/3‑148 (MQ=38)                                                                                                                                    gCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTa                                          >  7:193758/1‑125 (MQ=255)                                                                                                                                    gCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTa                                          <  8:193758/125‑1 (MQ=255)                                                                                                                                    gCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAag                  <  3:116902/149‑1 (MQ=255)                                                                                                                                      aa‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                >  3:284781/1‑149 (MQ=255)                                                                                                                                      aa‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                <  3:317892/149‑1 (MQ=255)                                                                                                                                      aa‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag                <  6:174132/149‑1 (MQ=255)                                                                                                                                         gCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg               <  1:75642/149‑1 (MQ=255)                                                                                                                                              tCGTCAAGGACTGACGAAA‑AGGTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGGGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaa          <  6:236547/148‑1 (MQ=255)                                                                                                                                                   aaGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa    <  2:290172/149‑1 (MQ=255)                                                                                                                                                     ggACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGAtt  <  2:53113/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                                CAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATT  >  NC_002947/5028770‑5029065

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