Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	3,725,237	A→G	intergenic (‑69/‑460)	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,237	0	A	G	78.3%	124.0 / 5.6	62	intergenic (‑69/‑460)	U712_19205/U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (34/13);  minor base C (3/10);  total (38/24)
Fisher's exact test for biased strand distribution p-value = 2.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAC  >  CP006881/3725119‑3725353                                                                                                                               |                                                                                                                     ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACCACCCTCCGCTTCTTTTTGGCAGGCAGGC‑A‑TTTT‑T‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCCCAAGACAAGCCGCCCGAtt                                                                                                   >  4:461684/1‑139 (MQ=255) ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                   >  4:117194/1‑139 (MQ=255) ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtt                                                                                                   >  5:236814/1‑139 (MQ=255)      caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                              >  1:396239/1‑139 (MQ=255)      caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                              >  6:400377/1‑139 (MQ=255)      caAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCGTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAgttcgtt                                                                                              >  1:402764/1‑139 (MQ=255)      caAAAACACCACCCACACCCGTAAACATACAAGCAAACG‑TTACATAGAAAAACACCCCTCCGCTTCTTTTTGGCCGGCA‑GC‑CTTTTT‑C‑CCAGCCCG‑TTTT‑CTCCCCCCCCCCCCGCAACAAAAGACGGCCGAtttcttt                                                                                              <  4:82895/139‑1 (MQ=255)           aGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg                                                                                         >  3:206701/1‑139 (MQ=255)           aGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg                                                                                         >  2:268793/1‑139 (MQ=255)             accaaccaCCCCCTTAAATAAACCAGCAACCG‑TTACATAGAAACACATCCCTCCGATTATTTTTGGCAGGCC‑GCAT‑TTTT‑A‑GCCGCCGT‑TTTT‑CTCCCCCGCACCCCTCAAAAAACGACGGCCGATTTCTTTTCTCGCa                                                                                       <  4:492794/139‑1 (MQ=255)                aTCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCGCGCAAAc                                                                                    >  7:2352/1‑139 (MQ=255)                 tCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAAc                                                                                    >  5:403082/1‑138 (MQ=255)                    acacCCTTAAAGAAACAAGCAACCG‑TTTAATACAAAAACATCCCTCCGTTTTTTTTTGCCCGCCA‑GCTT‑TTTT‑A‑CCAGCCCG‑TTTT‑CACCCCCCCCGCCCCAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                                <  4:193666/139‑1 (MQ=255)                      acCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCCCAAGAAAAGACGGCCGATTTCATTTATCGAAAACTATCgg                                                                              >  5:122788/1‑139 (MQ=255)                       cccGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCATCCGCAGCCCGCAAGATAAGACGGCCGATTTCTTTTCTCg                                                                                         >  8:262046/1‑127 (MQ=255)                       cccGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGATAAGACGGCCGATTTCTTTTCTCg                                                                                         <  7:262046/127‑1 (MQ=255)                           tAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGGA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCGTTTCTCGCAAAAGAGCGGcgg                                                                           >  4:299613/1‑135 (MQ=255)                           tAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCC‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAATTATTTTCTCGCAAACTAGCGGcgc                                                                           >  1:343922/1‑135 (MQ=255)                           tAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                           >  4:316512/1‑137 (MQ=255)                           tAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                           >  7:336418/1‑137 (MQ=255)                           tAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                         >  4:489286/1‑139 (MQ=255)                           tAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CGTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctc                                                                           >  5:94140/1‑137 (MQ=255)                              acaaaCAAGCAAACG‑TTCAAAACAAAAACAACCCCCCGCTTCTTTTTGGCCGCCC‑GC‑CTTTTT‑C‑GCCGCCGG‑TTTG‑CTCCGCCCCCCCCCGAAAAAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                      <  5:132830/135‑1 (MQ=255)                              acaTACAACCAAACG‑TTTCATATAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCCGCCCG‑TTTT‑CTCAGCCGCCGCCCGAAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                      <  5:123473/137‑1 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCC‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                   >  3:303367/1‑138 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCCGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                  >  6:197691/1‑139 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCCGCCCTCAAGAAAATACGGCCGAGTTGTTTTATAGGACACAAGCGGCAAGAGCAAAAc                                                                  >  1:435477/1‑139 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                  >  5:239357/1‑139 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGAAAGCAAAAc                                                                  >  5:266809/1‑139 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTATCGCAAACTAGCGGCTCTAGCaaaa                                                                   >  1:259418/1‑138 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTATTTTCTCGCAAACTAGCGGCTCTAGCAAAAc                                                                  >  6:123473/1‑139 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCGTTTCTCGCAAACTAGCGGCTCGAGCAAAAc                                                                  >  7:418780/1‑139 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGAA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCCCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                   >  7:294400/1‑138 (MQ=255)                                  aCAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCATTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGCCTGCCGATTTCTTTTCACGCAAACTAGCGGCTCGAGGAAAAc                                                                  >  1:383912/1‑139 (MQ=255)                                      caacaCG‑TTTCCTCGAAAAAACCCCCTCCCCTTCTTTTGGGCCGGCA‑GC‑CTTTTT‑A‑GCCGCCCG‑TTTA‑CCCCGCCGCCGCCCCCAATAAACCCCGGCCGATTTCTTTGCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                              <  6:130744/135‑1 (MQ=255)                                      caaaaCG‑TTTCAAAGAAAAACAACCCTCCGCTTCTTTTTGCCCGGCA‑CC‑CTTTTT‑A‑GCCGCCCG‑TTTT‑CTCCCCCCCAGCCCGCAACAAAACACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                              <  7:474695/137‑1 (MQ=255)                                       cAAACGTTTTAATATAAAAACACCCCGCGGTTTTTTTTGGGCAGGCC‑CC‑CTTTTT‑A‑GCCGCCGG‑TTTG‑CCCCCCCCCCGCCCCAAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCAt                                                              <  4:206701/139‑1 (MQ=255)                                        aaaCT‑TTTCTAACAAACACAGCCCTCCGCTTCTTTTGGGCAGCCA‑GC‑C‑TTTT‑ACGCACCCCG‑TATC‑CTCCGCCGCAGCCCGCAACAAAACACGGCCGTTTTCTTTTCTGGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                            <  6:198500/139‑1 (MQ=255)                                        aaaCG‑TTTAATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TTTT‑ACGCAGCCCG‑TTTT‑CTCAGCCCCCGCCCGCAAAAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                            <  1:268793/139‑1 (MQ=255)                                          gcG‑TTTCATACAAAACCAGCCCGCGGTTTGTTTTTGGCGGCCC‑GCGT‑TTTT‑T‑CCCCCCCG‑TTTT‑CTCCGCCCCCCCCCGAAACAAATCCCGGCCGATTTCTTTTCTCGCAAACAAGCGGCTCTAGCAAAACTCATGAct                                                          <  2:343922/138‑1 (MQ=255)                                          aCG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  3:489286/139‑1 (MQ=255)                                          aCG‑TTTCATACAAAACCACCCCTCCGGTTCTGTTGGGACGGCC‑GC‑C‑TTTT‑ACCCCGCCCG‑TTTC‑CTCTCCCCCCGCCCGCAACAAACGACGGCCGATGTCTTTTCTCGCAACCTAGCGGCTCTAGCAAAACTCATGAct                                                          <  6:94140/139‑1 (MQ=255)                                          aCG‑TTTCAAACAAAAACACCCCTCCGCTTATTTTTGGCCGCCC‑GC‑CTTTTT‑A‑GCA‑CCCG‑TTTTGCCCCGCCGCACCCCCAAACAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  5:235015/139‑1 (MQ=255)                                          aCG‑TTTAATAGAAAAACATCCCGCCGCTTGTTTTGGCCCGGCA‑GCTT‑TTTT‑C‑CCACCCCG‑TGTC‑CGACGCCCAACCCCCCAAGAAACCACGGCCGATTTCTTTTCTCGCAAACAAGCGGCTCTAGCAAAACTCATGAct                                                          <  8:29279/139‑1 (MQ=255)                                          aCG‑TTTAATACAAAAACATCCCTCGGCTTATTTTTGGCCGGCC‑GCGT‑TTTT‑A‑GCCGCCCG‑TTTT‑CCCCGCCCCCCCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  7:381092/139‑1 (MQ=255)                                          aCG‑ATTAATATAAAAACATCCCTCCGATTCTTTTTGGCCGCCA‑GC‑CTTTTT‑A‑CCCGCCCG‑TTTT‑CCCAGCCGCAGCCCCCAACAAACCACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  2:252293/139‑1 (MQ=255)                                           actttacATACAAAAACATCCCTCCGCTTCTTTTTGCCCGGCA‑GC‑CTTTTT‑C‑GCCCCCCG‑TTTT‑CTCAGCCCCAGCCCGAAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAct                                                          <  3:387060/133‑1 (MQ=255)                                                  acacctaaaCATCCCTCGGATTTTTTTTGGAAGGCA‑GC‑CTTTTT‑C‑GCAGCCGG‑TTTA‑CTCAGCGGCCCCCCGCAATAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAt                                                   <  6:468488/133‑1 (MQ=255)                                                     caaaaaCACCCCCCCGCTCTTTTTTGGCCGGCC‑GC‑T‑TTTT‑T‑ACCCCCCGTTTTT‑CACCGCCCCCCCCCCCAACAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGa                                                <  5:400377/138‑1 (MQ=255)                                                       aaaaCATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                                              >  2:40414/1‑139 (MQ=255)                                                           cATCCCTCCGTTTCTTTTTGCCCGTCT‑GCTT‑TTTT‑C‑GCCGCCCT‑TTTT‑CTCGGCCGCAGCCCCAAAAAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg                                          <  1:151304/139‑1 (MQ=255)                                                              cccTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TTTTAA‑GCCGCCGG‑TTTT‑CACAGCCGCACCCCCCAATAAACGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                                       <  7:16261/139‑1 (MQ=255)                                                              cccTCCGCGTCTTTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAGTTCTTTTCTCGCAAACGAGCCGCTCGATGAAAACTCATGCCTCTATGATACTCCGTggg                                         >  3:4134/1‑136 (MQ=255)                                                                     cttcttTTTTGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAACACGGCCGATTTATTTTATCGCAAACTAGCGGCTCTATCAAAACAAATGACTATATGAAATACCGTGGAAAGAGAAt                                 >  2:151304/1‑138 (MQ=255)                                                                        cttTTTGGCAGGCA‑GC‑CTTTTT‑A‑GCAGCCCG‑TTTT‑CTCACCCGCAGCCCGCAAGAAAAGACGGCCGAATTCTTTTATCCCAAAATAGCGCCTCTAGCAAAACACAAGACTAAATGATAGACACTCGAAAGaaaaataat                             >  8:410023/1‑134 (MQ=255)                                                                                                     agcCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTAGCAAACAAGCGGCTATAGCAAAACTCATGACTCTAAGATATACCGTGGAAAGAGAAATCTTTGTTACAATAAGAGTTAAAAAAca     >  4:279758/1‑137 (MQ=255)                                                                                                      gcCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGActc                                                         <  7:237469/85‑1 (MQ=255)                                                                                                      gcCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGActc                                                         >  8:237469/1‑85 (MQ=255)                                                                                                      gcCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa   >  3:118877/1‑139 (MQ=255)                                                                                                      gcCCG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACAAATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTAACAATATGAGTGAAAAAAACaa   >  3:227814/1‑139 (MQ=255)                                                                                                       cccG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGGCCCAATATGAGTGAAAAAACCAAc  >  1:439536/1‑139 (MQ=255)                                                                                                       cccG‑TTTT‑CTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAGAGACCGTGGAAAGAGAATTCTTTGTCACAATATGATTGAAAAAACCaa   >  7:429926/1‑138 (MQ=255)                                                                                                                               |                                                                                                                     TTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACG‑TTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TTTT‑A‑GCAGCCCG‑TTTT‑CTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAC  >  CP006881/3725119‑3725353

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