Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,921,177	0	A	C	41.2%	130.5 / 11.8	68	G308G (GGT→GGG)	gudP	putative D‑glucarate transporter
Reads supporting (aligned to +/- strand):  ref base A (33/7);  new base C (0/28);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.98e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCG  >  NC_000913/2921028‑2921324                                                                                                                                                      |                                                                                                                                                    aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccac                                                                                                                                                    <  1:882304/151‑1 (MQ=255) aGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccac                                                                                                                                                    >  1:601512/1‑151 (MQ=255)   cgcCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATcccaccc                                                                                                                                                  >  2:565504/1‑151 (MQ=255)        tAAAGCCGATGATCAAACACTCAACGTTAACGTAGTTTCAGAACACCAACACCAACGACAGCAACATGCCCATTACGATTGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGca                                                                                                                                              <  2:204585/150‑1 (MQ=255)         aaaGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                             <  1:879183/150‑1 (MQ=255)           aGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACAcc                                                                                                                                          <  1:606274/151‑1 (MQ=255)              cGATGATCATCCTCTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGcac                                                                                                                                             <  2:999488/145‑1 (MQ=255)                     cacccACTCAACGTGGACATAGTTGCAGAACACCATCCCCCTCGAAACCAACACGCCCATTACGATCGGTTTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCAGAAGCCAGTGGAAAATAATCCCCCCCAGCCCCCCGCCGATa                                                                                                                                   <  2:408130/145‑1 (MQ=255)                      atcAACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaac                                                                                                                               >  2:1184135/1‑151 (MQ=255)                        cccctaaaCGTGGCCGTATTTCCAGAACCCCACCACCTTCAACACAAACTCGCCCATTCCGACCGGTTTTTTCCGCGAAATTTCCAGCGACCCCGTGCGGCGCATCAGCCAGTGGGAAATAATCCCCCCAAGCCCACCGCCGATAaaaccg                                                                                                                             <  2:1141947/145‑1 (MQ=255)                        ccACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaacc                                                                                                                              >  1:980966/1‑150 (MQ=255)                         cACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                            >  2:140988/1‑151 (MQ=255)                         cACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgc                                                                                                                            >  2:1182806/1‑151 (MQ=255)                           cTCAACGTTGATTTATTTGCAGAACACCATCCCCATCGCCACCAACATGCCCATTACGATCGGTGTTTTACGCGCAAGTTTCAGCGATCCCGTGCGGCGCATCCGCCGGTCGGAAATAATCCCCCCCAGCACACCGCCGATAaaacc                                                                                                                              <  2:1050817/147‑1 (MQ=255)                             cAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaac                                                                                                                        >  1:75829/1‑151 (MQ=255)                               aCGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAaaaccgcaaaccg                                                                                                                      >  1:923659/1‑151 (MQ=255)                                     aCGTAGTTGCAGAACACCATCCCCATCGACAGCAACATGCCCTTTACGATGGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCcggaa                                                                                                                <  1:558200/151‑1 (MQ=255)                                        tAGTTACAAAACCCCATCACCTTCAACAGCAACATCCCCATTCCGATCGGTGTTTTACGCGAAATTTTCAGCGACCCCGGGCGGCGCACCCCCCAGTCGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCcggaacg                                                                                                              <  2:1001544/150‑1 (MQ=255)                                          gTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCcggaacgga                                                                                                            >  1:1167/1‑150 (MQ=255)                                               aGAACCCCATCCCCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTCCAGCGACCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCca                                                                                                       <  1:157647/150‑1 (MQ=255)                                                 aaCCCCACCCCCACCGACAGCAACATGCCCATTACGACCGGTTTTTTCCGCGAAATTTCCGGCGCCCCCGTGCGGCGCACCAGCCAGTCGAAAATAACCCCCCCCAGCCCACCGCCGATAAAACCGCAAACCGCcggaacggaa                                                                                                           <  2:248550/144‑1 (MQ=255)                                                      catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                               >  2:701433/1‑151 (MQ=255)                                                      catcaccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGccc                                                                                               >  1:785850/1‑151 (MQ=255)                                                         ccccatcGACAGCAACATGCCCATTACGATCGGTTTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCg                                                                                              <  1:653915/147‑1 (MQ=255)                                                           ccatcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCttt                                                                                          >  1:162240/1‑151 (MQ=255)                                                              tcGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                         >  2:793740/1‑149 (MQ=255)                                                               cGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTc                                                                                         >  2:1087403/1‑148 (MQ=255)                                                                      aaCATGCCCATTACGATGGGTTTTTTCCGCGAAATGTTCAGCGATCCCGGGCGGCGCATCAGCCAGTGGAAAATAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGaa                                                                                     <  2:902149/145‑1 (MQ=255)                                                                      aaCATCCCCATTACGATGGGTGTTTTCCGCGAAATGTTCAGCGCCCCCGGGGGGCGCATCAGCCAGTGGAAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCg                                                                                  <  2:1146954/148‑1 (MQ=255)                                                                          tGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATccc                                                                           >  1:1177942/1‑151 (MQ=255)                                                                            cccTTTCCGATCGGTGTTTTACGCGCATTGTCCAGCGACCCCGTGCGGCCCATCACCCAGTGGGAAATAACCCCCCCCCCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                          <  2:895510/150‑1 (MQ=255)                                                                            cccATTACGATGGGTGTTTTCCGCGCAATTTTCAGCGACCCCGTGCGGCCCATCAGCCAGTGGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCcg                                                                          <  2:1105876/150‑1 (MQ=255)                                                                                 tACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGcc                                                                    >  1:972629/1‑151 (MQ=255)                                                                                  aCGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                                   >  2:588475/1‑151 (MQ=255)                                                                                  aCGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCt                                                                   >  1:1061619/1‑151 (MQ=255)                                                                                     aTCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGc                                                                 >  1:50120/1‑150 (MQ=255)                                                                                      tCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAc                                                               >  2:356221/1‑151 (MQ=255)                                                                                        ggTGTTTAACGCGAAATTTTCACGGACCCGGTGCGGCGCACCGGCCAGGCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCa                                                                <  2:549172/148‑1 (MQ=255)                                                                                          tgtTTTACGCGCAATGTTCAGCGACCCCGTGCGGCCCATCAGCCAGTCGGAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAg                                                            <  1:365631/150‑1 (MQ=255)                                                                                           gtTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCAcc                                                              >  1:704903/1‑147 (MQ=255)                                                                                              ttccgcgCATTGTCCAGCGACCCCGTGCGGCGCATCAGCCAGGCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa                                                        <  2:139515/147‑1 (MQ=255)                                                                                              ttACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATaa                                                        >  1:292582/1‑150 (MQ=255)                                                                                                aCGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                                     >  1:275728/1‑151 (MQ=255)                                                                                                aCGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAAcc                                                     >  2:478844/1‑151 (MQ=255)                                                                                                    gaaaggTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATACCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCggg                                                  <  2:198438/145‑1 (MQ=255)                                                                                                     cAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGa                                                 >  1:345606/1‑150 (MQ=255)                                                                                                      aaTGTCCAGCGACCCCGGGCGGCGCATCAGCCAGTCGAAAATAACCCCACCCAGCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGaa                                                <  2:328932/150‑1 (MQ=255)                                                                                                                gATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGCGCCTGCACCAGATAAACCGGGAACCAGGtaata                                      >  1:836392/1‑150 (MQ=255)                                                                                                                  tCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGtaataa                                     >  1:317257/1‑149 (MQ=255)                                                                                                                  cccggtggggccaaTCACCCAGACGAAAAAAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaag                                   <  1:370590/138‑1 (MQ=255)                                                                                                                   cccGTGCGGCGCATCAGCCAGTCGGAAATAACCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAaaga                                  <  1:444782/151‑1 (MQ=255)                                                                                                                      gTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagt                               >  2:762579/1‑151 (MQ=255)                                                                                                                          ggcgCACCAGCCAGTCGGAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                           <  1:659212/151‑1 (MQ=255)                                                                                                                          ggcgCACCAGCCAGCCGGAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                           <  2:244497/151‑1 (MQ=255)                                                                                                                          ggccCACCAGCCAGTGGGAAAAAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGaagtaagt                           <  1:845648/151‑1 (MQ=255)                                                                                                                               aTCAGCCAGTCGAAAATAACCCCCCCCACCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAg                        <  1:77998/149‑1 (MQ=255)                                                                                                                               aTCACCCAGGCGAAAATAACCCCCCCCACCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAg                        <  2:401880/149‑1 (MQ=255)                                                                                                                                  aGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGt                   >  1:53618/1‑151 (MQ=255)                                                                                                                                  aGCCAGTCGAAAATAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTc                          <  2:479375/144‑1 (MQ=255)                                                                                                                                    ccAGTCGAAAATAATCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAggg                      <  1:388915/146‑1 (MQ=255)                                                                                                                                    ccAGTCGAAAAAAATCCCCCCCACCCCACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGtt                  <  1:143854/150‑1 (MQ=255)                                                                                                                                     cAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTg                 >  2:1182392/1‑150 (MQ=255)                                                                                                                                         cgaaaaTAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACa            <  1:253521/148‑1 (MQ=255)                                                                                                                                          gaaaaTAATCCCCCCCAGCACCCCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAg           <  1:128925/149‑1 (MQ=255)                                                                                                                                          gaaaaTAACCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAg           <  1:381827/149‑1 (MQ=255)                                                                                                                                            aaaaaaCCCCCCCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTa         <  1:1087431/151‑1 (MQ=255)                                                                                                                                             aataatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTAc        >  1:1016112/1‑151 (MQ=255)                                                                                                                                             aataatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTAc        >  2:739388/1‑151 (MQ=255)                                                                                                                                                aatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg      >  2:826904/1‑150 (MQ=255)                                                                                                                                                aatCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTg      >  2:119927/1‑150 (MQ=255)                                                                                                                                                 atCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAc    >  1:700157/1‑151 (MQ=255)                                                                                                                                                  ccccccccaCCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGAcc   <  2:601493/150‑1 (MQ=255)                                                                                                                                                   cccacccaGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCg  >  1:1081090/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AGCGCCATAAAGCCGATGATCATCCACTCAACGTTGACGTAGTTGCAGAACACCATCACCATCGACAGCAACATGCCCATTACGATCGGTGTTTTACGCGCAATGTTCAGCGATCCCGTGCGGCGCATCAGCCAGTCGGAAATAATCCCACCCAGCACACCGCCGATAAAACCGCAAACCGCCGGAACGGAAGCCACAAAGCCCGCTTTCAGAATCGACATCCCGCGTGCCTGCACCAGATAAACCGGGAACCAGGTAATAAAGAAGTAAGTCAGGGCGTTGATACAGTACTGACCG  >  NC_000913/2921028‑2921324

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