Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	816,710	A→T	5.9%	intergenic (‑63/‑334)	ybhK ← / → moaA	putative CofD superfamily transferase/molybdopterin biosynthesis protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	816,710	0	A	T	5.9%	244.1 / 5.3	68	intergenic (‑63/‑334)	ybhK/moaA	putative CofD superfamily transferase/molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGA  >  NC_000913/816564‑816858                                                                                                                                                   |                                                                                                                                                     aaGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAAc                                                                                                                                                    >  1:193003/1‑149 (MQ=255)           agaACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAatat                                                                                                                                          <  1:1665/149‑1 (MQ=255)            gaACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATAtt                                                                                                                                         >  1:506323/1‑149 (MQ=255)             aaCGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTg                                                                                                                                        >  3:277258/1‑149 (MQ=255)              aCGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGa                                                                                                                                       <  1:576930/149‑1 (MQ=255)               cgcgTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAg                                                                                                                                      <  3:89081/149‑1 (MQ=255)                       aGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACaaa                                                                                                                              <  2:121982/149‑1 (MQ=255)                        gTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAAc                                                                                                                             >  4:471009/1‑149 (MQ=255)                        gTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAAc                                                                                                                             >  1:402256/1‑149 (MQ=255)                         tCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACt                                                                                                                            >  3:66553/1‑149 (MQ=255)                            aTGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTAtt                                                                                                                         >  2:24527/1‑149 (MQ=255)                              gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTTACCTGCCAATATTGAGGATACAAACTAtttt                                                                                                                       >  1:267258/1‑149 (MQ=255)                              gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTTACCTGCCAATATTGAGGATACAAACTAtttt                                                                                                                       >  1:568601/1‑149 (MQ=255)                              gCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTAtttt                                                                                                                       >  4:522699/1‑149 (MQ=255)                               cccTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTttc                                                                                                                      <  1:359687/149‑1 (MQ=255)                                ccTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGc                                                                                                                                                <  2:303408/122‑1 (MQ=255)                                ccTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGc                                                                                                                                                >  1:303408/1‑122 (MQ=255)                                ccTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTttct                                                                                                                     >  4:104131/1‑149 (MQ=255)                                 cTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAg                                                                                                                                      <  3:296617/131‑1 (MQ=255)                                 cTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAg                                                                                                                                      >  4:296617/1‑131 (MQ=255)                                    cgccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTa                                                                                                                 >  2:616005/1‑149 (MQ=255)                                      ccgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTAcc                                                                                                               >  2:167215/1‑149 (MQ=255)                                       cgAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACct                                                                                                              >  4:126511/1‑149 (MQ=255)                                         agagCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACctct                                                                                                            >  2:235798/1‑149 (MQ=255)                                             cAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAg                                                                                                        >  1:285352/1‑149 (MQ=255)                                               aCGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                      >  2:214843/1‑149 (MQ=255)                                               aCGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGa                                                                                                      <  2:227531/149‑1 (MQ=255)                                                  acacGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGa                                                                                                   >  4:501013/1‑149 (MQ=255)                                                  acacGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGa                                                                                                   <  3:104131/149‑1 (MQ=255)                                                   cacGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGAt                                                                                                  <  1:521103/149‑1 (MQ=255)                                                      gATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                               >  1:289130/1‑149 (MQ=255)                                                      gATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                               <  2:216507/149‑1 (MQ=255)                                                      gATCAAGATCAGCCAGCGTACGATTGCGCATATAATTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCa                                                                                               <  1:510092/149‑1 (MQ=255)                                                       aTCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCAc                                                                                              >  3:522464/1‑149 (MQ=255)                                                       aTCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCAc                                                                                              >  2:78113/1‑149 (MQ=255)                                                           aGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTa                                                                                          <  2:412408/149‑1 (MQ=255)                                                                  ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTcc                                                                                   <  4:182081/149‑1 (MQ=255)                                                                  ccagcGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTcc                                                                                   <  3:178766/149‑1 (MQ=255)                                                                    agcGTACGATTGCGCATATACTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCt                                                                                 <  3:510060/149‑1 (MQ=255)                                                                      cGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtga                                                                               <  1:229940/149‑1 (MQ=255)                                                                        tACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCtgatg                                                                             >  3:305136/1‑149 (MQ=255)                                                                                   atgtaGTTTCCTGGAGTCGGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCGACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCaa                                                                  >  4:47302/4‑149 (MQ=255)                                                                                    tataGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCaaa                                                                 <  2:402256/149‑1 (MQ=255)                                                                                    tataGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCaaa                                                                 <  3:471009/149‑1 (MQ=255)                                                                                    tataGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCaaa                                                                 <  3:169382/149‑1 (MQ=255)                                                                                    tataGTGTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCaaa                                                                 <  2:340641/149‑1 (MQ=255)                                                                                               ggAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGcc                                                      <  3:126511/149‑1 (MQ=255)                                                                                                 aGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCAt                                                    <  1:271789/149‑1 (MQ=255)                                                                                                 aGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCAt                                                    >  3:515133/1‑149 (MQ=255)                                                                                                  gTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCAtt                                                   >  1:62962/1‑149 (MQ=255)                                                                                                   tCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAGACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTGTATCAAAGTAAAACCGCCAttt                                                  >  4:512038/1‑149 (MQ=255)                                                                                                      gATTATCCGCGCTATAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTcc                                               >  2:406356/1‑149 (MQ=255)                                                                                                          aTCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTa                                           <  1:24527/149‑1 (MQ=255)                                                                                                            ccGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTAtt                                         <  3:598383/149‑1 (MQ=255)                                                                                                                gcTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGt                                     <  2:506323/149‑1 (MQ=255)                                                                                                                 cTACAGTAGCGCAAATTGCGGGGAAACAGCAATTTACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTa                                    <  2:267258/149‑1 (MQ=255)                                                                                                                 cTACAGTAGCGCAAATTGCGGGGAAACAGCAATTTACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTa                                    <  2:568601/149‑1 (MQ=255)                                                                                                                 cTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTa                                    <  2:193003/149‑1 (MQ=255)                                                                                                                 cTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTa                                    <  2:301022/149‑1 (MQ=255)                                                                                                                 cTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTa                                    <  1:214843/149‑1 (MQ=255)                                                                                                                   aCAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTccc                                              >  3:158155/1‑137 (MQ=255)                                                                                                                   aCAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTccc                                              <  4:158155/137‑1 (MQ=255)                                                                                                                             aaaTTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCAc                        <  1:134092/149‑1 (MQ=255)                                                                                                                                    ggggAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATg                 >  1:564236/1‑149 (MQ=255)                                                                                                                                        aaaCAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGc             <  1:406356/149‑1 (MQ=255)                                                                                                                                             gCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCtatat        >  3:281054/1‑149 (MQ=255)                                                                                                                                               aaTTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATAc      <  1:616005/149‑1 (MQ=255)                                                                                                                                                 ttAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACAt    >  2:294288/1‑149 (MQ=255)                                                                                                                                                   aaCCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGa  <  4:515133/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AAGTGATGAGAGAACGCGTCCCAGTCCATGCCCTCCGCCGAGAGCAACGACACGATCAAGATCAGCCAGCGTACGATTGCGCATATAGTTTCCTGGAGTCAGATTATCCGCGCTACAGTAGCGCAAATTGCGGGGAAACAGCAATTAACCTGCCAATATTGAGGATACAAACTATTTTCTTCTACCTCTAAAGGACGATGCACGCTATGCCTCCCTGATGATGTATATCAAAGTAAAACCGCCATTTTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGA  >  NC_000913/816564‑816858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑