Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,919	A→C	6.3%	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.3%	196.5 / 2.9	63	D396A (GAC→GCC) ‡	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base C (2/2);  total (31/32)
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.25e-01

CGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGA  >  NC_002947/5028776‑5029063                                                                                                                                                |                                                                                                                                                  cGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGAc                                                                                                                                              >  2:98201/1‑149 (MQ=255)    gCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                                                                                                                            >  6:177844/1‑148 (MQ=255)        aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGt                                                                                                                                       >  2:274285/1‑149 (MQ=255)                     ccTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAg                                                                                                                                        <  4:224537/135‑1 (MQ=255)                         ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGttt                                                                                                                                     >  6:148878/1‑134 (MQ=255)                         ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGttt                                                                                                                                     <  5:148878/134‑1 (MQ=255)                         ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCCACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAg                                                                                                                      <  1:274285/148‑1 (MQ=255)                          cACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGtt                                                                                                                    <  4:272024/149‑1 (MQ=255)                           aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTta                                                                                                                   >  2:109388/1‑149 (MQ=255)                             ttCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTtata                                                                                                                 >  1:243198/1‑149 (MQ=255)                                aTGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGt                                                                                                               <  1:109388/148‑1 (MQ=255)                                     tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTT‑g                                                                                                                                   <  7:195318/123‑1 (MQ=38)                                     tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGCTT‑g                                                                                                                                   >  8:195318/1‑123 (MQ=38)                                             ggatggatGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCtt                                                                                                 >  6:272125/1‑149 (MQ=255)                                              gatggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCtt                                                                                                 >  7:213449/1‑148 (MQ=255)                                               atggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa                                                                                               <  2:243198/149‑1 (MQ=255)                                                 ggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCt                                                                                                                          <  2:123727/120‑1 (MQ=255)                                                  gatgatTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc                                                                                             <  8:295295/148‑1 (MQ=255)                                                      atTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACaa                                                                                        <  3:7547/149‑1 (MQ=255)                                                              cTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGAAGTTCTGCTCGACTACTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACTGCTATCTTCAGCGACAATGTGATgg                                                                                >  4:328371/1‑149 (MQ=255)                                                                     ttACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTCTCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCtt                                                                          <  1:98201/149‑1 (MQ=255)                                                                             agCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgtg                                                                                    <  5:40133/130‑1 (MQ=255)                                                                             agCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgtg                                                                                    >  6:40133/1‑130 (MQ=255)                                                                                 tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCAc                                                             <  5:272125/149‑1 (MQ=255)                                                                                    tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTc                                                           >  4:165350/1‑148 (MQ=255)                                                                                        cACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGc                                                       >  3:129754/1‑148 (MQ=255)                                                                                         aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCa                                                     <  3:165350/149‑1 (MQ=255)                                                                                         aCTGATACTGGGGTCAAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCa                                                     >  6:207328/1‑149 (MQ=255)                                                                                          cTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCAc                                                    >  1:66855/1‑149 (MQ=255)                                                                                             aTACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTcc                                                 <  5:177844/149‑1 (MQ=255)                                                                                                cTGGGGTCGAGACGTTCTGCTCAACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCttt                                              >  4:13210/1‑149 (MQ=255)                                                                                                 tGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                             >  8:155663/1‑149 (MQ=255)                                                                                                 tGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                             >  4:80978/1‑149 (MQ=255)                                                                                                    cgTCGAGACGTTCTGCTCGACTCCTAGCAAGCTTTTCTTCAAGGACTGACGAAAAGTTT‑GCAAGTCTCTGGAGTTATAGTGCCACGGCTATCTTCACCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtt                                           <  3:273501/147‑1 (MQ=255)                                                                                                      tcaagaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTa                                        <  5:207328/146‑1 (MQ=255)                                                                                                        gagaCGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTAc                                       >  3:295460/1‑148 (MQ=255)                                                                                                                  gCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAcc                                                                  <  2:109293/111‑1 (MQ=255)                                                                                                                  gCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAcc                                                                  >  1:109293/1‑111 (MQ=255)                                                                                                                   cTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgt                                                                                     >  1:334400/1‑91 (MQ=255)                                                                                                                   cTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgt                                                                                     <  2:334400/91‑1 (MQ=255)                                                                                                                   cTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCAc                            <  8:213449/148‑1 (MQ=255)                                                                                                                       aCTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGc                       >  2:132992/1‑149 (MQ=255)                                                                                                                         tcctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCt                     >  7:179656/1‑149 (MQ=255)                                                                                                                           cttgCAAACTTTTCGTCAGTCCTTGACGAAAAGCTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaaaaaaa                                >  7:95270/4‑136 (MQ=38)                                                                                                                           cttgCAAACTTTTCGTCAGTCCTTGACGAAAAGCTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACaaaaaaa                                <  8:95270/133‑1 (MQ=38)                                                                                                                           ctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCtt                    <  2:66855/148‑1 (MQ=255)                                                                                                                            tcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTcc                                                          <  3:319269/109‑1 (MQ=255)                                                                                                                            tcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTcc                                                          >  4:319269/1‑109 (MQ=255)                                                                                                                              gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAg                                                                                              >  1:95465/1‑71 (MQ=255)                                                                                                                              gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAg                                                                                              <  2:95465/71‑1 (MQ=255)                                                                                                                              gCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                   >  4:327636/1‑146 (MQ=255)                                                                                                                                aaGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag              >  4:104632/1‑149 (MQ=255)                                                                                                                                aaGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag              <  3:327636/149‑1 (MQ=255)                                                                                                                                 aGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg             <  3:13210/149‑1 (MQ=255)                                                                                                                                    ttttCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg             <  6:197426/146‑1 (MQ=255)                                                                                                                                    ttttCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg             <  6:192420/146‑1 (MQ=255)                                                                                                                                    ttttCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg             >  5:192420/1‑146 (MQ=255)                                                                                                                                    ttttCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg             >  5:197426/1‑146 (MQ=255)                                                                                                                                        cGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaaa       <  3:104632/148‑1 (MQ=255)                                                                                                                                          tCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGtt    <  7:155663/149‑1 (MQ=255)                                                                                                                                           cAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGtt    <  5:366886/148‑1 (MQ=255)                                                                                                                                            aaGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGaggttgaaggttga  >  3:169866/1‑149 (MQ=255)                                                                                                                                            aaGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa  <  4:169872/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  CGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTT‑GCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGA  >  NC_002947/5028776‑5029063

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