Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,227	0	A	G	44.3%	94.7 / 17.3	61	V7A (GTC→GCC)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base G (27/0);  total (42/19)
Fisher's exact test for biased strand distribution p-value = 6.79e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTG  >  NC_000913/1312083‑1312377                                                                                                                                                 |                                                                                                                                                        aCTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAAAAAATCGCCCTGAGGGGCg                                                                                                                                                       >  2:252909/1‑147 (MQ=255)    aCGCGACCGTGGGTAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCCGAGGGGCggt                                                                                                                                                     >  1:235856/1‑144 (MQ=255)        gACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                                                                                                                            <  1:715578/151‑1 (MQ=255)         aCCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑a                                                                                                                                          <  2:423144/151‑1 (MQ=255)          ccGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCGAAAGAAAGGAAAAAAAAAACGCCCTGCGGGGGGGTAtgtg                                                                                                                                                >  1:985354/1‑145 (MQ=255)                gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATgg                                                                                                                                   <  2:161956/151‑1 (MQ=255)                 cAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGt                                                                                                                                  <  2:43425/151‑1 (MQ=255)                  aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTa                                                                                                                                 <  1:1021378/151‑1 (MQ=255)                  aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGGACGTAAAAAAAGAAAGCCCCGGGGGGCGGTATGTGTTTT‑AGACATgg                                                                                                                                   >  2:823503/1‑149 (MQ=255)                      tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGGATGTAAAAAAAGATTGCCCCGCGGGGCGGTTTGTGTTGT‑TGACAAGGTaaaa                                                                                                                              >  1:99670/1‑150 (MQ=255)                      tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCTGAGGGGCGTTATGTGTTGT‑TGACATGGGaaaa                                                                                                                              >  2:1003911/1‑150 (MQ=255)                      tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCAAAATACGGAAAAAAAAATCGCCCTGAGGGGGGTTATGTGTTGT‑AGACAAGGGaaaa                                                                                                                              >  2:30613/1‑150 (MQ=255)                        ttttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCCGAGGGGCGTTTTGTGTTGT‑TGAAAAGGTaaaaac                                                                                                                            >  2:159321/1‑148 (MQ=255)                              cAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTaa                                                                                                                     <  2:605536/151‑1 (MQ=255)                                  aTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGAAAGTAAAAAAAGAACGCCCCGAGGGGCGTTATGGGTTGT‑TGACAAGGTAAAAACAACTTaaaa                                                                                                                   >  1:444304/1‑149 (MQ=255)                                   tAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATAAGTAAAAAAAAATCGCCCCGAGGGGCGGTATGTGTTGT‑AGACATGGTAAAAACGACTTaaaa                                                                                                                   >  2:528152/1‑148 (MQ=255)                                    aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAATACGTAAAAAAAAAACGCCCCGAGGGGCGGTATGTGGTGT‑TGAAATGGGAAAAACCACTTaaaa                                                                                                                   >  1:502802/1‑147 (MQ=255)                                    aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGaa                                                                                                               <  2:219165/151‑1 (MQ=255)                                    aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAGAAAGCCCTGAGGGGGGGTATGTGTTGT‑AGACATGGGAAAAACGACTTaaaaa                                                                                                                  >  1:295968/1‑147 (MQ=255)                                      cgccgcCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGACAATGCTAAAAGACGAAAAAAAAAAACGCCCTGAGGGGGGGTTTGTGGTTT‑AGACAAGGTAAAAACAAATTAAAAAGaa                                                                                                               >  1:386397/1‑149 (MQ=255)                                       gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGCAAAAAACGTAAAAAAAAATAGCCCCTAGGGGCGGTTTGTGTTGT‑GGACATGGGAAAAAACAATTAAAAAGaaaa                                                                                                             >  1:905334/1‑150 (MQ=255)                                        ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGGAAAAAAAGATTGGCCCCGGGGAGGCTCTGTTTTGT‑AGACATGGGAAAAACGACTTaaaaa                                                                                                                  >  1:458750/1‑143 (MQ=255)                                            cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGCAAAAAAAGTAAAAAAAAAACGCCCTTAGGGGGGGTTTGTGTTGT‑AGAAATGGTAAAAACCAATTTAAAAGAAAATaac                                                                                                         >  1:55067/1‑149 (MQ=255)                                            cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGAAAGGAAAAAAAAAACGCCCTGGGGGGGGGTATGTGTTGT‑AGACATGGGAAAAAACACCTAAAAAAAAAAAcaca                                                                                                        >  1:984690/1‑150 (MQ=255)                                             aaTATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCCGGGGAAAAAAAGTAAAAAAAAATTGCCCTGAGGGGCGGTTTTTGTTGTAAAACAGGGAAAAACGCCCTTAAAAGAaaaaaaaag                                                                                                        >  2:495479/1‑145 (MQ=255)                                                   cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTAAAAAACGTAAAAAAAGAACGCCCCGCGGGGCGGTATGTGTTTT‑AGAAAAGGGAAAAACGACCTTAAAAGAAAAACACAAGcgcg                                                                                                  >  1:803339/1‑149 (MQ=255)                                                          gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATAAGTAAAACAAAAAAGCCCCGAGGGGCGGTATTTGTTGT‑AGAAATGGGAAAAACCACTTAAAAAGAAAAAAACAGGcgcgggggta                                                                                            >  2:173640/1‑142 (MQ=255)                                                              ttAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGCCATCGCCAAAAAAAGTAAAAAAAAAAAGCCCTCGGGGGGGGTTTTTGTTGT‑AGAAAAGGGAAAAAACACTTTAAAAGAAAAacaccagggccgtggtac                                                                                           >  1:989304/1‑131 (MQ=255)                                                                 accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAAAAGATCGCCCCGAGGGGGGGTATGTGTTGT‑TGACATGGTAAAAACGACTTAAAAAGAAAAACACAGGCGCGGTGCTACTATGAtttt                                                                                  >  2:396644/1‑151 (MQ=255)                                                                   caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCTGGGGGGCGGTATGTGTTGT‑AGACATGGTAAAAACGACTTTAAAAGAAAAACACAGGCGCGATGCTACTATGAtttt                                                                                  >  2:455784/1‑149 (MQ=255)                                                                    aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGACATCGCCAAAGAACGAAAAAAAAAAATTCCCTGCGGGGGGGTTTGTGTTTT‑AGAAAAGGGGAAAACCACTTTAAAAAAAAATTACAAGGGCGATGCTACTATTAtttt                                                                                  >  1:565519/1‑148 (MQ=255)                                                                      ccgccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTga                                                                             <  2:190255/148‑1 (MQ=255)                                                                       caccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTgag                                                                            <  1:204800/151‑1 (MQ=255)                                                                           aaaGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATAAGTAAAAAAAAATCGCCCTGTGGGGGGGTATGTGTTGT‑AGACATGGGAAAACCGACTTAAAAAGAAAAACACAAGCGCGATGCTACTATGATTTTTCTgagggg                                                                         >  2:173006/1‑148 (MQ=255)                                                                           aaaGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAACAAAATCGCCCTGAGGGGCGGTATGTGTTGT‑AGAAATGGTAAAAACGACTTAAAAAGAAAATCACAGGCGCGATGCTACTATGATTTTTCTgagg                                                                           >  2:617024/1‑148 (MQ=255)                                                                               aTGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGGAAAAAAAAAACGGCCTGAGGGGGGGTATGTTTTTT‑AGACATGGGAAAAACGACTTAAAAAGAAAAATACAAGGGCGATGGTTCTATTATTTTTCTTCGGAGGACa                                                                     >  2:763937/1‑150 (MQ=255)                                                                                       attggcattggTATCGGCGGGCATTGCTAAAATACGGAAAAAAAAAAATCCCCTAGGGGGGTTATGTGTTGT‑TGACATGGTAAAAAACACTTAAAAAGAAAAAAACACGCGGGCTGCAACTATGATTTTTTTtaggaggaaaaa                                                                   >  1:148717/1‑138 (MQ=255)                                                                                          ggcattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAg                                                         <  1:864175/151‑1 (MQ=255)                                                                                                       gcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGaa                                            <  2:481971/151‑1 (MQ=255)                                                                                                       gcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGaa                                            <  1:656275/151‑1 (MQ=255)                                                                                                        cggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGa                                             >  2:372435/1‑149 (MQ=255)                                                                                                              tGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTgg                                     >  2:586691/1‑151 (MQ=255)                                                                                                              tGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTgg                                     >  2:529846/1‑151 (MQ=255)                                                                                                               ggCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTgg                                     <  2:853366/150‑1 (MQ=255)                                                                                                                   aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGt                                 >  1:834476/1‑150 (MQ=255)                                                                                                                    aaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCt                               >  1:636018/1‑151 (MQ=255)                                                                                                                       tACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGt                            >  1:914306/1‑151 (MQ=255)                                                                                                                       tACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGt                            >  2:409485/1‑151 (MQ=255)                                                                                                                        aCGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt                           >  2:701334/1‑151 (MQ=255)                                                                                                                           tAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAgg                        >  2:556020/1‑151 (MQ=255)                                                                                                                             aaaCAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAt                      <  1:675547/151‑1 (MQ=255)                                                                                                                             aaaCAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAt                      <  2:444348/151‑1 (MQ=255)                                                                                                                              aaCAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAtt                     <  1:71877/151‑1 (MQ=255)                                                                                                                               aCAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAttt                    >  1:78636/1‑151 (MQ=255)                                                                                                                                  aGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAGGGACTGGCAGTCTGGTTAGGATTTATc                 <  1:380674/151‑1 (MQ=255)                                                                                                                                   gATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCACCGAAGGACTGGCAGTCTGGTTAGGATTTATCt                >  1:291850/1‑151 (MQ=255)                                                                                                                                     tCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTc              >  1:624761/1‑151 (MQ=255)                                                                                                                                     tCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTc              >  1:986581/1‑151 (MQ=255)                                                                                                                                         ccTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGc          >  1:422917/1‑151 (MQ=255)                                                                                                                                              gggACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCAt     <  1:293588/151‑1 (MQ=255)                                                                                                                                              gggACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCAt     <  2:98651/151‑1 (MQ=255)                                                                                                                                                 aCGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCAtgtg  >  1:759136/1‑151 (MQ=255)                                                                                                                                                 aCGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCAtgtg  >  2:250274/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                        ACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGT‑AGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTG  >  NC_000913/1312083‑1312377

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