Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,733	0	T	G	42.5%	113.4 / 15.6	73	G20G (GGT→GGG)	yegD	Hsp70 chaperone family protein
Reads supporting (aligned to +/- strand):  ref base T (10/32);  new base G (31/0);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 3.90e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGG  >  NC_000913/2147586‑2147880                                                                                                                                                     |                                                                                                                                                    gCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGaa                                                                                                                                                   >  2:615509/1‑150 (MQ=255)   tGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TTCGTGACGGGaaa                                                                                                                                                  >  1:478031/1‑149 (MQ=255)    gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAAcc                                                                                                                                                <  1:1097145/150‑1 (MQ=255)    gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAAcc                                                                                                                                                <  2:856322/150‑1 (MQ=255)    gcgcTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGaaaac                                                                                                                                                >  1:1118222/1‑148 (MQ=255)      gctttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGGCGGGaaa                                                                                                                                                  >  1:682919/1‑146 (MQ=255)        tttgctttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGACA‑TGCGTGACGGGaaaaaccatt                                                                                                                                           >  2:256989/1‑144 (MQ=255)              ttCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGGCGGGAAAACGCATTTTCt                                                                                                                                       >  2:341762/1‑149 (MQ=255)                   atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAAt                                                                                                                                 <  1:139176/150‑1 (MQ=255)                   atTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAAt                                                                                                                                 <  1:651371/150‑1 (MQ=255)                        ggTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGaaaa                                                                                                                           <  1:563113/151‑1 (MQ=255)                         gTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAAc                                                                                                                          <  1:320349/151‑1 (MQ=255)                          ttGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAACCGCCTTTGCTAAAAAAGGaaaa                                                                                                                           >  1:192464/1‑149 (MQ=255)                           tGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACg                                                                                                                         <  2:983633/150‑1 (MQ=255)                               gcgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGc                                                                                                                    <  2:173683/151‑1 (MQ=255)                                cgACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGTGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCa                                                                                                                   <  2:273426/151‑1 (MQ=255)                                  aCGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCATGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACg                                                                                                                 <  2:1126543/151‑1 (MQ=255)                                    ggCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACgct                                                                                                               <  2:733958/151‑1 (MQ=255)                                       aaCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TTCCTTACGGGAAACCCCATTTTCTTAAAAATGAAAAcccacgcacccc                                                                                                               >  1:960194/1‑137 (MQ=255)                                       aaCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGGGGGACGGGAAAACGCATTTTTTTAAAAAGGAAAAAAACAGCACgccg                                                                                                              >  1:1037304/1‑147 (MQ=255)                                       aaCATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGGAAACCGCCTTTGCTAAAAATAGAAAAAGACAACACcctgc                                                                                                             >  1:838212/1‑150 (MQ=255)                                          aTTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGGAAAACCCCTTTGTTTAAAAAAGAAAAAAACAACACGCCGCTGcc                                                                                                         >  2:618613/1‑151 (MQ=255)                                            ttCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCt                                                                                                        <  1:569266/150‑1 (MQ=255)                                             tCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCtt                                                                                                       <  2:951959/150‑1 (MQ=255)                                                 tAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCA‑TGCGGGACGGTAAACAGCACTTTCTAAAAATGGAAAAAGACAGCAAGCCGCTGcc                                                                                                         >  2:75791/1‑144 (MQ=255)                                                  aaaGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAAt                                                                                                   <  2:326472/149‑1 (MQ=255)                                                     gTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAAACGCATTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCttt                                                                                              >  1:51734/1‑151 (MQ=255)                                                       cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCC‑TGCGTGACGGGAAAACGCCTTTTCTAAAAATGGAAAAAGACAACACGCTGCTGCCTTTAATGCttt                                                                                              >  1:805059/1‑149 (MQ=255)                                                       cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGGAAAACCCATTTTGTAAAAATTGAAAAAGACAGCCCCCCGCCGCCTTCCATGCttt                                                                                              >  2:633693/1‑149 (MQ=255)                                                       cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTgc                                                                                            <  2:452719/151‑1 (MQ=255)                                                       cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTg                                                                                             <  2:1019686/150‑1 (MQ=255)                                                       cACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGATTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TGGGGGGCGGTAAAACGCCGTTTTTTAAAAAGGAAAAAGACAACACGCCGCTGGCCTCAAACCttt                                                                                              >  2:958804/1‑149 (MQ=255)                                                        aCGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGAAAACGCCTTTGCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCCATGCttt                                                                                              >  1:849996/1‑148 (MQ=255)                                                         cGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTgcgc                                                                                          <  1:1132753/151‑1 (MQ=255)                                                           ccccTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TTCGTGACGGGAAACCGCAATTGCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTTAATGCttt                                                                                              >  1:402373/1‑145 (MQ=255)                                                            cccTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGGGACGGGAAAACGCATTTGCTTAAAAAGGGAAACGACAGCACGCTGCTGCCCTCAATGCTTTTCGCGcc                                                                                       >  2:724680/1‑151 (MQ=255)                                                                   tcACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGACT‑TTGCTGATGGGAAAACCCCTTTGTTACAAATAGAAAACGGCAACCCGCCGCTTCCCTCCACGCTTTGGGCGccccaccg                                                                                 >  1:286315/1‑144 (MQ=255)                                                                    cACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTg                                                                               <  1:287633/151‑1 (MQ=255)                                                                       cGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAg                                                                            <  1:419928/151‑1 (MQ=255)                                                                           aTGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TTCGTTACGGGAAAACGCCCTTGCTTAAAAAGGAAAACGACAGCACGCTGCTGCCCTCAATGCTTTTCGCGCCCACGCGGGaaaagg                                                                         >  2:863226/1‑146 (MQ=255)                                                                           aTGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAAACGCATTTGCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAAAGCTTTGCGCGCCAACGCGGGAAGCgg                                                                         >  1:638295/1‑150 (MQ=255)                                                                             ggggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTaa                                                                      <  2:403854/151‑1 (MQ=255)                                                                             ggggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAACGACAGCACGCTTCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCgg                                                                         >  1:763185/1‑148 (MQ=255)                                                                              gggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TGCGGGAAGGTAAAACCCATTTTGTAAAAATTGAAAAACACAAGAAGCGGCTGGCTTCAAATCTTTCTGGGCCAAACCGGGAAGAgg                                                                         >  2:934915/1‑147 (MQ=255)                                                                              gggATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATTGCG‑GAAGGGGAAACACATTTTGTAAAAATGGAAAAACAAAGGACGCTGCTGGCTTCAATGCTTTGGGCGCCAACGCGTGAAGCgg                                                                         >  1:33242/1‑147 (MQ=255)                                                                                gATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TGTGTGGAGGTAAACCGCATTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGGaa                                                                      >  2:461049/1‑148 (MQ=255)                                                                                gATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGGAAACCCCCTTTGCTAAAAAAGGAAAACGACAGCACCCTGCCGCCTTCAATGCTTTTCGCGCCCACGCGCGAAGCGGGAAGc                                                                    >  2:938259/1‑150 (MQ=255)                                                                                gATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAACGGGAAGCg                                                                   >  1:160042/1‑151 (MQ=255)                                                                                gATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAAACGCATTTGCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAACGGGaa                                                                      >  1:920820/1‑148 (MQ=255)                                                                                   tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTTA‑TGCGTTTCGGGAAAACCCCTTTTCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCCACGCGGGAAGCGGTAAGCGaa                                                                 >  2:696365/1‑150 (MQ=255)                                                                                   tattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGaa                                                                 <  2:448074/150‑1 (MQ=255)                                                                                    attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGGGGGACGGGAAAACGCATTTGCTTAAAATGGAAAAAGACAACACGCTGCTGCCTTCAATGCTTTTCGCGCCAACCCGTGAAGCGGGAAGCGGATg                                                               >  1:1025891/1‑151 (MQ=255)                                                                                        cGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGAAGGTAAAACCCATTTTCTAAAAAAGGGAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGGGCCAACCGGTGAAACGGGAAGCGGAAGGCTg                                                           >  2:138118/1‑151 (MQ=255)                                                                                              tATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACAGTAAAACGCATTTTCTTAAAAAGGGAAAAGACAGCACGCTGCTGGCCTCAAAGCTTTTCGCGCCCACGCGTGAAACGGGAAGCGGAAGGCTGTACCGc                                                     >  2:576915/1‑151 (MQ=255)                                                                                                ttGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGGAAAACGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAACGCGGGAAGCGGGAAGCGAATGGGTGTACCGcc                                                    >  1:487917/1‑150 (MQ=255)                                                                                                  ggTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcat                                                  <  1:428410/150‑1 (MQ=255)                                                                                                    ttttGATTACGGTACAGCAAACTGTTCAGTGGCGGCCA‑TTGGTGACGGTAAACCGCATTTGCTAAAAAAGGAAAACGAACGCACCGTGCTGCCTTTAATGCTTTTCGCGCCAACCCCGGAAGCGGGAAGCGAATGGCTGTACCGcc                                                    >  2:770216/1‑146 (MQ=255)                                                                                                     tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcatcat                                               <  1:413269/150‑1 (MQ=255)                                                                                                     tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCcatcat                                               <  2:1019804/150‑1 (MQ=255)                                                                                                        gATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TTCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAAAGGCTGTACCGCCCTCctg                                              >  1:966045/1‑148 (MQ=255)                                                                                                        gATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCatgat                                            <  2:174253/150‑1 (MQ=255)                                                                                                           tACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGGCGGGAAAACGCAATTTCTAAAAATGGGAAACGAGAGCACGCTTCGGCCTTCAATGCCTTGGGCGCCCAAGCCTGAAACGGGATACGAATGGGTGTACCGCcatca                                                >  1:506613/1‑143 (MQ=255)                                                                                                            aCGGTACAGCAAACTGTTCAGTGGCGGTCA‑TTCGGGGAGGGAAAACCCATTTTCTAAAAAAGGGAAAAGACCGCACGCCGCTGCCCTCAATGCTTTTCGCGCCAAAGCGGGAAGCGGTAAGCGGATGGCTGTAACCCCCTCATTATGt                                          >  2:1036932/1‑148 (MQ=255)                                                                                                            aCGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTc                                        <  1:149027/150‑1 (MQ=255)                                                                                                             cGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCg                                      <  1:927996/151‑1 (MQ=255)                                                                                                                 aCAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAg                                  <  2:106500/151‑1 (MQ=255)                                                                                                                      aaaCTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAAACGCATTTGCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGGAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAga                                 >  1:951113/1‑147 (MQ=255)                                                                                                                              cAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAAc                     <  2:559011/151‑1 (MQ=255)                                                                                                                               aGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACg                    <  2:477946/151‑1 (MQ=255)                                                                                                                               aGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACg                    <  1:372876/151‑1 (MQ=255)                                                                                                                                             cGTGACGGTAAACCGAATTTGCTAAAACTGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTAc       >  2:1044852/1‑151 (MQ=255)                                                                                                                                                gACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTc    >  2:1092233/1‑151 (MQ=255)                                                                                                                                                   ggTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCgg  <  2:622309/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCTGCGCTTTGCTTTCATATTCCGGTTGTCGCGACGGCAACATTTCGCTAAAGTCACGCCCCTTCTTCACCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCA‑TGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAAGCGAATGGCTGTACCGCCATCATGATGTTCCGGCAGACGACGATGAAACGCAGGCGCTGTTACGTCGG  >  NC_000913/2147586‑2147880

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