BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,028,905	G→A	5.2%	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	5.2%	349.6 / 4.0	113	K391K (AAG→AAA)	PP_4432	Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand): ref base G (57/50); new base A (3/3); total (60/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 114 total aligned reads displayed.

TGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGG  >  NC_002947/5028757‑5029052

tGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAg                                                                                                                                                       >  5:186246/1‑149 (MQ=255)

    cccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGGCTGGCGTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCttt                                                                                                                                                   <  6:70730/149‑1 (MQ=255)

      cGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTc                                                                                                                                                 >  2:22844/1‑149 (MQ=255)

           gACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATAATGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTc                                                                                                                                              >  8:100251/1‑147 (MQ=255)

              aacaGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGa                                                                                                                                        >  7:156554/1‑149 (MQ=255)

                   cGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGAc                                                                                                                                   <  2:406866/149‑1 (MQ=255)

                     tGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGa                                                                                                                                 <  5:179442/149‑1 (MQ=255)

                       cTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCaa                                                                                                                                            >  2:268996/1‑137 (MQ=255)

                       cTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCaa                                                                                                                                            <  1:268996/137‑1 (MQ=255)

                       cTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGaaa                                                                                                                               <  2:230612/149‑1 (MQ=255)

                       cTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGaa                                                                                                                                >  2:307486/1‑148 (MQ=255)

                        tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGaa                                                                                                                                >  6:412809/1‑147 (MQ=255)

                          aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGaaaa                                                                                                                              <  1:144752/147‑1 (MQ=255)

                          aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGt                                                                                                                            >  3:339150/1‑149 (MQ=255)

                            aCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGttt                                                                                                                          >  8:14768/1‑149 (MQ=255)

                            aCCTGGTATGACCTTCCACTTCATGCCTGGCCGATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACTAAAAGttt                                                                                                                          >  1:54711/1‑149 (MQ=255)

                              cTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGc                                                                                                                        <  4:339150/149‑1 (MQ=255)

                                 gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAg                                                                                                                     <  2:277467/149‑1 (MQ=255)

                                    tGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAACTTTTCGTCAG‑TCCTTGACGAAAAGCTTGCGAgga                                                                                                                   <  1:97408/148‑3 (MQ=38)

                                     gACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGc                                                                                                                 >  8:210724/1‑149 (MQ=255)

                                      aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGTGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCt                                                                                                                <  4:291124/149‑1 (MQ=255)

                                      aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGc                                                                                                                 <  3:49153/148‑1 (MQ=255)

                                        cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACATTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTg                                                                                                               >  4:281505/1‑148 (MQ=255)

                                          tCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAg                                                                                                            >  2:329050/1‑149 (MQ=255)

                                          tCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAATCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGGTTGCAAGTCGCTGGAg                                                                                                            >  3:319614/1‑149 (MQ=255)

                                           ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAg                                                                                                            <  6:186246/148‑1 (MQ=255)

                                             ccTTCATGCCTGGCCTATCGATCGATCATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTCACGAAAAGTTTGCAAGTAGCTGGAGTta                                                                                                         >  6:201170/2‑149 (MQ=255)

                                             aCTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTta                                                                                                         >  6:355258/1‑149 (MQ=255)

                                                     ccTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGGGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCa                                                                                                 <  3:228735/149‑1 (MQ=255)

                                                       tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGt                                                                                                     >  1:410351/1‑143 (MQ=255)

                                                       tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGt                                                                                                     <  2:410351/143‑1 (MQ=255)

                                                         gCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGc                                                                                             <  2:76413/149‑1 (MQ=255)

                                                           cTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTa                                                                                           >  1:42620/1‑149 (MQ=255)

                                                            tATGGATGGATGATTGGGGGCTCGAGATTACCAGAGAGCATCATGATCAATGATAATGGGGTCGAGACGTTCTGCTCGACTCATCGC‑AAGCTTTTCGTAAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTa                                                                                           >  1:53999/1‑149 (MQ=255)

                                                              tggatggaTGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAAACTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATc                                                                                         <  1:126016/149‑1 (MQ=255)

                                                              tggatggaTGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATc                                                                                         <  5:253943/148‑1 (MQ=255)

                                                              tggatggaTGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCt                                                                                        >  2:18423/1‑149 (MQ=255)

                                                                 atggatgATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa                                                                                     <  3:237973/149‑1 (MQ=255)

                                                                  tggatgATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACATTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa                                                                                     <  3:281505/148‑1 (MQ=255)

                                                                    gatgatTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGc                                                                                   >  5:406396/1‑148 (MQ=255)

                                                                     atgatTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGa                                                                                 <  2:51461/149‑1 (MQ=255)

                                                                      tgatTAGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGa                                                                                 >  8:155068/1‑148 (MQ=255)

                                                                       gatTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGAc                                                                                >  7:288015/1‑148 (MQ=255)

                                                                         ttGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAt                                                                             <  5:412809/149‑1 (MQ=255)

                                                                            ggggCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAAtgt                                                                           >  8:267401/1‑148 (MQ=255)

                                                                              ggCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGAt                                                                        <  8:156554/149‑1 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCg                                                                                                                  >  2:402463/1‑104 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCg                                                                                                                  <  1:402463/104‑1 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATg                                                                       <  6:226741/147‑1 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATg                                                                       >  5:226741/1‑147 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGt                                                                     >  7:299321/1‑149 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGAAAGTCGCTGGAGTTATAGTGAAACGGCTATCTTCAGCGACAATGTGATGGt                                                                     >  7:163670/1‑149 (MQ=255)

                                                                                 tCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGCGGTCGAGACGTTCTCCTCGAATCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCCCTGGAGTTAAAGTGCCACGGCTATCTTCAGCGACAATGTGATGGt                                                                     >  8:39486/1‑149 (MQ=255)

                                                                                   gagaTTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGt                                                                     <  5:355258/147‑1 (MQ=255)

                                                                                   gagaTTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTg                                                                    <  8:288015/148‑1 (MQ=255)

                                                                                   gagaTTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTg                                                                    <  7:267401/148‑1 (MQ=255)

                                                                                      aTTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCt                                                                 >  2:73510/1‑148 (MQ=255)

                                                                                          cc‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAg                                                            <  1:22844/149‑1 (MQ=255)

                                                                                           c‑GAGATCATCCTTATCATTGATACTGGGTGAGAGACGTTCTGATCGAGTCCTCGC‑AAGCTTTTCGTCAA‑TGACTGACGAAAAGTTTGCAAGTCGATGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGc                                                           <  4:319614/149‑1 (MQ=255)

                                                                                           c‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGAc                                                                                                                                   >  8:320480/1‑77 (MQ=255)

                                                                                           c‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGAc                                                                                                                                   <  7:320480/77‑1 (MQ=255)

                                                                                              agagCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCAc                                                         >  5:424898/1‑149 (MQ=255)

                                                                                              agagCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGACTATCTTCAGCGACAATGTGa                                                                         >  4:391888/1‑133 (MQ=255)

                                                                                              agagCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGACTATCTTCAGCGACAATGTGa                                                                         <  3:391888/133‑1 (MQ=255)

                                                                                                   aTCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCa                                                   <  2:42620/149‑1 (MQ=255)

                                                                                                    tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCAc                                                  >  2:121743/1‑149 (MQ=255)

                                                                                                    tCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCAc                                                  >  1:380090/1‑149 (MQ=255)

                                                                                                    tCCTGATCACTGATACTGGCGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCAc                                                  <  8:335278/149‑1 (MQ=255)

                                                                                                     ccTGATCACTGATACTGGGGTCGAGACTTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACC‑ATCAca                                                 <  1:86018/149‑2 (MQ=255)

                                                                                                      cTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACC‑ATCAca                                                 <  1:192947/148‑2 (MQ=255)

                                                                                                      cTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTc                                                <  5:334889/149‑1 (MQ=255)

                                                                                                       tGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTc                                                >  4:363278/1‑148 (MQ=255)

                                                                                                          tCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGc                                            <  2:380090/149‑1 (MQ=255)

                                                                                                            aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑a                                                      >  3:160974/1‑137 (MQ=255)

                                                                                                            aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑a                                                      <  8:288941/137‑1 (MQ=255)

                                                                                                            aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑a                                                      <  4:160974/137‑1 (MQ=255)

                                                                                                            aCTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑a                                                      >  7:288941/1‑137 (MQ=255)

                                                                                                             cTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTACAGCCAc                                         >  1:152387/1‑149 (MQ=255)

                                                                                                              tGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACt                                        >  4:44250/1‑149 (MQ=255)

                                                                                                                 tACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAAGGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTcc                                      >  1:320678/1‑149 (MQ=255)

                                                                                                                           gagaCGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTAc                            >  8:273442/1‑148 (MQ=255)

                                                                                                                                 ttCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACaaaaaaa                     <  1:426607/149‑1 (MQ=255)

                                                                                                                                   cTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAg                    >  8:113126/1‑148 (MQ=255)

                                                                                                                                      cTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGCCAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCATTTGTTTTACAAAAAAAGCagg                >  2:364012/1‑147 (MQ=255)

                                                                                                                                      cTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCa                  <  7:14768/147‑1 (MQ=255)

                                                                                                                                         gACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACC‑ATCAc                                                  <  7:181895/112‑1 (MQ=255)

                                                                                                                                         gACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCTTTAAGGCACC‑ATCAc                                                  >  8:181895/1‑112 (MQ=255)

                                                                                                                                          aCTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGc            >  4:205680/1‑149 (MQ=255)

                                                                                                                                            tcctcGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCt          >  6:149109/1‑149 (MQ=255)

                                                                                                                                            tcctcGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCt          >  7:237548/1‑149 (MQ=255)

                                                                                                                                              ctcGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg        <  1:329050/149‑1 (MQ=255)

                                                                                                                                                cGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg        >  2:395674/1‑147 (MQ=255)

                                                                                                                                                 gC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGa       <  2:30223/147‑1 (MQ=255)

                                                                                                                                                 gC‑AAGCTTTTAGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTTGAGTTATAGTGCCACGGCTATCTTCAGCTACAATGTGATGGTGACTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGACTTGaa      <  1:360876/148‑1 (MQ=255)

                                                                                                                                                 gC‑AAGATTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATTGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGaa      <  5:201170/148‑1 (MQ=255)

                                                                                                                                                 gC‑AAACTTTTCGTCAG‑TCCTTGACGAAAAGCTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACCTAAAAGTACAGCCCCTCCTTTGTTTTACAAAAAAAGCACGAGGCCTTg        >  1:146397/1‑147 (MQ=25)

                                                                                                                                                 gC‑AAACTTTTCGTCAG‑TCCTTGACGAAAAGCTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACA‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTg        >  5:105299/1‑146 (MQ=38)

                                                                                                                                                    aaGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAagag   <  7:273442/149‑1 (MQ=255)

                                                                                                                                                    aaGCTTTTCGTAAA‑GGACTGACGAAAAGTTTGCAAGTCGCTTGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTTATGGAGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCATCCTTGaa      <  2:54711/146‑1 (MQ=255)

                                                                                                                                                     aGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAgg  <  2:152387/149‑1 (MQ=255)

TGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACC‑GAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGC‑AAGCTTTTCGTCAA‑GGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAGCGACAATGTGATGGTGCCTTAAAGCACC‑ATCACTCCAGCCACTCCTTTGTTTTACAAAAAAAGCACGCAGCCTTGAAGAGG  >  NC_002947/5028757‑5029052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`