Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,237,936	T→A	6.0%	K156* (AAG→TAG)	PP_4615 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,237,936	0	T	A	6.0%	295.3 / 4.7	95	K156* (AAG→TAG)	PP_4615	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (44/45);  new base A (2/4);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTTGTGGT  >  NC_002947/5237794‑5238082                                                                                                                                               |                                                                                                                                                   cGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGTCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCcctc                                                                                                                                              <  5:83202/149‑1 (MQ=255)      cgcCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCa                                                                                                                                         <  1:435518/149‑1 (MQ=255)      cgcCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCAAGTAGCcctcaattg                                                                                                                                         <  5:362355/149‑6 (MQ=255)      cgcCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCAAGTAGCcctcaatt                                                                                                                                          <  8:276728/148‑5 (MQ=255)      cgcCTGCAACGGGTGATAGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCAAGTAGCcctcaattg                                                                                                                                         <  6:102706/149‑6 (MQ=255)        ccTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACt                                                                                                                                       <  1:416854/149‑1 (MQ=255)        ccTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCAAGTAGCcctcaattggc                                                                                                                                       <  1:376222/149‑8 (MQ=255)           gCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGcc                                                                                                                                    >  4:84951/1‑149 (MQ=255)            cAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGc                                                                                                                                     <  7:153337/147‑1 (MQ=255)            cAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCt                                                                                                                                   >  5:89404/1‑149 (MQ=255)              aCGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCTGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTg                                                                                                                                 >  1:194824/1‑149 (MQ=255)              aCGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTg                                                                                                                                 >  1:188444/1‑149 (MQ=255)              aCGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTCCTGGCCGCGCTGACCCGTGACGCACGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTAGCACTCCTCCACTGCCttt                                                                                                                                 >  5:74062/1‑148 (MQ=255)                 ggTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCg                                                                                                                              >  8:67286/1‑149 (MQ=255)                     aTCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTgg                                                                                                                          >  4:179516/1‑149 (MQ=255)                     aTCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTg                                                                                                                           <  5:358108/148‑1 (MQ=255)                     aTCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTg                                                                                                                           <  5:168457/148‑1 (MQ=255)                        gACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCg                                                                                                                        <  3:16014/148‑1 (MQ=255)                           aGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTc                                                                                                                      >  6:346044/1‑147 (MQ=255)                             cATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAg                                                                                                                   >  3:12923/1‑148 (MQ=255)                                ccGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCAcc                                                                                                               >  4:104933/1‑149 (MQ=255)                                 cGCCTGGGCCTTGCGATCGAACTCAGCCAGATCCCCCTGCAAGTGGTGGCCGCGCGGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCt                                                                                                              <  7:369075/149‑1 (MQ=255)                                 cGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCt                                                                                                              <  5:346044/149‑1 (MQ=255)                                   ccTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTg                                                                                                            <  3:84951/149‑1 (MQ=255)                                   ccTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTg                                                                                                            >  2:305831/1‑149 (MQ=255)                                       ggCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGt                                                                                                          <  4:177018/147‑1 (MQ=255)                                        gCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGa                                                                                                       <  4:52918/149‑1 (MQ=255)                                         ccTTGCGATCGAACTCAGCCATAACCCCCTGCTAGTGCTGGCCGCGCTGAGGCTTGTCGCCCGCCACGAGCCGCGCCATTTTGTTGGCAATTGAGGGCTTCTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAg                                                                                                      <  6:74062/149‑1 (MQ=255)                                          cTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAgg                                                                                                     <  2:428395/149‑1 (MQ=255)                                             gCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAgggg                                                                                                   >  4:328703/1‑148 (MQ=255)                                               gATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAg                                                                                                >  8:223753/1‑149 (MQ=255)                                               gATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAg                                                                                                >  2:35295/1‑149 (MQ=255)                                                   gAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGt                                                                                            >  4:382815/1‑149 (MQ=255)                                                   gAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGt                                                                                            <  5:181812/149‑1 (MQ=255)                                                   gAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGt                                                                                            <  5:188398/149‑1 (MQ=255)                                                    aaCTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGtt                                                                                                                            <  4:173737/116‑1 (MQ=255)                                                    aaCTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGtt                                                                                                                            >  3:173737/1‑116 (MQ=255)                                                    aaCTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAg                                                                                                      <  4:91510/138‑1 (MQ=255)                                                    aaCTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAg                                                                                                      >  3:91510/1‑138 (MQ=255)                                                    aaCTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTa                                                                                           <  4:253171/149‑1 (MQ=255)                                                         aGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCTCTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGaa                                                                                       <  1:109368/148‑1 (MQ=255)                                                               aaCCCCCTGCACGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAgc                                                                                <  4:230297/149‑1 (MQ=255)                                                               aaCCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGATGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAgc                                                                                <  1:305831/149‑1 (MQ=255)                                                                  ccccTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCtt                                                                                                             >  8:219301/1‑117 (MQ=255)                                                                  ccccTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCtt                                                                                                             <  7:219301/117‑1 (MQ=255)                                                                  ccccTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGca                                                                             >  5:357270/1‑149 (MQ=255)                                                                  ccccTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGca                                                                             >  6:383803/1‑149 (MQ=255)                                                                    ccTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGcac                                                                            >  6:434949/1‑148 (MQ=255)                                                                       gCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGc                                                                                                                  <  4:233875/107‑1 (MQ=255)                                                                       gCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGc                                                                                                                  >  3:233875/1‑107 (MQ=255)                                                                       gCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCaccacc                                                                        >  3:10730/1‑149 (MQ=255)                                                                        ccaGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCAccacca                                                                       <  4:370351/147‑1 (MQ=255)                                                                              cTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATc                                                                  <  4:373330/148‑1 (MQ=255)                                                                              cTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATc                                                                  <  3:328703/148‑1 (MQ=255)                                                                                 gCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCgag                                                               >  2:330483/1‑148 (MQ=255)                                                                                 gCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCgag                                                               <  4:12923/148‑1 (MQ=255)                                                                                 gCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCgag                                                               >  6:192665/1‑148 (MQ=255)                                                                                 gCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCgag                                                               >  7:232899/1‑148 (MQ=255)                                                                                 gCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCga                                                                <  8:172952/147‑1 (MQ=255)                                                                                 gCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGGCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCga                                                                <  4:10730/147‑1 (MQ=255)                                                                                   cgcgCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAAc                                                            <  6:357270/149‑1 (MQ=255)                                                                                            ccGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCaa                                                    >  5:31374/1‑148 (MQ=255)                                                                                            ccGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCaa                                                    >  3:261058/1‑148 (MQ=255)                                                                                            ccGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGAATCGAGAACGCCAGCAAc                                                   <  2:352908/149‑1 (MQ=255)                                                                                                 aCGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGAt                                              <  1:35295/149‑1 (MQ=255)                                                                                                 aCGCCCGCAACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCTCAACACCAGGATCGAGAACGCCAGCAACAGGAt                                              >  2:249439/1‑149 (MQ=255)                                                                                                  cgcccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATg                                             >  3:294491/1‑149 (MQ=255)                                                                                                   gcccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATgc                                            >  6:185918/1‑149 (MQ=255)                                                                                                    cccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCATCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATgcc                                           >  8:265048/1‑149 (MQ=255)                                                                                                     ccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGt                                                                                                          >  2:48926/1‑85 (MQ=255)                                                                                                     ccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGt                                                                                                          <  1:48926/85‑1 (MQ=255)                                                                                                     ccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCa                                                        <  8:117963/135‑1 (MQ=255)                                                                                                     ccgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCa                                                        >  7:117963/1‑135 (MQ=255)                                                                                                      cgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATc                                                                  >  3:240178/1‑124 (MQ=255)                                                                                                      cgccACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATc                                                                  <  4:240178/124‑1 (MQ=255)                                                                                                           cGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGca                                     >  5:211727/1‑148 (MQ=255)                                                                                                               ccGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGt                                <  6:89404/149‑1 (MQ=255)                                                                                                               ccGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTA‑‑GCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTa                               >  4:93194/1‑148 (MQ=255)                                                                                                                cgcgCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTa                               <  7:265048/149‑1 (MQ=255)                                                                                                                cgcgCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTa                               <  5:383803/149‑1 (MQ=255)                                                                                                                            ttgCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCtt                    <  1:204960/148‑1 (MQ=255)                                                                                                                              gCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCAt                 >  8:343661/1‑149 (MQ=255)                                                                                                                               ccAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATc                <  8:232899/149‑1 (MQ=255)                                                                                                                                  gaggAGGGCAAGTAGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCggg             >  6:149809/4‑149 (MQ=255)                                                                                                                                   ttGAGGGCTACTTGCCCTCGTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGt            <  6:31374/149‑1 (MQ=255)                                                                                                                                        ggCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCATCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCtt       >  2:339115/1‑149 (MQ=255)                                                                                                                                        ggCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCtt       >  4:293834/1‑149 (MQ=255)                                                                                                                                        ggCTACTTGCCCTCCTCCACTGCCTTGCCGTAGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCtt       >  1:266825/1‑149 (MQ=255)                                                                                                                                         gCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTtg      >  1:177489/1‑149 (MQ=255)                                                                                                                                         gCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGAGTGCCTtg      >  3:229749/1‑149 (MQ=255)                                                                                                                                           tACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTtgtg    >  6:119103/1‑149 (MQ=255)                                                                                                                                             cTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTtgtg    <  7:306520/147‑1 (MQ=255)                                                                                                                                             cTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTtgtg    <  3:179516/147‑1 (MQ=255)                                                                                                                                             cTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTtgtg    <  4:294491/147‑1 (MQ=255)                                                                                                                                             cTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTTgtggt  >  7:7906/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                   CGGTGCGCCTGCAACGGGTGATCGACAGCATCCGCCTGGGCCTTGCGATCGAACTCAGCCAGAACCCCCTGCAAGTGCTGGCCGCGCTGACCCGTGACGCCCGCCACGAGCCGCGCCATTTTGTTGCCAATTGAGGGCTACTTGCCCTCCTCCACTGCCTTGCCGTTGGCGTCCAGCACCTTGGTCGAGGGGTAGCGGTACGAAGCGTAGCGCACCACCAGGATCGAGAACGCCAGCAACAGGATGCCGCCGCACAGGTACAGCAGCCCTTCATCGGGTGCCTTGTGGT  >  NC_002947/5237794‑5238082

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