Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,173	0	T	C	40.7%	129.3 / 25.0	81	I310V (ATT→GTT)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base T (40/8);  new base C (0/33);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 2.12e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACC  >  NC_000913/2921028‑2921323                                                                                                                                                   |                                                                                                                                                       aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATcccac                                                                                                                                                   >  1:958565/1‑151 (MQ=255)     ccATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                                 >  1:1135341/1‑149 (MQ=255)       aTAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGca                                                                                                                                             <  1:711829/151‑1 (MQ=255)         aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                            <  1:642490/150‑1 (MQ=255)          aaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCagac                                                                                                                                          <  1:103639/151‑1 (MQ=255)           aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACAcc                                                                                                                                         <  2:1287144/151‑1 (MQ=255)             ccGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgc                                                                                                                                       <  2:640685/151‑1 (MQ=255)              cGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgcc                                                                                                                                      >  2:341347/1‑151 (MQ=255)               gatgatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACAccgccg                                                                                                                                     >  2:25968/1‑151 (MQ=255)                  gatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGAt                                                                                                                                   >  2:411675/1‑150 (MQ=255)                  gatCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGAt                                                                                                                                   >  1:329903/1‑150 (MQ=255)                    tcatcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATaa                                                                                                                                 >  2:465475/1‑150 (MQ=255)                      atcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaac                                                                                                                              >  2:615067/1‑151 (MQ=255)                       tcCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaacc                                                                                                                             >  1:835657/1‑151 (MQ=255)                         cACTCAACGTTGACGTAGTTGCAGAACACCATCCCCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGACCCCGTGCGGCGC‑ATCAGCCGGTCGGAAATAACCCCCCCCAGCACACCGCCGATAaaacc                                                                                                                             <  1:334151/149‑1 (MQ=255)                              aaCGTTGACGTAGTTGCAGAACCCCACCCCCACCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTCCAGCGACCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAACCCCACCCAGCACACCGCCGATAaaaccgcaaa                                                                                                                        <  1:501596/149‑1 (MQ=255)                                       gTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTCCGCGAAATGTTCAGCGACCCCGTGCGGCGC‑ATCAGCCAGTCGAAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaa                                                                                                               <  2:257179/149‑1 (MQ=255)                                       gTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                             >  1:1260988/1‑151 (MQ=255)                                         aGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgg                                                                                                            >  2:1130280/1‑150 (MQ=255)                                           ttacAGAACCCCATCCCCCTGGACAGCAACATGCCCATTCCGATCGGTTTTTCCCGCGAAGGGGCCAGCGCCCCGGAGCGGCGC‑ATCAACCAGTAGAAAATAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                          <  1:649341/147‑1 (MQ=255)                                            tGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGc                                                                                                        >  2:277886/1‑151 (MQ=255)                                              cAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGcc                                                                                                       >  2:311483/1‑150 (MQ=255)                                                   tcccacccccatcGCCAGCAACTTCCCCATTCCGATGGGTGTTTTCCCCGCAATGTGCAGATACCCCGTGGGGCCC‑ATCGCCCTGGCGAAAATAACCCCCCCCAGCACTCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaaa                                                                                                  <  2:236855/142‑1 (MQ=255)                                                     ccaccccccccGACAACAACATGCCCATTAGGGCCGGTGTTTTCCGCCCAATGTCCACCGCCCCCGTGCGGCGC‑ACCAGCCAGCGGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGcc                                                                                               <  1:1235054/141‑1 (MQ=255)                                                     ccacccccaccGACACAAACATCCCCATTACGACCGGTGTTTTACGCGCAATGTCCAGCGACCCCGGGCGGCGC‑ACCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGcc                                                                                               <  1:807550/145‑1 (MQ=255)                                                     ccaccaccatcGACAGCAACATGCCCATTACGATCGGTTTTTTCCGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGcc                                                                                               <  2:685734/151‑1 (MQ=255)                                                      catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACaaa                                                                                                  >  1:236823/1‑147 (MQ=255)                                                        tcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGc                                                                                            >  2:638310/1‑151 (MQ=255)                                                         ccaccccaaCAGAAACAGGCCCATTACTATCGGTGTTTGACGCAAATTGTCCAGCGATCCCGTGCGGCCC‑ATCCCCCAGGCGGAAATAATCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGc                                                                                                <  1:1063899/138‑1 (MQ=255)                                                           ccatcGACAGCAACATCCCCATCCCGATCGGTGTTTCCCGCACAATGTCCAGCGACCCCGGGCGGCGC‑ACCACCCAGTCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                         <  1:768177/151‑1 (MQ=255)                                                           ccatcAACAGCAACATCCCCATTCCGACCGGTGTTTTGCGCGAAATTTCCAGCGCCCCCGTGCGGCGC‑ATCACCCATTCGGAAATAACCCCACCCACCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCt                                                                                           <  2:801251/149‑1 (MQ=255)                                                            catcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                        >  1:521119/1‑151 (MQ=255)                                                              tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                        >  1:735928/1‑149 (MQ=255)                                                               cGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                        >  2:1197498/1‑148 (MQ=255)                                                                 aCAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGaa                                                                                    >  1:372465/1‑150 (MQ=255)                                                                  cAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                  >  2:63320/1‑151 (MQ=255)                                                                   aGCAACATCCCCTTTGCGACCGGTTTTTTCCGCGCAATGTACCCCGACCCGGTG‑GGCACTGTCACCCGTCCAAAAATAACCCCCCCCGGCACTCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                  <  2:646897/150‑1 (MQ=255)                                                                     caccATGCCCATGACGATCGGTGTTTTACGCGAAAGTTCCAGCGATCCGGTGCGGCCC‑ATCAGCAAGCCGGAAAAAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                  <  1:8726/145‑1 (MQ=255)                                                                      accagccccATCCCGATCGGTGTTTTCCGCCCAATGTCCAGCGACCCCGGGCGGCGC‑ACCACCCAGCCGAAAATAACCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATc                                                                                  <  1:1165354/141‑1 (MQ=255)                                                                      aaCATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                               >  2:253737/1‑150 (MQ=255)                                                                      aaCATGCCCATGACGTAGGGTGTTTTACGCGCAGTGTACAGCGACCCAGTCGGGCAC‑TTAAGCCAGTCGGAAAAAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGAc                                                                               <  2:302529/150‑1 (MQ=255)                                                                       acacGCCCATTACGATCGGTGTTTTACGCGAAATTTACAGCGACCCCGTGGGGCGC‑ACCAGCCAGTCGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCg                                                                                 <  1:776270/147‑1 (MQ=255)                                                                         ctgcctactgcGATGGGTTTTTTACGCGAAATGTACAGCGACCCCGCGCGGCGC‑ACCAGCCACTGGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATc                                                                            <  1:529306/140‑1 (MQ=255)                                                                            cccATTACGATCGGTTTTTTCCGCGCAATGTCCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc                                                                          <  2:789387/149‑1 (MQ=255)                                                                            cccATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcgc                                                                        >  2:1199749/1‑151 (MQ=255)                                                                            cccATTACGACCGTTGTTTTACGCGAATTGTCCAGCGACCCCGGGGGGCCC‑ACCAGCCAGCCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                         <  1:870891/150‑1 (MQ=255)                                                                              cATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGt                                                                      >  1:33880/1‑151 (MQ=255)                                                                                ttACGATCGGTGTTTTACGCGAATTTTTCAGCGCTCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTg                                                                     <  2:892263/150‑1 (MQ=255)                                                                                    gttCGGTGTTTTCCGCGCAATGTCCAGCGACCCCGGGCGGCGC‑ACCAGCCAGCCGAAAAAAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                                  <  1:447169/147‑1 (MQ=255)                                                                                    gATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                                >  1:1203160/1‑151 (MQ=255)                                                                                     aTCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                                >  2:226024/1‑150 (MQ=255)                                                                                     aTCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                                >  1:125817/1‑150 (MQ=255)                                                                                             tttACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATa                                                        >  2:1170072/1‑150 (MQ=255)                                                                                             tttACGCGAAATGTTCAGCGACCCCGTGCGGCGC‑ATCACCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATa                                                        <  1:890380/150‑1 (MQ=255)                                                                                             tttACGCGAAAGTTTCAGCGACCCGGTGCGGCGC‑ATCAGCCAGTGGAAAATAACCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGAt                                                         <  1:469928/149‑1 (MQ=255)                                                                                               tACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAc                                                     >  1:1178292/1‑151 (MQ=255)                                                                                               tACGCGAAATGTTCAGCGATCCCGGGCGGCGC‑ATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaaa                                                      <  2:453578/150‑1 (MQ=255)                                                                                                  gcgcAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCgg                                                  >  1:1235567/1‑151 (MQ=255)                                                                                                      aaTGTTCAGCGACCCCGGGCGGCCC‑ATCACCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAc                                              <  2:491746/151‑1 (MQ=255)                                                                                                         gTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATATAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAAcc                                             >  1:100791/1‑149 (MQ=255)                                                                                                           tCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                         >  1:54306/1‑151 (MQ=255)                                                                                                            cAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGt                                         >  1:8358/1‑150 (MQ=255)                                                                                                               cGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                      >  2:579726/1‑150 (MQ=255)                                                                                                                gaccccTTGCGGCGC‑ACCAGCCGGTCAAAAAAAACCCCCCCCAACACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaat                                      <  1:373122/146‑1 (MQ=255)                                                                                                                   cccGTGCGGCGC‑TTCAGCCAGTCGAAAATACCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                                     <  1:91746/147‑1 (MQ=255)                                                                                                                   cccGGGCGGCCC‑ATCAGCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaata                                     <  1:614966/147‑1 (MQ=255)                                                                                                                        gcggcgC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                             >  1:377573/1‑150 (MQ=255)                                                                                                                            cgc‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                         >  2:626612/1‑150 (MQ=255)                                                                                                                             gc‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                         >  1:90167/1‑149 (MQ=255)                                                                                                                                accagccggggggAAATAACCCCCCCAACCCCGCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAgg                      <  1:64282/139‑1 (MQ=255)                                                                                                                                accagcacatCGAAAAAAACCCCCCCCGGCACACCGCCGATAAACCCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagta                             <  1:1282845/134‑1 (MQ=255)                                                                                                                                aTCAGCCAGTCGAAAATAACCCCCCCCAGCCCCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAgg                      <  2:868795/150‑1 (MQ=255)                                                                                                                                    gCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt                 >  2:495366/1‑151 (MQ=255)                                                                                                                                    gCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt                 >  2:312686/1‑151 (MQ=255)                                                                                                                                     ccAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt                 <  2:299194/150‑1 (MQ=255)                                                                                                                                       aGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGAt              >  1:852632/1‑151 (MQ=255)                                                                                                                                           gaaaaTAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc            <  2:644283/147‑1 (MQ=255)                                                                                                                                             aaaaaaCCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTa        <  1:341293/151‑1 (MQ=255)                                                                                                                                             aaaTACCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATAc            <  1:372277/147‑1 (MQ=255)                                                                                                                                             aaaTAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt         >  2:1289439/1‑150 (MQ=255)                                                                                                                                             aaaTAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGt         <  1:253705/150‑1 (MQ=255)                                                                                                                                                   tCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAcc  >  2:1195472/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGC‑ATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACC  >  NC_000913/2921028‑2921323

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