Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,392,310	T→C	12.9%	I121M (ATA→ATG)	PP_5710 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,392,310	0	T	C	12.9%	211.8 / 20.4	85	I121M (ATA→ATG)	PP_5710	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/38);  new base C (6/5);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

AGTTCCCCCTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTGCAAGTTG  >  NC_002947/5392163‑5392450                                                                                                                                                    |                                                                                                                                              aGTTCCCCCTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTa                                                                                                                                              <  8:109325/149‑1 (MQ=255)  gTTCCCCCTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCat                                                                                                                                             >  6:213836/1‑149 (MQ=255)      ccccTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTgcatgcat                                                                                                                                         <  3:169665/149‑1 (MQ=255)       cccTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTatgcatg                                                                                                                                        >  6:175205/1‑149 (MQ=255)         cTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTtgcatgcatgca                                                                                                                                      <  7:45703/149‑1 (MQ=255)             aaaaCTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTa                                                                                                                                  >  5:149320/1‑149 (MQ=255)             aaaaCTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTa                                                                                                                                  >  6:41004/1‑149 (MQ=255)               aaCTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTgcatgcat                                                                                                                                         >  5:48512/1‑140 (MQ=255)               aaCTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTgcatgcat                                                                                                                                         <  6:48512/140‑1 (MQ=255)                aCTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGa                                                                                                                              <  5:291439/149‑1 (MQ=255)                    tGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCt                                                                                                                           <  4:56529/148‑1 (MQ=255)                      aaaCGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCAt                                                                                                                        >  7:243825/1‑149 (MQ=255)                        aCGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCAttt                                                                                                                      <  6:149320/149‑1 (MQ=255)                          gccgcATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCt                                                                                                                    >  2:43205/1‑149 (MQ=255)                           ccgcATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   <  3:331562/149‑1 (MQ=255)                           ccgcATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   >  2:275422/1‑149 (MQ=255)                           ccgcATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   <  4:174614/149‑1 (MQ=255)                           ccgcATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   <  8:41077/149‑1 (MQ=255)                            cgcATTTACTTTTAACCGTTCGAGAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCTTTTCTTg                                                                                                                  <  5:256192/149‑1 (MQ=255)                                tttACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   <  7:304357/144‑1 (MQ=255)                                tttACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   >  8:304357/1‑144 (MQ=255)                                tttACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCtt                                                                                                              >  2:1843/1‑149 (MQ=255)                                tttACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCtt                                                                                                              >  5:168520/1‑149 (MQ=255)                                   aCTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTggg                                                                                                           >  4:174923/1‑149 (MQ=255)                                   aCTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTggg                                                                                                           >  8:46273/1‑149 (MQ=255)                                   aCTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTgg                                                                                                            <  3:252994/148‑1 (MQ=255)                                   aCTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTAGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTggg                                                                                                           >  4:80634/1‑149 (MQ=255)                                    cTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGt                                                                                                          <  5:127632/149‑1 (MQ=255)                                     ttttAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTc                                                                                                         >  3:17537/1‑149 (MQ=255)                                         aaCCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTat                                                                                                                                             <  3:214027/110‑1 (MQ=255)                                         aaCCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTat                                                                                                                                             >  4:214027/1‑110 (MQ=255)                                          aCCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGt                                                                                                      >  2:260676/1‑147 (MQ=255)                                          aCCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGt                                                                                                      >  5:363863/1‑147 (MQ=255)                                            cGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGc                                                                                                  <  5:213836/149‑1 (MQ=255)                                                 gACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGc                                                                                             <  3:174923/149‑1 (MQ=255)                                                 gACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGc                                                                                             <  1:1843/149‑1 (MQ=255)                                                 gACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   >  6:103115/1‑127 (MQ=255)                                                 gACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCtt                                                                                                                   <  5:103115/127‑1 (MQ=255)                                                          aTCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTa                                                                                    >  3:136386/1‑149 (MQ=255)                                                             gcgcTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATg                                                                                  <  1:260676/148‑1 (MQ=255)                                                                   ctcGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGt                                                                           >  1:105430/1‑149 (MQ=255)                                                                   ctagCTATCTACGATGCATTCTCTGACCGTGTTCGTCGACTGGATGTAGCCTGCACCGGAAACTCTTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCa                                                                             <  3:49963/144‑1 (MQ=255)                                                                      gCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTTCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTc                                                                          <  3:80634/147‑1 (MQ=255)                                                                      gCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTc                                                                          <  6:168520/147‑1 (MQ=255)                                                                      gCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCt                                                                        <  2:105430/149‑1 (MQ=255)                                                                              cGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGt                                                                    <  4:17537/145‑1 (MQ=255)                                                                                 tCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCaa                                                              >  8:34760/1‑148 (MQ=255)                                                                                   cATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCaaaa                                                            >  4:83998/1‑148 (MQ=255)                                                                                     ttCTCTGACCGTGTTCGTCGAATGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGATCTCATTTCTTGCCTTGGGTCTGTATGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAaataa                                                         <  5:175205/149‑1 (MQ=255)                                                                                        tctGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAAtt                                                       >  3:358901/1‑148 (MQ=255)                                                                                         ctGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑c                                                                                                                                >  6:56340/1‑74 (MQ=255)                                                                                         ctGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑c                                                                                                                                <  5:56340/74‑1 (MQ=255)                                                                                         ctGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAAttt                                                      >  5:62685/1‑148 (MQ=255)                                                                                            aCCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCa                                                    <  2:177455/147‑1 (MQ=255)                                                                                            aCCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGc                                                  <  4:358901/149‑1 (MQ=255)                                                                                            aCCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATAAAGCGAGATCATTTATTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTGTTCAAAATAATTTAAg                                                   >  1:293478/1‑149 (MQ=255)                                                                                            aCCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGc                                                  <  4:136386/149‑1 (MQ=255)                                                                                                  ttCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGt                                                                                     >  2:19597/1‑108 (MQ=255)                                                                                                  ttCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGt                                                                                     <  1:19597/108‑1 (MQ=255)                                                                                                        gACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCttt                                        <  3:20935/147‑1 (MQ=255)                                                                                                        gACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCttt                                        >  7:92026/1‑147 (MQ=255)                                                                                                        gACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCttt                                        >  4:20935/1‑147 (MQ=255)                                                                                                                    ccTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          <  7:46273/149‑1 (MQ=255)                                                                                                                    ccTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          <  1:43205/149‑1 (MQ=255)                                                                                                                    ccTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          <  6:62685/149‑1 (MQ=255)                                                                                                                    ccTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          <  1:275422/149‑1 (MQ=255)                                                                                                                    ccTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          <  7:34760/149‑1 (MQ=255)                                                                                                                      tGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          <  3:266802/147‑1 (MQ=255)                                                                                                                      tGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCttt                          >  4:266802/1‑147 (MQ=255)                                                                                                                      tGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTg                         >  6:172241/1‑148 (MQ=255)                                                                                                                       gCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCt                                          <  8:30809/130‑1 (MQ=255)                                                                                                                       gCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCt                                          >  7:30809/1‑130 (MQ=255)                                                                                                                       gCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTgg                        >  8:78652/1‑148 (MQ=255)                                                                                                                           cGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCttt                                        >  7:55651/1‑128 (MQ=255)                                                                                                                           cGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCttt                                        <  8:55651/128‑1 (MQ=255)                                                                                                                            ggAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTACTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCa                                                                             >  2:113893/1‑90 (MQ=255)                                                                                                                            ggAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTACTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCa                                                                             <  1:113893/90‑1 (MQ=255)                                                                                                                             gAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCca                  >  6:238231/1‑148 (MQ=255)                                                                                                                               aaagggTCTCTGGAAATCTTGTATGCATTCATATA‑CGATCTCATTTCTTGAATTGGGTCTGTATGCGTCGCTAAGACGTATGACGCTGTCCTGGGTCTTCAAAATAATTTCTGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCcatca               <  7:306870/144‑1 (MQ=255)                                                                                                                                aCTCGTCTCTGCAAATCTTGCATGCATGCATTTA‑CGAGCTCAGTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTAAGCACGCTCCTTTGAGCCCAGCTCTTTGGAAAGcatca               >  1:124766/1‑148 (MQ=255)                                                                                                                                  tCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCg            >  7:42437/1‑149 (MQ=255)                                                                                                                                      ctctGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTATTGCCTTGGGTCTGTCAGCGTCGCTAAGCCGTATGACTCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTACTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTg         <  8:243825/148‑1 (MQ=255)                                                                                                                                       tctGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCTGCTCTTTGGAATCCATCATCGTTGc        >  1:126510/1‑148 (MQ=255)                                                                                                                                         tGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTGCAAg     >  5:98909/1‑149 (MQ=255)                                                                                                                                           cAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTGCAAGtt   <  2:293478/149‑1 (MQ=255)                                                                                                                                           cAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTGCAAGtt   >  8:334287/1‑149 (MQ=255)                                                                                                                                            aaaTCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTGCAAGTTg  <  8:42437/149‑1 (MQ=255)                                                                                                                                                cTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAaata                                                          >  2:230936/1‑89 (MQ=255)                                                                                                                                                cTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAaata                                                          <  1:230936/89‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              AGTTCCCCCTCCAAAACTATGAAACGCCGCATTTACTTTTAACCGTTCGACAGCGCCATCGCGCTCCTCGCTATCTACGATCCATTCTCTGACCGTGTTCGTCGACTGGATGGAGCCTGCACCGGAAACTCGTCTCTGCAAATCTTGTATGCATGCATTTA‑CGAGCTCATTTCTTGCCTTGGGTCTGTAAGCGTCGCTAAGCCGTATGACGCAGTCCTGGGTCTTCAAAATAATTTCAGCACGCTCCTTTGAGCCCAGCTCTTTGGAATCCATCATCGTTGCAAGTTG  >  NC_002947/5392163‑5392450

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