Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,920	C→T	6.1%	D396D (GAC→GAT) ‡	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,920	0	C	T	6.1%	249.4 / 3.9	82	D396D (GAC→GAT) ‡	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base T (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

CAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGA  >  NC_002947/5028773‑5029063                                                                                                                                                      |                                                                                                                                                caGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑GTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACt                                                                                                                                                >  5:223253/1‑149 (MQ=255)      tGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACg                                                                                                                                            >  6:145718/1‑148 (MQ=255)         tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCTAGA‑TTACAGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGATCGA‑CTCCTCGCAAGCTTTTCGTAAAGGACTGACGaa                                                                                                                                          >  6:100616/1‑147 (MQ=255)         tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGaa                                                                                                                                          >  1:252667/1‑148 (MQ=255)           aaaCCTTGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGATACCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAg                                                                                                                                       <  1:100043/148‑1 (MQ=255)             aCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGttt                                                                                                                                    <  1:112268/149‑1 (MQ=255)             aCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGctt                                                                                                                                    <  7:33832/149‑1 (MQ=255)             aCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGctt                                                                                                                                    <  7:28956/149‑1 (MQ=255)              ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTg                                                                                                                                   <  8:13664/149‑1 (MQ=255)               cTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTg                                                                                                                                   >  6:238272/1‑148 (MQ=255)               cTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGc                                                                                                                                  <  4:75233/149‑1 (MQ=255)                       aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGc                                                                                                                           <  5:17220/148‑1 (MQ=255)                       aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                          >  6:210375/1‑149 (MQ=255)                          ttCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                                                                                                                         >  6:52140/1‑147 (MQ=255)                              aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                          <  8:186299/142‑1 (MQ=255)                              aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                          >  7:186299/1‑142 (MQ=255)                              aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTta                                                                                                                   <  5:145718/149‑1 (MQ=255)                                ttCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTtata                                                                                                                 >  8:2916/1‑149 (MQ=255)                                   aTGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGt                                                                                                               <  5:45830/148‑1 (MQ=255)                                   aTGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTg                                                                                                              >  3:72237/1‑149 (MQ=255)                                      ccTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGcc                                                                                                            <  2:257038/148‑1 (MQ=255)                                                 gatggatgATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCtt                                                                                                 <  2:252667/148‑1 (MQ=255)                                                  atggatgATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCtt                                                                                                 >  8:284561/1‑147 (MQ=255)                                                    ggatgATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc                                                                                             >  6:253900/1‑149 (MQ=255)                                                      atgatTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                          <  3:46543/118‑1 (MQ=255)                                                      atgatTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                          >  4:46543/1‑118 (MQ=255)                                                         atTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACaa                                                                                        >  4:9487/1‑149 (MQ=255)                                                          ttGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAt                                                                                       <  3:9487/149‑1 (MQ=255)                                                          ttGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAt                                                                                       <  5:52140/149‑1 (MQ=255)                                                          ttGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAt                                                                                       <  5:238272/149‑1 (MQ=255)                                                             ggggCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                                                                                                     >  2:12373/1‑100 (MQ=255)                                                             ggggCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                                                                                                                     <  1:12373/100‑1 (MQ=255)                                                                gCTCGAGT‑TTACCGAGAGCATCCTGATCACTTATACTGGGGTCTAGACGTTCTGCTCGA‑CTCCTCGCAATCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATg                                                                                 <  5:100616/149‑1 (MQ=255)                                                                gCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTTCCACGGCTATCTTCAGCGACAATGTGATg                                                                                 <  1:204071/149‑1 (MQ=255)                                                                  tCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTg                                                                                                              <  8:84455/118‑1 (MQ=255)                                                                  tCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTg                                                                                                              >  7:84455/1‑118 (MQ=255)                                                                  tCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGt                                                                               >  4:104587/1‑149 (MQ=255)                                                                    gaga‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGAc                                                                                          >  4:201741/1‑136 (MQ=255)                                                                    gaga‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGAc                                                                                          <  3:201741/136‑1 (MQ=255)                                                                    gaga‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGc                                                                             >  8:197080/1‑149 (MQ=255)                                                                    gaga‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGc                                                                             <  5:210375/149‑1 (MQ=255)                                                                          tACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTa                                                                         <  4:116803/148‑1 (MQ=255)                                                                             cGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaaa                                                                                                                                        >  2:153505/1‑82 (MQ=255)                                                                             cGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaaa                                                                                                                                        <  1:153505/82‑1 (MQ=255)                                                                               agagCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCAc                                                                   <  6:223253/149‑1 (MQ=255)                                                                                        tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTAAGCGACAATTTGATGGTGACTTAAAGCACAATCACtac                                                          >  7:219426/1‑147 (MQ=255)                                                                                             aCTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCa                                                     >  7:36061/1‑149 (MQ=255)                                                                                               tGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAg                                                                                              <  7:77120/106‑1 (MQ=255)                                                                                               tGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAg                                                                                              >  8:77120/1‑106 (MQ=255)                                                                                               tGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACt                                                   >  6:286136/1‑149 (MQ=255)                                                                                                 aTACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTcc                                                 >  1:149719/1‑149 (MQ=255)                                                                                                  tACTGGGGTCGAGACGTTCTGCTCGAGCTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTcc                                                 >  1:164073/1‑149 (MQ=255)                                                                                                       gggTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtt                                           <  3:104587/149‑1 (MQ=255)                                                                                                         gTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGtttt                                         <  2:164073/149‑1 (MQ=255)                                                                                                          tCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTa                                        >  2:221424/1‑149 (MQ=255)                                                                                                            gagaCGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCtt                                               <  8:84272/140‑1 (MQ=255)                                                                                                            gagaCGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCtt                                               >  7:84272/1‑140 (MQ=255)                                                                                                                      gCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCAc                            >  7:40231/1‑149 (MQ=255)                                                                                                                       cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc                                                  <  6:251667/126‑1 (MQ=255)                                                                                                                       cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc                                                  >  5:251667/1‑126 (MQ=255)                                                                                                                       cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGc                              <  5:253900/146‑1 (MQ=255)                                                                                                                       cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCa                             <  8:36061/147‑1 (MQ=255)                                                                                                                       cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACg                           >  3:80276/1‑149 (MQ=255)                                                                                                                        tCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACg                           >  8:137562/1‑148 (MQ=255)                                                                                                                          gA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTGAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAg                        >  7:170479/1‑149 (MQ=255)                                                                                                                                ctcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCtt                    >  4:123564/1‑148 (MQ=255)                                                                                                                                 tcGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCTAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                   >  1:263448/1‑148 (MQ=255)                                                                                                                                   gCAAACTTTTCGTCAGTCCTTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                             >  7:106173/1‑120 (MQ=38)                                                                                                                                   gCAAACTTTTCGTCAGTCCTTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                             <  8:106173/120‑1 (MQ=38)                                                                                                                                   gCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg                   >  2:197425/1‑146 (MQ=39)                                                                                                                                    cAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAACATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa                 >  5:31935/1‑147 (MQ=255)                                                                                                                                       gCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg             >  8:265573/1‑148 (MQ=255)                                                                                                                                        cTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGtt           <  4:72237/149‑1 (MQ=255)                                                                                                                                        cTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGtt           <  3:123564/149‑1 (MQ=255)                                                                                                                                             cGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGaaa       <  2:149719/148‑1 (MQ=255)                                                                                                                                              gTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAg      <  8:40231/148‑1 (MQ=255)                                                                                                                                              gTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGt     <  7:2916/149‑1 (MQ=255)                                                                                                                                               tCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGtt    >  2:195203/1‑149 (MQ=255)                                                                                                                                                cAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGtt    >  6:206667/1‑148 (MQ=255)                                                                                                                                                 aaGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTTTCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa  >  4:282529/1‑149 (MQ=255)                                                                                                                                                 aaGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa  <  2:263448/149‑1 (MQ=255)                                                                                                                                                 aaGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGa  <  5:286136/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                CAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGA  >  NC_002947/5028773‑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: ‑