Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,725,233	G→C	14.2%	intergenic (‑65/‑464)	U712_19205 ← / → U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,725,233	0	G	C	14.2%	158.6 / 3.6	65	intergenic (‑65/‑464)	U712_19205/U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base G (40/15);  new base C (2/8);  total (42/23)
Fisher's exact test for biased strand distribution p-value = 2.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑TTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTT  >  CP006881/3725101‑3725369                                                                                                                                          |                                                                                                                                         cTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGcccg                                                                                                                                      <  5:310645/138‑2 (MQ=255)      cGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGcccgcaaga                                                                                                                                 >  2:266343/1‑132 (MQ=255)         cGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTTGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCCCAAGaaaa                                                                                                                              >  5:122742/1‑138 (MQ=255)              cccAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTATAAGCCGCCGCCCGCAACAAAAGACGGc                                                                                                                        >  2:80637/1‑139 (MQ=255)                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                                   >  4:378708/1‑139 (MQ=255)                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCAGCAAGAAAAGACGGCCGaat                                                                                                                   >  8:429851/1‑137 (MQ=255)                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CGTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                                   >  1:29366/1‑139 (MQ=255)                   ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCAGC‑CGTTT‑T‑A‑TCCGCCCGGTTTCTCAGCCGCACCCCCCAAGAAAAGAAGGCCgcgg                                                                                                                   >  1:160090/1‑136 (MQ=255)                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCCGGCAGC‑C‑‑TT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttc                                                                                                            >  3:100804/1‑139 (MQ=255)                        caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                                              >  8:126288/1‑139 (MQ=255)                         aaaaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttctttt                                                                                                             >  4:382248/1‑139 (MQ=255)                             aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTATTTTTGGCAGGCAGC‑CTTTT‑T‑A‑ACAGCCCGTTTTCTCATCCCCACCCCGCAACAAAAGACCGCCGATTTATTTTCACg                                                                                                         >  7:160644/1‑139 (MQ=255)                                         cccGAAAAGAAACAAGCAACCGTTTCAAAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                                                             <  1:408660/139‑1 (MQ=255)                                         agcttAAAGAAACAAGCAATGGTTTCATAGAAAAACATCCCAGCGATTATTTTTGGCAGGCAGC‑T‑TTT‑T‑ATGCCCCCCTTTTTCCCCGCCGCCGCCCGCACGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGc                                                                                             <  8:123947/135‑1 (MQ=255)                                          ccGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGccc                                                                                                                                       >  4:367072/1‑96 (MQ=255)                                          ccGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGccc                                                                                                                                       <  3:367072/96‑1 (MQ=255)                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCCGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGGCGAGTTCTTTTGTCGGAAACTAGGGGctc                                                                                           >  1:279595/1‑137 (MQ=255)                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                           >  7:136224/1‑137 (MQ=255)                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCggggc                                                                                           >  7:397391/1‑134 (MQ=255)                                             tAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                                           >  8:242218/1‑137 (MQ=255)                                                acaTACAAGCAAACTTTTCATACAAAAACATCCCTCCGCTTCTTTTTGGCCGGCTGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCCGCCGCCGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                      <  7:429851/137‑1 (MQ=255)                                                acaTACAACCAACCTTTGCATAAAAAAACACCCCACCGATTATTTTGGCCCGGCTCC‑A‑TTT‑TAA‑GCAGCCCGTTTCCTCCCCCGCACCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                      <  5:140466/137‑1 (MQ=255)                                                aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                      <  1:80637/139‑1 (MQ=255)                                                aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                      <  7:126288/139‑1 (MQ=255)                                                  aTACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCaa                                                                                     >  1:49152/1‑138 (MQ=255)                                                   tACAAGCAAACGTTTCATAGAAAAAAATCCCTCCCCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GAAGCCCGTTTTCTCAGCAGCAGCCCGCAAGAAAAGACGGCCGATTTCTTATCTAGCAAACAATCGGCTCTAGCaac                                                                                    >  7:22241/1‑137 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  3:124657/1‑139 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  7:345503/1‑139 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  3:299563/1‑139 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  3:118354/1‑139 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCCGCTCAAGCaaaa                                                                                   >  8:282559/1‑138 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTAATTTATAGCAAACAAGCGGCTCTAGCaaaa                                                                                   >  8:175919/1‑138 (MQ=255)                                                    aCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCGTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                                  >  6:120380/1‑139 (MQ=255)                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTGTTTTTTGCAGGCATCGT‑TTT‑T‑A‑GCAGCCCGTTTTATAAGCCGCAGCCCGCAAGAAAAGACGGACGATTTCTTGTGTGTCAAACATGCGGTTCGAACaaaaat                                                                                 >  6:353776/1‑136 (MQ=255)                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACTAGCGGCTCAAGCCAAACt                                                                                 >  4:103100/1‑138 (MQ=255)                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTGTAGCCAACTATCGGCGCTAGCAAAAca                                                                                 >  1:230933/1‑137 (MQ=255)                                                      aaGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAATACGGCCGAGTTCTTTTCTCGCCAACTAGCGGCTCTAGCCAAAcct                                                                                >  7:123947/1‑137 (MQ=255)                                                        gCAACCGTTTCATACAAAAAAAGCCCCCCGCTTCTTTTTGGCCGCCACC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCCCAGCCCCAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  3:398258/139‑1 (MQ=255)                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑C‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  1:1625/139‑1 (MQ=255)                                                        gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                                              <  8:320780/139‑1 (MQ=255)                                                          aaaCTTTTCATATAAAAACATCCCTCCGCTTCTTTTTGCCAGGCACCAT‑TTT‑T‑C‑CCAGCCCTTTTTCACCGCCGCAGCCCGAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                            <  8:345503/139‑1 (MQ=255)                                                           aaCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCGAGACAAACTCATGa                                                                            >  6:216563/1‑138 (MQ=255)                                                            ccttttCATCCAAAAACATCCCGCCGCTTCTTTTTGGCAGGACGC‑C‑TTT‑T‑ACGCACCCCGTTTTCGCAGCCCCACCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  7:197823/136‑1 (MQ=255)                                                            aCGTTTCATATAAAAAGATCCCGCCGCTGCTTTTTGGCAGGCAGCGT‑TTT‑T‑A‑GCCGCCCGTTTTCTCAGCCCCAGCCCGCAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  1:266343/139‑1 (MQ=255)                                                            aCGTTTCATATAAAAACATCCCTCCGTTTCTTTTTGGCAGGCAGCTT‑TTT‑T‑T‑CCCGCCGGTTTTCTCAGCCGCAGCGCGAAAAAAATGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  7:282559/139‑1 (MQ=255)                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          <  3:202184/139‑1 (MQ=255)                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          >  3:109979/1‑139 (MQ=255)                                                            aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTGTTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          >  2:90411/1‑139 (MQ=255)                                                             cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                                          >  5:347465/1‑138 (MQ=255)                                                                          aacaacCCTCCGGTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCTTTTTCGCCGCCCCCCCCCCAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGt                                                            <  8:136224/139‑1 (MQ=255)                                                                                    cGCTTCTTTTTGGCAGGCAGC‑CTTTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGaa                                                  >  5:398337/1‑139 (MQ=255)                                                                                         cttTTTGGCAGGCCCC‑C‑TTTAT‑A‑CCACCCCGTTTACTAACCCGCCCCCCGCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                              <  1:23235/138‑1 (MQ=255)                                                                                             ttgcctggccgC‑C‑TTT‑T‑ACGCCGCCGGTTTTCGCCGCCCCACCCCCCAACAAACCACGGCCGATTTCTTTTCTTGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                          <  3:378708/132‑1 (MQ=255)                                                                                             ttGCCAGGCCGC‑C‑TTT‑TAA‑TCAGCCCGTTTTCTCCGCGTCCCCCCGAAATAACCCACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                          <  4:108444/138‑1 (MQ=255)                                                                                             tgggcagccaGC‑GTTTT‑T‑A‑GCCGCCCGTTTTCTACGCCCCACCCCCCAACAAAACCCGCCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                          <  8:345811/136‑1 (MQ=255)                                                                                               ggcaggcaGCAT‑TTT‑T‑A‑GCAGCCCGTTTTATCCGCCGCACCCCGAAATAAAAGACGGCCTATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTca                                        <  4:222290/138‑1 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTTTGACAATATGAGTGAAAAAACCaa                   >  1:435171/1‑138 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc                  >  7:354598/1‑139 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc                  >  3:421628/1‑139 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc                  >  4:120136/1‑139 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc                  <  6:398337/139‑1 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACACATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATTAGAGTGAAAAAAGCaa                   >  4:43658/1‑138 (MQ=255)                                                                                                                        cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTGTTTGTAAAGATATTTGTGAAAAAACCAAc                  >  7:165177/1‑139 (MQ=255)                                                                                                                                        cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATTGAACGTGGAAAGAGAATTCATTTTCACAAAATGAGTGAAAAAACCAACTAGTTTTTAGAAGttt  >  8:54516/1‑139 (MQ=255)                                                                                                                                        cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTATAAGttt  >  2:178638/1‑139 (MQ=255)                                                                                                                                          |                                                                                                                                         CTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAGC‑C‑TTT‑T‑A‑GCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTT  >  CP006881/3725101‑3725369

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