Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,787,070	T→G	7.7%	intergenic (+37/‑145)	tsf → / → pyrH	elongation factor Ts/uridylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,787,070	0	T	G	7.7%	198.4 / 7.5	65	intergenic (+37/‑145)	tsf/pyrH	elongation factor Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (3/2);  total (33/32)
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.91e-01

AGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCG  >  NC_002947/1786928‑1787212                                                                                                                                                |                                                                                                                                               aGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcg                                                                                                                                          >  3:268906/1‑149 (MQ=255) aGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcg                                                                                                                                          >  4:19968/1‑149 (MQ=255)    aaTCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcgcg                                                                                                                                        >  3:169453/1‑148 (MQ=255)     aTCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCa                                                                                                                                      <  7:185671/149‑1 (MQ=255)             cTTCACCTACTTCAACGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCtt                                                                                                                               <  7:308555/148‑1 (MQ=255)             cTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCctc                                                                                                                                 >  4:355454/1‑146 (MQ=255)             cTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCctc                                                                                                                                 <  3:355454/146‑1 (MQ=255)             cTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                                                                                                                              <  8:108929/149‑1 (MQ=255)             cTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                                                                                                                              <  2:65645/149‑1 (MQ=255)             cTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                                                                                                                              <  4:268906/149‑1 (MQ=255)                  ccTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCctct                                                                                                                                <  8:160102/142‑1 (MQ=255)                  ccTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCctct                                                                                                                                >  7:160102/1‑142 (MQ=255)                   cTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTATGCCGATCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCa                                                                                                                          >  1:106652/1‑148 (MQ=255)                   cTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCaaa                                                                                                                        <  3:340949/149‑1 (MQ=255)                      cTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACg                                                                                                                     >  7:318268/1‑149 (MQ=255)                       ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACgg                                                                                                                    >  8:187548/1‑149 (MQ=255)                       ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACgg                                                                                                                    >  4:172751/1‑149 (MQ=255)                                ggCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGtt                                                                                                           >  7:281924/1‑149 (MQ=255)                                 gCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGttt                                                                                                          >  5:296201/1‑149 (MQ=255)                                  cGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTta                                                                                                         <  3:19968/149‑1 (MQ=255)                                   gAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTtat                                                                                                        <  1:339020/149‑1 (MQ=255)                                   gAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTta                                                                                                         >  1:214294/1‑148 (MQ=255)                                   gAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTta                                                                                                         >  1:97685/1‑148 (MQ=255)                                    aaGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTtata                                                                                                       <  2:106652/149‑1 (MQ=255)                                            gagaAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTcc                                                                                                <  2:214294/148‑1 (MQ=255)                                              gaAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCctc                                                                                             <  8:281924/149‑1 (MQ=255)                                              gaAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCctc                                                                                             >  2:1258/1‑149 (MQ=255)                                              gaAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCctc                                                                                             >  2:6923/1‑149 (MQ=255)                                                  ccAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCg                                                                                         >  1:93908/1‑149 (MQ=255)                                                    aGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCt                                                                                       >  5:328311/1‑149 (MQ=255)                                                              ttCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAg                                                                             >  5:121465/1‑149 (MQ=255)                                                                 gctgctGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCg                                                                           <  1:189136/148‑1 (MQ=255)                                                                   tgctgAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTc                                                                                          >  3:199319/1‑131 (MQ=255)                                                                   tgctgAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTc                                                                                          <  4:199319/131‑1 (MQ=255)                                                                   tgctgAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTg                                                                        >  4:315018/1‑149 (MQ=255)                                                                   tgctgAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTg                                                                        <  8:318268/149‑1 (MQ=255)                                                                                 gCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGtt                                                           <  6:328311/149‑1 (MQ=255)                                                                                 gCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGt                                                            >  5:80102/1‑148 (MQ=255)                                                                                       gTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGtt                                                           <  2:253923/143‑1 (MQ=255)                                                                                       gTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGtt                                                           >  1:253923/1‑143 (MQ=255)                                                                                          gctgctGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTg                                                  >  2:246468/1‑149 (MQ=255)                                                                                            tgctgcCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGc                                                 >  6:50280/1‑148 (MQ=255)                                                                                                  cAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGATGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGc                                          <  8:16867/149‑1 (MQ=255)                                                                                                   aaGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCctct                                                                                                                                >  1:296229/1‑62 (MQ=17)                                                                                                   aaGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCGCGCGTGAGCGGGCAGCctc                                                                                                                                 <  2:296229/61‑1 (MQ=17)                                                                                                    aGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGccc                                        <  2:97685/149‑1 (MQ=255)                                                                                                    aGCAGTAAGACAGCCCCGTCTGTCGCACCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGccc                                        >  3:226518/1‑149 (MQ=255)                                                                                                     gCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCg                                       <  5:50280/149‑1 (MQ=255)                                                                                                           aGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGaa                                 >  3:295628/1‑149 (MQ=255)                                                                                                             aCAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAAtt                               >  1:192701/1‑149 (MQ=255)                                                                                                                   ccGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTaa                         <  1:1258/149‑1 (MQ=255)                                                                                                                   ccGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTaa                         <  1:6923/149‑1 (MQ=255)                                                                                                                     gTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTaaaa                       <  6:121465/149‑1 (MQ=255)                                                                                                                     gTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTaaaa                       <  7:187548/149‑1 (MQ=255)                                                                                                                            gCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTAcgcc                <  4:169453/149‑1 (MQ=255)                                                                                                                                cAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAg            >  5:119824/1‑149 (MQ=255)                                                                                                                                  aaGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAg            >  3:361708/1‑147 (MQ=255)                                                                                                                                  aaGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGga          <  4:295628/149‑1 (MQ=255)                                                                                                                                   agagGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgag         <  6:119824/149‑1 (MQ=255)                                                                                                                                    gagGCTGCGCGCGTGAGCGGGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGga          >  7:137153/1‑147 (MQ=255)                                                                                                                                    gagGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgaga        <  2:192701/149‑1 (MQ=255)                                                                                                                                            ccGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgaga        <  7:64871/141‑1 (MQ=255)                                                                                                                                            ccGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGgaga        >  8:64871/1‑141 (MQ=255)                                                                                                                                              gCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCg  <  6:83982/145‑1 (MQ=255)                                                                                                                                              gCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCg  >  5:83982/1‑145 (MQ=255)                                                                                                                                                tCACGCGCGCAGCCTCTTTGTCAAAACTGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCACACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCg  <  4:25337/143‑1 (MQ=255)                                                                                                                                                tCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCACACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCg  >  3:25337/1‑143 (MQ=255)                                                                                                                                                |                                                                                                                                               AGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGT‑TGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTATACGACTCCCTCTTCGCTGGCACCGAAGCGGTGCCACGCTACAGTTAGCAGGCTGCAAACAGCCCGCACGAATTTTCTAAAAGTACGCCGCAGGAGAGACTCG  >  NC_002947/1786928‑1787212

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