Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,905	G→A	6.7%	K391K (AAG→AAA)	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,905	0	G	A	6.7%	224.1 / 4.3	75	K391K (AAG→AAA)	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (2/3);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGG  >  NC_002947/5028764‑5029052                                                                                                                                                |                                                                                                                                                    gAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCg                                                                                                                                              <  5:19669/149‑1 (MQ=255)   gCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTc                                                                                                                                            <  5:289627/149‑1 (MQ=255)        aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGa                                                                                                                                       >  7:13664/1‑149 (MQ=255)        aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGa                                                                                                                                       >  6:17220/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                                                                                                               <  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                                                                            <  5:238272/149‑1 (MQ=255)                                                                   ttGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAt                                                                            <  5:52140/149‑1 (MQ=255)                                                                   ttGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAt                                                                            <  3:9487/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                                                                                                   >  7:84455/1‑118 (MQ=255)                                                                           tCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTg                                                                                                   <  8:84455/118‑1 (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                                                                                                                             <  1:153505/82‑1 (MQ=255)                                                                                      cGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaaa                                                                                                                             >  2:153505/1‑82 (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                                                                                   >  8:77120/1‑106 (MQ=255)                                                                                                        tGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAg                                                                                   <  7:77120/106‑1 (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                                    >  7:84272/1‑140 (MQ=255)                                                                                                                     gagaCGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCtt                                    <  8:84272/140‑1 (MQ=255)                                                                                                                               gCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCAc                 >  7:40231/1‑149 (MQ=255)                                                                                                                                cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc                                       >  5:251667/1‑126 (MQ=255)                                                                                                                                cTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTc                                       <  6:251667/126‑1 (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                                  <  8:106173/120‑1 (MQ=38)                                                                                                                                            gCAAACTTTTCGTCAGTCCTTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTg                                  >  7:106173/1‑120 (MQ=38)                                                                                                                                            gCAAACTTTTCGTCAGTCCTTGACGAAAAGCTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg        >  2:197425/1‑146 (MQ=39)                                                                                                                                             cAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAACATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa      >  5:31935/1‑147 (MQ=255)                                                                                                                                                gCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg  >  8:265573/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                    GAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGA‑TTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGA‑CTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGG  >  NC_002947/5028764‑5029052

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