Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,757,453	A→G	5.2%	D92G (GAC→GGC)	PP_3321 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,757,453	0	A	G	5.2%	228.7 / 2.3	77	D92G (GAC→GGC)	PP_3321	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base G (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

AACTGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGAG  >  NC_002947/3757306‑3757598                                                                                                                                                    |                                                                                                                                                  aaCTGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGcgacga                                                                                                                                                   >  1:307032/1‑148 (MQ=255)    tGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGc                                                                                                                                               >  7:302678/1‑149 (MQ=255)      gcctgcctGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCg                                                                                                                                             >  7:101279/1‑149 (MQ=255)      gcctgcctGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATAAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCg                                                                                                                                             >  8:113046/1‑149 (MQ=255)               gCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGGCTGTCGTCGCagtc                                                                                                                                     <  7:168633/148‑4 (MQ=255)                      ccAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCgag                                                                                                                             <  8:101279/149‑1 (MQ=255)                      ccAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCga                                                                                                                              <  2:307032/148‑1 (MQ=255)                            cggcgCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGc                                                                                                                                         <  4:394021/131‑1 (MQ=255)                            cggcgCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGc                                                                                                                                         >  3:394021/1‑131 (MQ=255)                            cggcgCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTc                                                                                                                       <  5:258175/149‑1 (MQ=255)                               cgcTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTc                                                                                                                    <  8:302678/149‑1 (MQ=255)                                     aCCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAac                                                                                                              >  3:284813/1‑149 (MQ=255)                                      ccGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGTGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTc                                                                                                                                >  8:180230/1‑130 (MQ=255)                                      ccGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGTGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTc                                                                                                                                <  7:180230/130‑1 (MQ=255)                                              ccGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTg                                                                                                     <  4:250743/149‑1 (MQ=255)                                                 cgACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTc                                                                                                  <  4:284813/149‑1 (MQ=255)                                                    ccGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                               <  2:97983/149‑1 (MQ=255)                                                    ccGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                               <  5:121151/149‑1 (MQ=255)                                                       tGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGAc                                                                                            <  7:194251/149‑1 (MQ=255)                                                          ccACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATc                                                                                         <  2:114669/149‑1 (MQ=255)                                                           cacaCGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATc                                                                                         >  5:65647/1‑148 (MQ=255)                                                              acGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGTGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAgt                                                                                     >  5:357605/1‑149 (MQ=255)                                                                  gATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAgttgtt                                                                                 >  5:82475/1‑149 (MQ=255)                                                                   aTGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAgttgtt                                                                                 <  2:10171/148‑1 (MQ=255)                                                                   aTGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTc                                                                                <  6:357605/149‑1 (MQ=255)                                                                     gtgGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                               >  2:183994/1‑132 (MQ=255)                                                                     gtgGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                               <  1:183994/132‑1 (MQ=255)                                                                        gATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCaag                                                                             >  5:384868/1‑147 (MQ=255)                                                                          tCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAacg                                                                                                             <  8:17575/113‑1 (MQ=255)                                                                          tCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAacg                                                                                                             >  7:17575/1‑113 (MQ=255)                                                                             ggTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGcgcg                                                                      >  3:233293/1‑149 (MQ=255)                                                                               tGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGt                                                                     >  8:30532/1‑148 (MQ=255)                                                                                  tGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTc                                                                    <  2:403145/146‑1 (MQ=255)                                                                                  tGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTc                                                                    >  1:403145/1‑146 (MQ=255)                                                                                  tGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGc                                                                 >  7:279658/1‑149 (MQ=255)                                                                                  tGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTCCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGc                                                                 >  6:95199/1‑149 (MQ=255)                                                                                    gCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCg                                                               >  6:377843/1‑149 (MQ=255)                                                                                          tACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCt                                                         >  1:27591/1‑149 (MQ=255)                                                                                            cGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTgg                                                       <  7:113046/149‑1 (MQ=255)                                                                                               aCCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                               <  4:99595/106‑1 (MQ=255)                                                                                               aCCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCtgt                                                                                               >  3:99595/1‑106 (MQ=255)                                                                                               aCCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGa                                                    >  8:115950/1‑149 (MQ=255)                                                                                               aCCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGa                                                    <  5:73587/149‑1 (MQ=255)                                                                                                 cTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGAcc                                                  >  4:403861/1‑149 (MQ=255)                                                                                                 cTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGAcc                                                  >  6:231285/1‑149 (MQ=255)                                                                                                  tGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCg                                                 >  1:340357/1‑149 (MQ=255)                                                                                                   gATGTGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGt                                                >  6:409931/1‑149 (MQ=255)                                                                                                    aTGCGATCCTGTTTCCGACCGTGTCCAACCAAGTTCGTGACTGCGACGACAGCCGTCGCAACTATTTGGAGAATTTCCTCTCTAACGACGGGCTGGGCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGTGGACCGTg                                               <  2:245220/149‑1 (MQ=255)                                                                                                    aTGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTg                                               <  1:380738/149‑1 (MQ=255)                                                                                                     tGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGTGACCGTgg                                              <  7:30532/149‑1 (MQ=255)                                                                                                     tGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTgg                                              >  4:282597/1‑149 (MQ=255)                                                                                                         aTCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCAACGACAGCCATCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGt                                          >  1:303456/1‑149 (MQ=255)                                                                                                          tCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTTCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGtt                                         >  1:8710/1‑149 (MQ=255)                                                                                                          tCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGt                                                                     <  8:111999/121‑1 (MQ=255)                                                                                                          tCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGt                                                                     >  7:111999/1‑121 (MQ=255)                                                                                                            cTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCa                                                                   <  7:82066/121‑1 (MQ=255)                                                                                                            cTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCa                                                                   >  8:82066/1‑121 (MQ=255)                                                                                                               gTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCt                                                                                                 <  8:44830/88‑1 (MQ=255)                                                                                                               gTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCt                                                                                                 >  7:44830/1‑88 (MQ=255)                                                                                                               gTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGa                                      <  6:351935/147‑1 (MQ=255)                                                                                                               gTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAg                                     <  6:82475/148‑1 (MQ=255)                                                                                                                ttCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGc                                   <  2:27591/149‑1 (MQ=255)                                                                                                                   cGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCa                                 <  6:384868/148‑1 (MQ=255)                                                                                                                   cGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCa                                 <  5:409931/148‑1 (MQ=255)                                                                                                                   cGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAg                                >  1:153625/1‑149 (MQ=255)                                                                                                                    gACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAg                                >  1:174914/1‑148 (MQ=255)                                                                                                                       ccTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTaca                            >  4:90527/1‑149 (MQ=255)                                                                                                                        cTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTacac                           <  3:90527/149‑1 (MQ=255)                                                                                                                           tCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGtt                                         >  2:283062/1‑132 (MQ=255)                                                                                                                           tCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGtt                                         <  1:283062/132‑1 (MQ=255)                                                                                                                           tCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCtt                        >  1:361040/1‑149 (MQ=255)                                                                                                                             caaccaaGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTc                       <  4:233293/148‑1 (MQ=255)                                                                                                                                       cGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTga            <  2:340357/149‑1 (MQ=255)                                                                                                                                        gTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTgac           >  6:386384/1‑149 (MQ=255)                                                                                                                                           gttGCGACGGCTGTCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCaaga                                                                            >  7:180882/3‑81 (MQ=255)                                                                                                                                           gttGCGACGGCTGTCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCaaga                                                                            <  8:180882/79‑1 (MQ=255)                                                                                                                                            cTGCGACGACAGCAGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGacgacg       >  8:129445/1‑149 (MQ=255)                                                                                                                                             tGCGACGGCTGTCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACgg      >  2:378318/1‑149 (MQ=255)                                                                                                                                               cgacgaCAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGc     >  5:34519/1‑148 (MQ=255)                                                                                                                                                gacgacAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGc     >  8:299841/1‑147 (MQ=255)                                                                                                                                                 acgacAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCATTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGag  <  7:115950/149‑1 (MQ=255)                                                                                                                                                   gacAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGtt                                         >  5:213355/1‑108 (MQ=255)                                                                                                                                                   gacAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGtt                                         <  6:213355/108‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AACTGGCCTGCCTGGCAGGTGCCAAAGCGGCGCTCTACCGTTGGGCCGCGACCGTGGCCACACGTGATGTGGATCAGGTGCTGGCGTTGTACGCACCTGATGCGATCCTGGTTCCGACCCTGTCCAACCAAGTGCGTGACTGCGACGACAGCCGTCGCAACTATTTCGAGAATTTCCTCGCTAACGACGGGCTGGTCTGTGACATCCAGTTGTTCAAGAAGCGCGTCAGCCGCAAGCTGGGGACCGTGGTCGTTGGAGGCCAGTACACCTTCGTCTTTCGTGACGACGGCGAG  >  NC_002947/3757306‑3757598

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