Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,181,603	0	C→T	11.7%	11.5 0.0 ‑0.3	61	intergenic (+22/+16)	ygiC/ygiD	ATP‑Grasp family ATPase/4,5‑ DOPA‑extradiol‑dioxygenase
Reads supporting (aligned to +/- strand):  new base (4/3):  ref base (27/26):  total (31/30)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.30e-01

GTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTT  >  NC_000913/3181397‑3181812                                                                                                                                                                                                               |                                                                                                                                                                                                                  gTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGcc                                                                                                                                                                                                                >  1:535832/1‑210 (MQ=255)  tCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTCTAATAATCCACGATACCGGACGGCACTCGCCa                                                                                                                                                                                                               >  2:185423/1‑210 (MQ=255)     gTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATcc                                                                                                                                                                                                            <  1:297313/210‑1 (MQ=255)      tATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCg                                                                                                                                                                                                           <  2:535804/210‑1 (MQ=255)       acacgacgctcttccgatctGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                                                                                          <  2:580390/191‑1 (MQ=255)       aaggCGGAGAAGGGATGACTGTTCAGCAAGTCCACCCGTGACCGATATTCGGCGACAGCTATATGCTGATTGGTGGCGGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                                                                                          <  2:635522/208‑1 (MQ=255)       aTGGCGAAGAAGGTATGATTGTTCAGCAATTCCGCCCGTTACAGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCgg                                                                                                                                                                                                          <  2:127994/210‑1 (MQ=255)       aTGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCgg                                                                                                                                                                                                          <  1:216464/210‑1 (MQ=255)          gCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTaa                                                                                                                                                                                                       <  1:52343/210‑1 (MQ=255)                           gTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCGAGTGCCATCCggagatcggaagagcacacgtc                                                                                                                                                                                      >  1:580422/1‑190 (MQ=255)                             tCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCAcc                                                                                                                                                                                    <  1:289726/210‑1 (MQ=255)                                gCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGAc                                                                                                                                                                                 >  1:118816/1‑210 (MQ=255)                                      ccACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATa                                                                                                                                                                           >  1:422810/1‑210 (MQ=255)                                         cccGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTc                                                                                                                                                                        >  1:100557/1‑210 (MQ=255)                                               aCCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTg                                                                                                                                                                  <  2:508678/210‑1 (MQ=255)                                                         ggCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCta                                                                                                                                                        >  1:420706/1‑210 (MQ=255)                                                              cAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATAccc                                                                                                                                                   <  1:250081/210‑1 (MQ=255)                                                                 cTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCa                                                                                                                                                <  2:118812/210‑1 (MQ=255)                                                                 cTATATGCTAATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCa                                                                                                                                                <  1:8307/210‑1 (MQ=255)                                                                                tAGCTGGCTGGTGACCGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAAtt                                                                                                                                 <  2:115321/210‑1 (MQ=255)                                                                                     ggctggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTc                                                                                                                            >  1:522269/1‑210 (MQ=255)                                                                                       ntggTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTCTACCCTCAACCGGAATGGTAATTGGCTCCt                                                                                                                          >  2:101479/2‑210 (MQ=255)                                                                                          gTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCCGcc                                                                                                                       >  2:533811/1‑210 (MQ=255)                                                                                            gAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATCCTCAGGCTGCCCATTTCTATACCCTCAACCGGACTGGTAATTGGCTCCTGCCCa                                                                                                                     >  2:409636/1‑210 (MQ=255)                                                                                              aCGATCAACCCGCAGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATc                                                                                                                   <  2:602861/210‑1 (MQ=255)                                                                                                gATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATccc                                                                                                                 <  2:100553/210‑1 (MQ=255)                                                                                                 aTCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGTCTCCTGCCCATCCCa                                                                                                                >  2:176048/1‑210 (MQ=255)                                                                                                 aTCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCCCCTGCCCATCCCa                                                                                                                >  2:77981/1‑210 (MQ=255)                                                                                                   cAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCAGTGATCACCCAGGATATTTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACg                                                                                                              <  2:71491/210‑1 (MQ=255)                                                                                                     gtcgtccgGTATTGGCATTCGTGAAGCCCGTGCATTGATCACCCAGGATATGTCTCGCTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTAGCCATCCGGTAATTGTTAGCCTATCTGCACAGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCa                                                                                                            <  1:6958/205‑1 (MQ=255)                                                                                                      cccGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCAc                                                                                                           >  1:395051/1‑210 (MQ=255)                                                                                                           cgGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCCCCTAac                                                                                                      >  2:414366/1‑210 (MQ=255)                                                                                                                        cGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACg                                                                                         <  1:658135/210‑1 (MQ=255)                                                                                                                           gAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCa                                                                                      >  1:661986/1‑210 (MQ=255)                                                                                                                               aCCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCCCGCACCTAACACATACAACAACGGCAaata                                                                                  >  1:342355/1‑210 (MQ=255)                                                                                                                                       ttGGTCACCCATGGTATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGGTGGCCCTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTg                                                                          <  2:422792/210‑1 (MQ=255)                                                                                                                                        tGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAATACATACAACAACGGCAAATAATGCTCTgg                                                                         <  2:265761/210‑1 (MQ=255)                                                                                                                                         gATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATGGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCCCTGGc                                                                        >  1:322171/1‑210 (MQ=255)                                                                                                                                             aCCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTg                                                                    <  2:599614/210‑1 (MQ=255)                                                                                                                                              cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTgg                                                                   <  2:666996/210‑1 (MQ=255)                                                                                                                                              cccAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTgg                                                                   >  1:605419/1‑210 (MQ=255)                                                                                                                                                         tctcGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACg                                                        <  1:62804/210‑1 (MQ=255)                                                                                                                                                                tttATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCAcc                                                 <  1:330549/210‑1 (MQ=255)                                                                                                                                                                      cacaTATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATg                                           >  2:470213/1‑210 (MQ=255)                                                                                                                                                                            tttttGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAg                                     >  1:486317/1‑210 (MQ=255)                                                                                                                                                                                   tgAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTc                              >  2:657215/1‑210 (MQ=255)                                                                                                                                                                                    gAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCa                             <  2:402178/210‑1 (MQ=255)                                                                                                                                                                                      ctAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCAcc                           <  2:415722/209‑1 (MQ=255)                                                                                                                                                                                       tAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCa                          >  1:491502/1‑210 (MQ=255)                                                                                                                                                                                          gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCagag                       >  2:464823/1‑210 (MQ=255)                                                                                                                                                                                          gCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCagag                       <  1:414383/210‑1 (MQ=255)                                                                                                                                                                                                aTACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATCCACCAGAGGAtgtt                 >  1:596280/1‑210 (MQ=255)                                                                                                                                                                                                  aCCGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACACCAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGAtgttgt               >  2:39754/1‑210 (MQ=255)                                                                                                                                                                                                  aCCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGAtgttgt               >  2:218979/1‑210 (MQ=255)                                                                                                                                                                                                   ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATTCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGGAATTAACCAGCGGCTgttgtt              >  2:137990/1‑210 (MQ=255)                                                                                                                                                                                                   ccGGATGGCGAGTGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATgttgtt              >  2:84783/1‑210 (MQ=255)                                                                                                                                                                                                     ggATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGAGGTTGTTcc            >  1:22350/1‑210 (MQ=255)                                                                                                                                                                                                       aTGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCAc          <  2:368111/210‑1 (MQ=255)                                                                                                                                                                                                        tGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACt         <  2:22349/210‑1 (MQ=255)                                                                                                                                                                                                           cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGc      <  2:605384/210‑1 (MQ=255)                                                                                                                                                                                                           cACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTGGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGc      >  2:149212/1‑210 (MQ=255)                                                                                                                                                                                                               cGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCAATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACCACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCtt  <  1:657252/210‑1 (MQ=255)                                                                                                                                                                                                               |                                                                                                                                                                                                                  GTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCGAAATTCGGCGACAGCTATATGCTGATTGGTAGCTGGCTGGTGAACGATCAACCCGCCGGAATTGGCATTCGTGAAGACCGTGCATTGATCACCCAGGATATGTCTCGGTTTTATCCACATATTTTTGTTGAATAAGCCACGATACCGGATGGCACTCGCCATCCGGTAATTGTTAGCCTATCTGCACCGACAGCATACTCAGGCTGCCCATTTCTATACCCTCAACCGGAATGGTAATTGGCTCCTGCCCATCCCACGCACCTAACACATACAACAACGGCAAATAATGCTCTGGCGTTGGGTTCGATAACGTGCCACCTTCATGGTCGAGGTAATTCACCAGAGGATGTTGTTCCACTGGCCCTT  >  NC_000913/3181397‑3181812

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.