Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,907	C→G	8.9%	G150G (GGG→GGC)	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,907	0	C	G	8.9%	217.7 / 3.8	72	G150G (GGG→GGC)	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base G (2/5);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGAT  >  NC_002947/1622769‑1623053                                                                                                                                            |                                                                                                                                                   ccAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAgg                                                                                                                                          >  2:4524/1‑149 (MQ=255)   aGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTg                                                                                                                                        <  4:22574/149‑1 (MQ=255)    ggCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTgg                                                                                                                                       <  7:259876/149‑1 (MQ=255)     gCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGc                                                                                                                                      >  7:95805/1‑149 (MQ=255)             gATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTg                                                                                                                                   <  4:320972/144‑1 (MQ=255)             gATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTg                                                                                                                                   >  3:320972/1‑144 (MQ=255)              aTCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCggg                                                                                                                             >  6:264825/1‑149 (MQ=255)                 ttGTCCCCCAGCACCGGCTGGAGTTGCTAGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGATCTTCGGCATAGATGAGTACCGGCAGCGCTCCAATCTTGTTGACC‑TCGCCCGCCAACGTGGACGGGTCgggagt                                                                                                                          >  4:244659/1‑146 (MQ=255)                   gTCCCCCAGCACCGTCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGGCCACCTTGGGTGGGCTGGTcaac                                                                                                                             <  7:225011/144‑4 (MQ=255)                   gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGGCCACCTTGGGTGGGCTGGTcaac                                                                                                                             <  8:19655/144‑4 (MQ=255)                   gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGGCCACCTTGGGTGGGCTGGTcaac                                                                                                                             <  5:279796/144‑4 (MQ=255)                   gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGttt                                                                                                                          >  1:308664/1‑147 (MQ=255)                   gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGtt                                                                                                                           <  1:181254/146‑1 (MQ=255)                   gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGATCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGGCCACCTTGGGTGGGCTGGTcaacca                                                                                                                           <  5:26959/146‑6 (MQ=255)                   gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGCCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGGCCACCTTGGGTGGGCTGGTcaacca                                                                                                                           <  5:31817/146‑6 (MQ=255)                    tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTgag                                                                                                                       <  5:292412/149‑1 (MQ=255)                      ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                      <  5:264825/148‑1 (MQ=255)                      ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                     <  2:308664/149‑1 (MQ=255)                        ccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgt                                                                                                                   <  6:156725/149‑1 (MQ=255)                        ccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgt                                                                                                                   <  1:4524/149‑1 (MQ=255)                        ccAGCACCGGATGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgt                                                                                                                   >  5:157632/1‑149 (MQ=255)                                     gAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACg                                                                                                        <  5:241748/147‑1 (MQ=255)                                      aGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                       >  7:117891/1‑147 (MQ=255)                                           cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGAttt                                                                                                >  2:199107/1‑149 (MQ=255)                                           cTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGAttt                                                                                                <  8:117891/149‑1 (MQ=255)                                               aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                            <  5:126936/149‑1 (MQ=255)                                                tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCg                                                                                              >  3:195968/1‑146 (MQ=255)                                                 cacaTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGt                                                                                             <  5:331293/146‑1 (MQ=255)                                                 cacaTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGt                                                                                             >  6:331293/1‑146 (MQ=255)                                                    aTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGa                                                                                       <  6:157632/149‑1 (MQ=255)                                                        tGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGca                                                                                   >  1:288601/1‑148 (MQ=255)                                                           agtctaaCGGGGTAGGACAACAGTTCGTTGTGCTCTTCGGCATAGATTAGTACCGGCAGCGCACCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                <  8:138901/144‑1 (MQ=255)                                                            ggCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGc                                                                                >  2:77423/1‑148 (MQ=255)                                                                   cGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCa                                                                        <  1:199107/149‑1 (MQ=255)                                                                    gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGa                                                                                                   <  1:218922/121‑1 (MQ=255)                                                                    gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGa                                                                                                   >  2:218922/1‑121 (MQ=255)                                                                    gggTCGGACAACAGTTCGTTGGGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                       <  5:329745/149‑1 (MQ=255)                                                                     ggTCGGAGAACAGTTAGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTTGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAgg                                                                      <  6:259484/149‑1 (MQ=255)                                                                     ggTCGGACAAGAGTTCGTTGTTCTCTTCGGCATAGATGAGTACGGGCCGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAgg                                                                      <  3:244659/149‑1 (MQ=255)                                                                           acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCgg                                                                                                                              <  8:140575/87‑1 (MQ=255)                                                                           acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCgg                                                                                                                              >  7:140575/1‑87 (MQ=255)                                                                             aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                              >  8:75490/1‑149 (MQ=255)                                                                                  ttCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGc                                                                                    <  4:24572/122‑1 (MQ=255)                                                                                       tgtgCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTAc                                                                  >  4:25287/1‑135 (MQ=255)                                                                                       tgtgCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTAc                                                                  <  3:25287/135‑1 (MQ=255)                                                                                           ctctTCGGCATAGAGGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCAACCAAGGTGGCGTGGTCGGGTTTGAGAGGTGTGAAGCGACGTCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCATATATCGCGCa                                                <  5:282167/149‑1 (MQ=255)                                                                                            tctTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCTAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGGGTGAAGCGACGGCACATTTCGTAGATGATGCCCGCCCCTTCCAGGCTACTGGCCCAGATATCGCGCa                                                >  1:240430/1‑149 (MQ=255)                                                                                                 ggCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGttt                                            >  1:282788/1‑147 (MQ=255)                                                                                                 ggCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGttt                                            <  2:282788/147‑1 (MQ=255)                                                                                                   cATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCgcgc                                                 <  2:89492/140‑1 (MQ=255)                                                                                                   cATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCgcgc                                                 >  1:89492/1‑140 (MQ=255)                                                                                                      aGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                      >  5:31925/1‑148 (MQ=255)                                                                                                              aCCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCg                              <  6:175112/148‑1 (MQ=255)                                                                                                               ccGGCAGCGCGCCAATCTGGTTGACC‑AGCCCCCCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                            >  2:20154/1‑149 (MQ=255)                                                                                                               ccGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                            <  1:77423/149‑1 (MQ=255)                                                                                                                    aGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCg                       <  7:75490/149‑1 (MQ=255)                                                                                                                        cgcCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCAGGCa                   >  3:60089/1‑149 (MQ=255)                                                                                                                              tCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                      <  2:89384/124‑1 (MQ=255)                                                                                                                              tCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                      >  1:89384/1‑124 (MQ=255)                                                                                                                                 ggTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaa              <  6:328060/145‑1 (MQ=255)                                                                                                                                 ggTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaa              >  5:328060/1‑145 (MQ=255)                                                                                                                                 ggTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATc          >  8:206677/1‑149 (MQ=255)                                                                                                                                  gTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCt         >  4:100111/1‑149 (MQ=255)                                                                                                                                  gTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCAGGCATGGAAAATCt         <  4:60089/149‑1 (MQ=255)                                                                                                                                   ttGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg        >  2:323777/1‑149 (MQ=255)                                                                                                                                   ttGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTg        <  2:240430/149‑1 (MQ=255)                                                                                                                                    cgACC‑AGGCCACCTTGGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc       >  3:41997/2‑149 (MQ=255)                                                                                                                                      aCC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCt      <  1:323777/148‑1 (MQ=255)                                                                                                                                      aCC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTc     <  3:100111/149‑1 (MQ=255)                                                                                                                                       cc‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTc     >  3:251316/1‑148 (MQ=255)                                                                                                                                        c‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGa   >  6:254385/1‑149 (MQ=255)                                                                                                                                           gCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGAt  >  6:7792/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                                   CCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACC‑AGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGAT  >  NC_002947/1622769‑1623053

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