Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,857,372	G→T	5.2%	E433D (GAG→GAT)	PP_5133 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,857,372	0	G	T	5.2%	247.0 / 2.4	77	E433D (GAG→GAT)	PP_5133	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base T (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CCCACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTCGC  >  NC_002947/5857225‑5857515                                                                                                                                                     |                                                                                                                                                cccACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAg                                                                                                                                                 <  5:179685/148‑1 (MQ=255)  ccACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAgcg                                                                                                                                               >  4:338763/1‑149 (MQ=255)  ccACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAgcg                                                                                                                                               >  5:157080/1‑149 (MQ=255)  ccACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAgc                                                                                                                                                <  3:42681/148‑1 (MQ=255)   cACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCCGGCTGCACCGGGCCAGCCCCTCGATGCCGAgcgc                                                                                                                                              >  1:23997/1‑149 (MQ=255)    aCCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGcc                                                                                                                                             >  3:239761/1‑149 (MQ=255)    aCCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGcc                                                                                                                                             >  7:83309/1‑149 (MQ=255)    aCCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCAATCGcc                                                                                                                                             >  8:81756/1‑149 (MQ=255)     ccGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCt                                                                                                                                            >  6:249730/1‑149 (MQ=255)          aCGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATATGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGc                                                                                                                                        >  2:321428/1‑148 (MQ=255)           cGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGa                                                                                                                                      <  3:338763/149‑1 (MQ=255)              cTGGCCATCGCCTGCGACCCACGCCGCTCGCCCGGACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGAtt                                                                                                                                    >  1:238840/1‑149 (MQ=255)                 gCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCa                                                                                                                                <  5:309659/149‑1 (MQ=255)                 gCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCAATCGCCCAGGCGCTCGGc                                                                                                                                 <  7:81756/148‑1 (MQ=255)                  ccATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCAc                                                                                                                               <  2:211625/149‑1 (MQ=255)                    aTCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTgg                                                                                                                                  >  7:122777/1‑144 (MQ=255)                    aTCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTgg                                                                                                                                  <  8:122777/144‑1 (MQ=255)                    aTCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACg                                                                                                                              >  5:366173/1‑148 (MQ=255)                       gCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAg                                                                                                                           <  8:83309/148‑1 (MQ=255)                        ccTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGAAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGcc                                                                                                                         >  1:115527/1‑149 (MQ=255)                         cTGCGACCCACGCCGCTCGGGC‑GAGCGCGGCTGCTTGGGGCTGCTCGCCGATCTGGCGCGCAATGCCGGCGCAAGCCGTATCTGGCTGTTGCAGGCTGCACGGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGTTTGGCACGAAGccc                                                                                                                        <  6:366173/149‑1 (MQ=255)                         cTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGccc                                                                                                                        <  6:157080/149‑1 (MQ=255)                         cTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGccc                                                                                                                        <  4:178297/149‑1 (MQ=255)                            cGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCtt                                                                                                                      <  3:122582/148‑1 (MQ=255)                            cGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCtt                                                                                                                      <  2:42120/148‑1 (MQ=255)                             gACCCACGCCGCGCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGa                                                                                                                    <  1:146543/149‑1 (MQ=255)                                    gccgcTCGCGC‑GACCGCGGCAGCTTGGCGCTGCTCGACGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCtt                                                                                                             >  1:61800/1‑149 (MQ=255)                                    gccgcTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCtt                                                                                                             >  5:17817/1‑149 (MQ=255)                                     ccgcTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTg                                                                                                            <  7:267329/149‑1 (MQ=255)                                        cTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCAc                                                                                                                               <  3:8881/127‑1 (MQ=255)                                        cTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCAc                                                                                                                               >  4:8881/1‑127 (MQ=255)                                         tCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCt                                                                                                        >  1:141620/1‑149 (MQ=255)                                         tCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCAATCGCCCAGGCGCTCGGCa                                                                                                                                >  6:120821/1‑125 (MQ=255)                                         tCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCAATCGCCCAGGCGCTCGGCa                                                                                                                                <  5:120821/125‑1 (MQ=255)                                           gCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTg                                                                                                       >  5:327831/1‑148 (MQ=255)                                            ccc‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGc                                                                                                      >  4:168971/1‑148 (MQ=255)                                                   cgcgGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGc                                                                                               <  2:141620/149‑1 (MQ=255)                                                   cgcgGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGc                                                                                               <  2:238840/149‑1 (MQ=255)                                                     cgGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCg                                                                                             >  5:185506/1‑149 (MQ=255)                                                      ggCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCaa                                                                                                    >  4:102444/1‑141 (MQ=255)                                                      ggCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCaa                                                                                                    <  3:102444/141‑1 (MQ=255)                                                                     ctCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGa                                                                              <  1:319454/148‑1 (MQ=255)                                                                       cGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGCCCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACct                                                                           >  2:251049/1‑149 (MQ=255)                                                                        gCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGAcc                                                                            <  3:168971/147‑1 (MQ=255)                                                                        gCCGAACTGGCGCGCAATGCCGGCACAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTg                                                                                                            <  8:100463/115‑1 (MQ=255)                                                                        gCCGAACTGGCGCGCAATGCCGGCACAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTg                                                                                                            >  7:100463/1‑115 (MQ=255)                                                                         ccGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTTCACAGGTCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACctgg                                                                         <  5:357942/149‑1 (MQ=255)                                                                           gAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACctggct                                                                       >  3:20927/1‑149 (MQ=255)                                                                                 gcgcgcAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGCCCGCCTTGGCCTGCAACATGCCGACACCTCACAGCTGACCTGGCTGGAGCa                                                                 >  6:292015/1‑149 (MQ=255)                                                                                   gcgcAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATg                                                               <  6:224727/149‑1 (MQ=255)                                                                                      cAATACCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGcc                                                            <  3:113790/149‑1 (MQ=255)                                                                                            cGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACt                                                      >  1:364913/1‑149 (MQ=255)                                                                                             ggCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTg                                                     <  8:227354/149‑1 (MQ=255)                                                                                                        tATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGaaa                                          <  6:17817/149‑1 (MQ=255)                                                                                                            tGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTgg                                      >  5:457/1‑149 (MQ=255)                                                                                                              gctgctGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTgg                                      <  2:115527/147‑1 (MQ=255)                                                                                                                tgctgcTGGCGGCACCGGGCCAGCCCCTCGATGGCGTGCGCCTGGGCGATTGGCAGGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGGCCTGGCTGGAGCATGGc                                                             <  6:185469/122‑1 (MQ=255)                                                                                                                tgctgcAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGGCCTGGCTGGAGCATGGc                                                             >  5:185469/1‑122 (MQ=255)                                                                                                                  ctgcAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGt                                  <  4:239761/147‑1 (MQ=255)                                                                                                                     cAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCg                             <  5:130186/149‑1 (MQ=255)                                                                                                                       ggCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGc                           >  5:200864/1‑149 (MQ=255)                                                                                                                        gCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGcc                          <  2:61800/149‑1 (MQ=255)                                                                                                                        gCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGcc                          <  7:7011/149‑1 (MQ=255)                                                                                                                          tGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCgg                            <  6:19130/145‑1 (MQ=255)                                                                                                                          tGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCgg                            >  5:19130/1‑145 (MQ=255)                                                                                                                          tGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCa                         >  6:368889/1‑148 (MQ=255)                                                                                                                            cACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATAc                      >  3:63935/1‑149 (MQ=255)                                                                                                                               cGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCgg                            <  5:261003/140‑1 (MQ=255)                                                                                                                               cGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCgg                            >  6:261003/1‑140 (MQ=255)                                                                                                                               cGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCaa                   >  3:254855/1‑149 (MQ=255)                                                                                                                                gggCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCACACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAAt                  <  5:249730/149‑1 (MQ=255)                                                                                                                                   ccAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGtt               <  2:364913/149‑1 (MQ=255)                                                                                                                                     aGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTgg             <  4:63935/149‑1 (MQ=255)                                                                                                                                        cccTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCaa          <  6:457/149‑1 (MQ=255)                                                                                                                                          cTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCaaaa        <  6:200864/149‑1 (MQ=255)                                                                                                                                            cGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAAcc      >  2:123171/1‑149 (MQ=255)                                                                                                                                              aTGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTc    >  7:334328/1‑149 (MQ=255)                                                                                                                                              aTGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTc    >  5:285711/1‑149 (MQ=255)                                                                                                                                               tGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTCg   >  7:113641/1‑149 (MQ=255)                                                                                                                                                 ccGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTCGc  <  4:20927/148‑1 (MQ=255)                                                                                                                                                   gAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTCg   <  8:102186/145‑1 (MQ=255)                                                                                                                                                   gAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTCg   >  7:102186/1‑145 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCCACCGGCACGCCTGGCCATCGCCTGCGACCCACGCCGCTCGCCC‑GACCGCGGCAGCTTGGCGCTGCTCGCCGAACTGGCGCGCAATGCCGGCGCAACCCGTATCTGGCTGCTGCAGGCTGCACCGGGCCAGCCCCTCGATGCCGAGCGCCTGGGCGATTGGCACGAAGCCCTTGACCGCCTTGGCCTGCAACATGCCGACACCTCACCGCTGACCTGGCTGGAGCATGGCCATGACTGAGCCACTGAAACTGGCCGTGGTCGGCCATACCAATGTTGGCAAAACCTCGC  >  NC_002947/5857225‑5857515

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