Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	43.9%	85.8 / 12.9	59	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (29/3);  new base C (0/25);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786710‑787005                                                                                                                                                  |                                                                                                                                                        agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑Tg                                                                                                                                                   >  2:641455/1‑151 (MQ=255)  gagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑Tg                                                                                                                                                   <  1:764240/150‑1 (MQ=255)  gagGGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTGGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAGGCCGACTCACTTGCAGTCGGCTTTCACATTTTAAACGAATGACGTTTGCTTCGCTTTCCCC‑Tg                                                                                                                                                   <  2:706511/150‑1 (MQ=255)   agTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGt                                                                                                                                                 >  2:927652/1‑151 (MQ=255)          tATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAAc                                                                                                                                          >  1:548569/1‑151 (MQ=255)                      caccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGACGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCtt                                                                                                                               >  1:955055/1‑150 (MQ=255)                        ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCg                                                                                                                            >  1:98215/1‑151 (MQ=255)                          acaGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCt                                                                                                                          >  1:1126471/1‑151 (MQ=255)                           caGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                         >  2:1174775/1‑151 (MQ=255)                           caGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                         >  2:1032474/1‑151 (MQ=255)                                  gATAGGTCGCGCGGCGTCGCATCAGGCATTGTCCTCAATTTTTAGAAAAAAAAAAGCCGCCTCACTTGCAGTCGCCTTTCTCATTTGAAACGAATGCCTTTTACTTCGTTTTCCCC‑TGTTTTCAAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                  <  2:855273/151‑1 (MQ=255)                                         cggggaGCGTCGCATCAGGAAATGTGCTCCATTGTAAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGTTTTCCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                              <  1:896575/143‑1 (MQ=255)                                           cggggCGTCGCAACAGGAAAGGTCCCCCATTTTTACCAACAAAAAAGCCGCCCTACTTGCAGCCGGTTTCTTCATTTTAAACGAATGTCGTTTCCTTCGTTTTCCCC‑TGTTTTGAATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                                          <  2:384828/146‑1 (MQ=255)                                              aGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                      >  1:1054657/1‑151 (MQ=255)                                                   cGCAGCAGGCATTTTGCCCGTTTTTTAGCAAAAAAAACTCCTACCCACTCCCAATGGGTTTTCTTTTTTAAACCGAATGAC‑TTTTTTTCCCTTTTCCCGTGGTTGCCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCAt                                                                                                          <  2:955939/142‑1 (MQ=255)                                                     cATCAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTCCTTCGCTTTCCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGata                                                                                                <  1:920759/150‑1 (MQ=255)                                                      aTCAGGCAATGTCCCCCATTGTTAGAAACAAAAAACCCGACTCCCTTGCAGTCGGCTTTCTCATTTAAAACGAATGACGTTTCCTTCTCTTTCCCC‑GGTTTGGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              <  2:1106454/151‑1 (MQ=255)                                                        cAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              >  2:152544/1‑149 (MQ=255)                                                        cAGGCAATGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACAAATGACGTTTACTTCGCTTTCCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataata                                                                                             <  1:251462/150‑1 (MQ=255)                                                         aGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                           >  2:793063/1‑151 (MQ=255)                                                           gCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                         >  2:903022/1‑151 (MQ=255)                                                            cAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGtt                                                                                        >  1:642971/1‑151 (MQ=255)                                                             aaTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                       >  1:121393/1‑151 (MQ=255)                                                             aaTGTGCTCCATTGTTAACAAAAAAAAAGCCGACCCCCTGGCGGTGGGCTTTCTCATTTAAAACAAAGGCCGTTTCCTTCGCTTTCCCC‑GGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                       <  1:549226/151‑1 (MQ=255)                                                                   cTCCATTGTTAGCAAAAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAACAAATGACGTTTACTTCGCTTTCCCC‑TGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                                  <  2:758177/150‑1 (MQ=255)                                                                        tttttaaaaaaaaaaaaGCCGCCTCACTGGCAGTCGGTTTTCTTTTTTTAAACGAATGACTTTCACTTCGCTTTCCCC‑TGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                                <  2:917258/136‑1 (MQ=255)                                                                        ttGTTAGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              <  2:535648/149‑1 (MQ=255)                                                                          gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                           >  1:715874/1‑150 (MQ=255)                                                                          gTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGa                                                                            >  2:814454/1‑149 (MQ=255)                                                                             aGCAACAACAAAGACTAATCACTTGCAGTCGGCTTTCTCATTTTACACGAATCACGTTTACTTCGCTTTACTC‑TGGTTTGCAACTGCCGCCGCATTCGCTGCGATCTCGTCAGCATTACGCATATAATATCTCTTCATCGGTTTGAAAt                                                                         >  2:238895/1‑149 (MQ=255)                                                                                aaaaaaaaaCCCGACTCCGTTGCAGTCGGTTTTCTTTTTTAAAACGAATGACTTTTACTTCGCTTTACCC‑GGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTTCCCAGATAATAGCGTTTCAGCGGTTTGaa                                                                           <  1:109844/141‑1 (MQ=255)                                                                                   aaaaaaGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtcg                                                                 >  2:543634/1‑151 (MQ=255)                                                                                         gccgccTCACTTGCTGTCGGTTTTCTCATTTAAAACAAAGGACGTTTCCTTCTTTTTCCCC‑TGTTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTc                                                            <  2:35861/150‑1 (MQ=255)                                                                                                cccTTGGAGGCGCCTTTCTCATTTAAAACGAATGACGTTTACTGCGCTTCCCCC‑GTTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATaca                                                       <  2:384036/146‑1 (MQ=255)                                                                                                cccTTGCAGTCGGTTTTCTCATTTTAAACGAATGTCGTTTCCTCCCTTTTCCCC‑TGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     <  2:739002/148‑1 (MQ=255)                                                                                                 aCTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  1:660985/1‑149 (MQ=255)                                                                                                    tGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                                  <  1:148284/149‑1 (MQ=255)                                                                                                    tGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                >  1:909169/1‑151 (MQ=255)                                                                                                              tttCTCATTTTAAACGAATAACGTTTACTTCGCTTTCCCC‑GGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAg                                        <  1:1202532/149‑1 (MQ=255)                                                                                                               ttCTCATTTTAACCAAATGACTTTTCCTCCTTTTTCCCC‑GGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      <  2:1122816/150‑1 (MQ=255)                                                                                                                tctcATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCgg                                    >  2:905119/1‑151 (MQ=255)                                                                                                                      tttAAACGAATGACTTTTACTTCGTTTTCCCC‑TGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                               <  1:433438/150‑1 (MQ=255)                                                                                                                        tAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             >  1:364587/1‑150 (MQ=255)                                                                                                                         aaaCGAATGAGTTTTACTTCTCTTTCCCC‑TGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             <  1:337155/149‑1 (MQ=255)                                                                                                                         aaaCGAATGACTTTTACTTGGCTTTCCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTaa                            <  2:395589/150‑1 (MQ=255)                                                                                                                         aaaCGAATGACTTTTACTTCGCTTTCCCC‑TGTTTGACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTa                             <  2:691777/149‑1 (MQ=255)                                                                                                                          aaCGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGc                          >  2:54737/1‑151 (MQ=255)                                                                                                                            cGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                        >  1:1074733/1‑151 (MQ=255)                                                                                                                            cGAATGACGTATACTTCGCTTTACCC‑AGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAACCAGATAATAGCGTTTCAGAGGTTTGAAATTCTCGTCGAATTCATACAGCAGCGGAACGCCAGTCGATATACTAAgtt                         >  2:46699/1‑148 (MQ=255)                                                                                                                              aaTGCCCTTTACTCCGCTGTCCCC‑TGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAGACACCAGCGGCACGCCAGTCGGGATATTAAGCTCaa                      <  2:416116/151‑1 (MQ=255)                                                                                                                              aaTGACTTTTACTTCTCTTTCCCC‑TGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCaa                      <  2:648693/151‑1 (MQ=255)                                                                                                                              aaTGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGAGGCACGCCAGTCGGGATATTAAGATCa                       >  1:1168758/1‑150 (MQ=255)                                                                                                                                     tttCCTTCGCTTTCCCC‑TGGTTTACACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGaa                   <  2:437236/147‑1 (MQ=255)                                                                                                                                       ttcTTCCCTTTGCCC‑TGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct                <  1:1066984/146‑1 (MQ=255)                                                                                                                                        cctccgttttACCC‑GGTTTAGCAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc            <  2:98210/144‑1 (MQ=255)                                                                                                                                        aCTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc            >  1:857733/1‑151 (MQ=255)                                                                                                                                         cTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct           >  1:900885/1‑151 (MQ=255)                                                                                                                                               tttCCCT‑TGTTTTGAACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt      <  2:603106/150‑1 (MQ=255)                                                                                                                                               tttACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc     >  1:914628/1‑151 (MQ=255)                                                                                                                                                  aCCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  1:254721/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        AGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCC‑TGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786710‑787005

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