Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,917	G→T	5.9%	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	5.9%	259.6 / 4.1	84	K395N (AAG→AAT) ‡	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base T (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

AACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGA  >  NC_002947/5028771‑5029063                                                                                                                                                     |                                                                                                                                                   aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGa                                                                                                                                                  >  6:17220/1‑149 (MQ=255) aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGa                                                                                                                                                  >  7:13664/1‑149 (MQ=255)   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                                                                                                                          >  7:186299/1‑142 (MQ=255)                                aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                          <  8:186299/142‑1 (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                                                                             <  5:210375/149‑1 (MQ=255)                                                                      gaga‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGc                                                                             >  8:197080/1‑149 (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           <  3:123564/149‑1 (MQ=255)                                                                                                                                          cTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGtt           <  4:72237/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)                                                                                                                                                     |                                                                                                                                                   AACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGGTTGAAAGTTGA  >  NC_002947/5028771‑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: ‑