Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,919	A→C	6.0%	D396A (GAC→GCC) ‡	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,919	0	A	C	6.0%	206.6 / 2.9	66	D396A (GAC→GCC) ‡	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (26/36);  new base C (2/2);  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 = 9.69e-01

AACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTG  >  NC_002947/5028771‑5029066                                                                                                                                                      |                                                                                                                                                               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                                                                                                     <  7:87949/120‑1 (MQ=255)                                                                               gagaGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCtt                                                                                                     >  8:87949/1‑120 (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                                          >  4:116902/1‑148 (MQ=255)                                                                                                             gagaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTAc                                          >  5:174132/1‑148 (MQ=255)                                                                                                               gaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCCTCTCCCTTGCCCCCGCCTAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCttt                                                 <  1:268680/148‑1 (MQ=255)                                                                                                               gaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCa                                                            >  8:132806/1‑128 (MQ=255)                                                                                                               gaCGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCa                                                            <  7:132806/128‑1 (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)                                                                                                                                                     gACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTg  >  7:177147/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                               AACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAA‑GCTTTTCGTCAAGGACTGACGAAA‑AGTTTGCAAGTCGCTGG‑‑‑‑‑‑‑‑‑AGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAC‑TCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGATTG  >  NC_002947/5028771‑5029066

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