Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	474,812	A→T	5.5%	F58Y (TTC→TAC) ‡	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	474,812	0	A	T	5.5%	272.5 / 4.8	91	F58Y (TTC→TAC) ‡	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base A (45/41);  new base T (2/3);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

TTCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCGTCAATG  >  NC_002947/474665‑474958                                                                                                                                                     |                                                                                                                                                   ttCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGaa                                                                                                                                                   >  6:281857/1‑149 (MQ=255)      gTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGc                                                                                                                                               <  7:318283/148‑1 (MQ=255)        cgcgGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTTCTACTACTGCttc                                                                                                                                            <  3:227802/149‑2 (MQ=255)         gcggcgGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTTCTACTACTGCttc                                                                                                                                            <  3:260386/148‑2 (MQ=255)          cggcggCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTTCTACTACTGCttcgg                                                                                                                                          <  6:305621/149‑4 (MQ=255)             cggcCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCggg                                                                                                                                       <  4:343816/149‑1 (MQ=255)              ggcCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCgggg                                                                                                                                      <  6:83030/149‑1 (MQ=255)                 cTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCAGGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAGCTGGTTGTCGGGGCTg                                                                                                                                   >  6:373049/1‑149 (MQ=255)                      tCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGt                                                                                                                              >  4:162097/1‑149 (MQ=255)                        gcgGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTg                                                                                                                             <  4:171801/148‑1 (MQ=255)                        gcgGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGa                                                                                                                            <  2:88425/149‑1 (MQ=255)                                cTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGAcc                                                                                                                                <  5:357106/137‑1 (MQ=255)                                cTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGAcc                                                                                                                                >  6:357106/1‑137 (MQ=255)                                 ttCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACgggg                                                                                                                   >  8:292468/1‑149 (MQ=255)                                 ttCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACgggg                                                                                                                   >  3:262912/1‑149 (MQ=255)                                     aTGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCTTGACAAAGCCCAAGGCGTTGCCACCGGCGCCAGAGCCGAATCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGtttt                                                                                                               <  3:148370/149‑1 (MQ=255)                                     aTGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGtttt                                                                                                               >  7:2703/1‑149 (MQ=255)                                           gcggcgCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCttt                                                                                                         <  5:281857/149‑1 (MQ=255)                                            cggcgCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGGTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtg                                                                                                        <  7:313820/149‑1 (MQ=255)                                             ggcgCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtgt                                                                                                       <  2:70633/149‑1 (MQ=255)                                                gcgGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGCCAAAGCCCAAGGCGTTGCCACCGGCGCCACCGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTTACCGTGAAGGACGGCGTTTTTTCTTTGTgg                                                                                                     >  7:276128/1‑148 (MQ=255)                                                     ccAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGCCGGGGGTTTTTCTTTGTGGAACggg                                                                                               >  6:252831/1‑149 (MQ=255)                                                       aGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCGTTGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACggg                                                                                               >  1:254365/1‑148 (MQ=255)                                                         ttCCTCGACGGCCTGGGGGAAGTCCAGTTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGgcagc                                                                                           <  1:80230/149‑1 (MQ=255)                                                            cTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGcagcag                                                                                         <  1:357046/148‑1 (MQ=255)                                                            cTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAgg                                                                                        >  1:324513/1‑149 (MQ=255)                                                               gACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGcc                                                                                      >  5:384401/1‑148 (MQ=255)                                                               gACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGcc                                                                                      >  5:157755/1‑148 (MQ=255)                                                                aCGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGa                                                                                    <  3:280531/149‑1 (MQ=255)                                                                   gCCTGGGTGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagta                                                                                 <  7:292468/149‑1 (MQ=255)                                                                   gCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagta                                                                                 <  1:323979/149‑1 (MQ=255)                                                                    ccTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagtag                                                                                >  4:68754/1‑149 (MQ=255)                                                                     cTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagtagt                                                                               >  2:201257/1‑149 (MQ=255)                                                                     cTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagtagt                                                                               >  6:60558/1‑149 (MQ=255)                                                                          ggAAGTCCAGGTTGTCG‑TGGTCCGTGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTg                                                                          >  4:62016/1‑149 (MQ=255)                                                                           gAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGc                                                                         >  5:145680/1‑149 (MQ=255)                                                                               tCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGt                                                                     <  2:254365/149‑1 (MQ=255)                                                                                ccAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGtt                                                                    >  8:182791/1‑149 (MQ=255)                                                                                 cAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtgt                                                                                                       >  2:324519/1‑113 (MQ=255)                                                                                 cAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTtgt                                                                                                       <  1:324519/113‑1 (MQ=255)                                                                                 cAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGttt                                                                   >  3:216495/1‑149 (MQ=255)                                                                                  aGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGtttt                                                                  >  7:257060/1‑149 (MQ=255)                                                                                       gTCG‑TGGTCCCTGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTc                                                              <  3:62016/148‑1 (MQ=255)                                                                                        tCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAg                                                            >  7:12149/1‑149 (MQ=255)                                                                                         cG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCCGc                                                           >  5:214903/1‑149 (MQ=255)                                                                                          g‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGc                                                           <  6:183943/148‑1 (MQ=255)                                                                                              gTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTgg                                                        <  2:324513/148‑1 (MQ=255)                                                                                              gTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGa                                                       <  6:265859/149‑1 (MQ=255)                                                                                              gTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCAACGCCGAAGAAGTAGTAGAAATGCTTGTGGGGGGTGTCCGTGATGGCAGGGGTTTTTTGTTTGTGTAAGGGGCAGAAGGGCGAGTAGTTCTTGCAGGTTTTTTTGAGCTGGa                                                       >  8:161180/1‑149 (MQ=255)                                                                                                ccATGACAATGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagt                                                                                  <  5:271345/120‑1 (MQ=255)                                                                                                ccATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGgc                                                                                              <  4:111967/108‑1 (MQ=255)                                                                                                ccATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGgc                                                                                              >  3:111967/1‑108 (MQ=255)                                                                                                ccATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGagt                                                                                  >  6:271345/1‑120 (MQ=255)                                                                                                 cATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCTAGTAGTTCTAGACGGTTTTTTTCAGCTGGAcgc                                                    >  4:99460/1‑149 (MQ=255)                                                                                                  aTGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcg                                                   >  3:261014/1‑149 (MQ=255)                                                                                                   tGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcgc                                                  <  6:157755/149‑1 (MQ=255)                                                                                                    gACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGa                                                                                    >  4:345166/1‑114 (MQ=255)                                                                                                    gACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGa                                                                                    <  3:345166/114‑1 (MQ=255)                                                                                                     aCAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcgcg                                                 <  5:60558/148‑1 (MQ=255)                                                                                                      cAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGttttttt                                                               >  5:277651/1‑133 (MQ=255)                                                                                                      cAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGttttttt                                                               <  6:277651/133‑1 (MQ=255)                                                                                                      cAAAGCCCAAGGCGTTGCCACCGGCGCAACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGAGAAGGACGGCGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGAcgcgcca                                                >  8:101760/1‑146 (MQ=255)                                                                                                       aaaGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAAc                                              >  8:267358/1‑149 (MQ=255)                                                                                                       aaaGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAACGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAAc                                              >  6:331927/1‑149 (MQ=255)                                                                                                           cccAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTcac                                          <  8:12149/149‑1 (MQ=255)                                                                                                             cAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTcacca                                        >  1:218211/1‑149 (MQ=255)                                                                                                             cAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTcacca                                        >  1:333388/1‑149 (MQ=255)                                                                                                                 gCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTc                                    >  7:30248/1‑149 (MQ=255)                                                                                                                   gTTGCCACCGGCGCCACATCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTcga                                  <  7:331586/149‑1 (MQ=255)                                                                                                                    ttGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTcgac                                 <  6:37547/149‑1 (MQ=255)                                                                                                                       ccACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGACCTCACCACGTCGACGAt                              >  4:273872/1‑149 (MQ=255)                                                                                                                       ccACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGACGGTTTTTTTCAGCTGGAAGCGCGAACTCACCACGTCGACGAt                              >  1:282094/1‑149 (MQ=255)                                                                                                                        cACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATg                             >  5:177230/1‑149 (MQ=255)                                                                                                                         aCCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGt                            >  3:67078/1‑149 (MQ=255)                                                                                                                            ggCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGag                                                                                   <  3:185687/91‑1 (MQ=255)                                                                                                                            ggCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGag                                                                                   >  4:185687/1‑91 (MQ=255)                                                                                                                             gcgcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAg                        >  5:77424/1‑149 (MQ=255)                                                                                                                             gcgcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAg                        <  4:262912/149‑1 (MQ=255)                                                                                                                             gcgcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAg                        >  6:274253/1‑149 (MQ=255)                                                                                                                               gcCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGc                      <  8:257060/149‑1 (MQ=255)                                                                                                                                ccACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCATGTCGACGATGTCGAGGCg                     <  1:201257/149‑1 (MQ=255)                                                                                                                                     cgacaaGCAGTTCTACTACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTGGa                >  2:312218/5‑149 (MQ=255)                                                                                                                                     cgacaaGCAGTTCTACTACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTATCTTTGTGGAACGCGCAGCAGGCCGAGTACTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGCCGATGTCGAGGCGGTTgc                >  5:283401/5‑148 (MQ=255)                                                                                                                                       cGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATc              <  7:267358/149‑1 (MQ=255)                                                                                                                                         aaGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAg            <  3:68754/149‑1 (MQ=255)                                                                                                                                          aGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTCCGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAgg           <  4:76930/149‑1 (MQ=255)                                                                                                                                            cAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGgtc         >  2:55667/1‑149 (MQ=255)                                                                                                                                             agtagtAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGgtcg        <  4:67078/149‑1 (MQ=255)                                                                                                                                                  tagAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGCCGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCGTCaa    <  7:182791/148‑1 (MQ=255)                                                                                                                                                  tagAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCGTCAAt   <  3:99460/149‑1 (MQ=255)                                                                                                                                                   agAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCGTCAAt   >  3:261145/1‑148 (MQ=255)                                                                                                                                                    gAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCGTCAATg  <  3:162097/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TTCTGGTCGCGGCGGCCTTGCTCGCGGGGCACTTCCATGCCGGCGGCGCGGGCCAGTTCCTCGACGGCCTGGGGGAAGTCCAGGTTGTCG‑TGGTCCATGACAAAGCCCAAGGCGTTGCCACCGGCGCCACAGCCGAAGCAGTAGTAGAACTGCTTGTCGGGGCTGACCGTGAAGGACGGGGTTTTTTCTTTGTGGAACGGGCAGCAGGCCGAGTAGTTCTTGCCGGTTTTTTTCAGCTGGACGCGCGAACTCACCACGTCGACGATGTCGAGGCGGTTGATCAGGTCGTCAATG  >  NC_002947/474665‑474958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑