Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	52.1%	86.8 / 129.0	71	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (29/5);  new base C (0/37);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 4.06e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-13

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589                                                                                                                                                      |                                                                                                                                                     gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                     >  2:531038/1‑151 (MQ=255) gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCacac                                                                                                                                                     >  2:604530/1‑151 (MQ=255)  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAcc                                                                                                                                                    <  1:421545/151‑1 (MQ=255)     atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                                                                 <  2:455406/151‑1 (MQ=255)      tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                >  2:776659/1‑151 (MQ=255)       gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                <  1:658160/150‑1 (MQ=255)       gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                                                                               >  1:286384/1‑151 (MQ=255)            gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccg                                                                                                                                          >  1:1212794/1‑151 (MQ=255)             agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgc                                                                                                                                         <  2:878679/151‑1 (MQ=255)              gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgcc                                                                                                                                        <  1:786924/151‑1 (MQ=255)               aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCccgccg                                                                                                                                       >  2:764971/1‑151 (MQ=255)                    aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                   >  2:191507/1‑150 (MQ=255)                          gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                            >  2:1114721/1‑151 (MQ=255)                            gggAGCCCAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTg                                                                                                                          <  1:693921/151‑1 (MQ=255)                                  ccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                    >  1:959645/1‑151 (MQ=255)                                    aGAAAATGTATGTCGCCTTGCCACCCATCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTc                                                                                                                    <  2:770266/149‑1 (MQ=255)                                      aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                >  2:951703/1‑151 (MQ=255)                                        aaTGTATGTCGACTTGCCAGCCATCATTATCAAACAAATTAAGGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                              <  1:471986/151‑1 (MQ=255)                                          tgtatgtCGCCTTGCCCGCCATCATTATCAAAAAGATTAAGGCCGGCAGATTTTTATGGATAAGGAGAAAAGGCAGATCTTTTCTCCCCTTTAAAGGGACAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                             <  1:1266769/150‑1 (MQ=255)                                            tatgtCGTCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGAAGAGCTTTTCGTCACTTTAAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                            <  2:353026/149‑1 (MQ=255)                                            tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACTTTAAATGGTCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                            <  1:179675/149‑1 (MQ=255)                                                     ttcccacccaTCTTTATTAAAAAAATTAATGCCGGAAAATTTTAAGGGGTAGGGAGAAAAGGCAGAGTTTTTCGTCCCTTTAAATAGTCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  1:435004/147‑1 (MQ=255)                                                     ttcccacccaTATTTATTAAACAGATTAATGCCGGAATATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAAGGGTAAGCCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                    <  1:558957/145‑1 (MQ=255)                                                     ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACGTTAAATGTTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                    <  1:666555/148‑1 (MQ=255)                                                     ttGCCAGCCATAATTATTAAACAGATAAATGACGAGATATTTTTATGGGAAAGGAGAAAATGCAGATCTTTTCGTCACTTTAAAGGGTAAAGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                    <  1:463874/148‑1 (MQ=255)                                                       cccacccaTCATTATTAAAAAGATAAATGACGGCAGATTTTTATCGGAAAGGAGAAAATGCAGAGTTTTTCGTCCCGTTCAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                  <  1:951693/147‑1 (MQ=255)                                                        ccaccctTCTTTATTCAACAGTTAAATGCCGGAAAATTTTTATGGGAAAGGAGAAAATGCAGATCTTTTCGTCCCGTTAATTGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                 <  2:479127/148‑1 (MQ=255)                                                        ccaaccaTAATTATTCAACAGATTAATGACGGAATATATTTATCGGAAAGGAGGAAAGGAAGGGTTTTTCGTCACGTTAAATGGTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                               <  1:1010087/150‑1 (MQ=255)                                                         cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                             >  2:1288431/1‑151 (MQ=255)                                                          agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                             >  1:1066732/1‑150 (MQ=255)                                                          agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                             >  1:222368/1‑150 (MQ=255)                                                            ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                          >  2:836303/1‑151 (MQ=255)                                                               tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                        >  2:1305959/1‑150 (MQ=255)                                                                  ttattCAACAGATAAATGCCGGAAGATTTTTATCGGAAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                         <  1:302022/146‑1 (MQ=255)                                                                   ttttCAACAGTTTAATGGCGGCAGATATTTATGGGAAAGGAGAAATGGCAGAGCTTTTCGTCACTTTCAAGGGTAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                      <  2:837116/146‑1 (MQ=255)                                                                   tattCAACAGATTAATGACGGAAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTTAAATGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                       <  1:1277685/147‑1 (MQ=255)                                                                    attCAACAGATTAATGACGGCAGATATTTATCGTTAAGGAGAAAATGCAGAGCTTTTCGCCACGTTCAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                  <  2:1055152/151‑1 (MQ=255)                                                                        aaaaaaTTAATGACGGCAAATTTTTATCGGAAAGGAGAAAATGCAGAGTTATTCGTCACGTTAAGGGGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                <  1:292141/144‑1 (MQ=255)                                                                        aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                 >  2:202976/1‑148 (MQ=255)                                                                        aaCAGATTAAGGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGGGTTATTCGTCACTTTAAATGGTCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                  <  2:349342/147‑1 (MQ=255)                                                                                aaggACGGAAGATTTTTAGCGGTAAGGAGGAAAGGCAGAGCTATTCGTCACGTTAAATGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                       <  1:1127808/147‑1 (MQ=255)                                                                                aaTGCCGGCAGATTTTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCTTTAAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                       <  2:1279641/150‑1 (MQ=255)                                                                                aaTGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTATTCGTCCCTTTAAATGTTCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                       <  2:281051/150‑1 (MQ=255)                                                                                 aTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                     >  1:575150/1‑151 (MQ=255)                                                                                  tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                     >  1:200809/1‑150 (MQ=255)                                                                                    aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                  >  2:1167820/1‑151 (MQ=255)                                                                                         aGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCTCCACTTTAAAGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                               <  1:1219583/149‑1 (MQ=255)                                                                                            tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                          >  2:127105/1‑151 (MQ=255)                                                                                            tatTTATCGGAAAGGAGAAAATGAAGAGTTTTTCTTCACTTTAAATGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                             <  1:937664/148‑1 (MQ=255)                                                                                               ttgtgggTAAGGGGGAAAGGCAGAGCTTTTCGCCACGTTAAAGGGTCAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                        <  1:295661/145‑1 (MQ=255)                                                                                                 accGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAAGGGTCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                      <  1:459429/148‑1 (MQ=255)                                                                                                 aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                       >  2:314711/1‑149 (MQ=255)                                                                                                   cGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                   >  1:852085/1‑151 (MQ=255)                                                                                                         ggaggaAATGAAGAGCTTTTCGTCACGTAAAAGGATCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                    <  1:1003481/144‑1 (MQ=255)                                                                                                           agaaaaTGCAGGGCTTTTCGCCACGTTAAATGGTCAAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                             <  1:939508/146‑1 (MQ=255)                                                                                                            ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                          >  1:365931/1‑151 (MQ=255)                                                                                                               aaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                        >  1:969042/1‑150 (MQ=255)                                                                                                                 tGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                     >  2:825054/1‑151 (MQ=255)                                                                                                                      agCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                >  2:232057/1‑151 (MQ=255)                                                                                                                        cTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                              >  1:1220483/1‑151 (MQ=255)                                                                                                                           ttcgccagtttaaaTGGTAAAGCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                     <  1:67466/129‑1 (MQ=255)                                                                                                                           ttCGTCACGTTAAAGGGTAAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                             <  2:652393/149‑1 (MQ=255)                                                                                                                              gTCACTTTAAATGGTCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                           <  2:418971/148‑1 (MQ=255)                                                                                                                                 aCGTTAAAGGGACAAGGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                           <  2:222371/145‑1 (MQ=255)                                                                                                                                   tttCAAGGGTCAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca                     <  1:660069/148‑1 (MQ=255)                                                                                                                                    ttaaaTGGTCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                       <  1:397120/143‑1 (MQ=255)                                                                                                                                    ttCAATGGTAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGctggtttctctgg                   <  1:202973/150‑13 (MQ=255)                                                                                                                                     tCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                 >  2:502161/1‑151 (MQ=255)                                                                                                                                           gTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat           >  2:359697/1‑151 (MQ=255)                                                                                                                                              aaTGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         <  2:1211418/150‑1 (MQ=255)                                                                                                                                              aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         >  1:304533/1‑150 (MQ=255)                                                                                                                                              aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg         >  2:793853/1‑150 (MQ=255)                                                                                                                                                  ccACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg    <  1:86159/151‑1 (MQ=255)                                                                                                                                                   caccCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa   <  1:610032/151‑1 (MQ=255)                                                                                                                                                    acacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAc  >  1:972756/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGAC  >  NC_000913/2984292‑2984589

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