Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,866,487	G→C	7.0%	intergenic (+22/+72)	ygbN → / ← rpoS	putative transporter YgbN/RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,866,487	0	G	C	7.0%	196.1 / 5.6	71	intergenic (+22/+72)	ygbN/rpoS	putative transporter YgbN/RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base C (2/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

TTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGC  >  NC_000913/2866347‑2866626                                                                                                                                             |                                                                                                                                              ttGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGAAAGACAAACAGGATAGGTTGAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGGATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGc                                                                                                                                             >  1:388145/1‑143 (MQ=255)      attatatatTTGTGGTTCTCGGTATCCGTCGTTCTGAAGTCCTGGACTGTGTTAACGACCATTCTCGGTTCTACTGGATTCTTAATTACCTGGTGCGTATGGGCTGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCgg                                                                                                                                        <  2:436167/138‑1 (MQ=255)        cAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAc                                                                                                                                      <  1:255803/143‑1 (MQ=255)         aagtatCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                     <  2:259275/140‑1 (MQ=255)         aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                     >  2:405738/1‑143 (MQ=255)         aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                     >  1:89172/1‑143 (MQ=255)         aaaTATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGAcc                                                                                                                                     >  1:891841/1‑143 (MQ=255)           ttatCTGGGGTTGGCTGTAGCAGACGGTCTGATAACCTGGACAGTGTTTACGGCCATTCTCGGTTTTACCGGATTCTTAATTACCTGGGGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGCCCATCGGACC‑tt                                                                                                                                  <  2:634765/142‑1 (MQ=255)           ggatCTGGGGTGGTGGGTAGCGGACGGGCTGGAAGCCTGGACGGGGGTAACGACCATTCTCGGGTGTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAAGAGAAAGGTCCGATGCCCATCGGACC‑tt                                                                                                                                  <  2:209753/141‑1 (MQ=255)                   ggTTTTCGGTAACAGCCGGTCTGCAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGGCCGATGGGCATCGGACC‑TTTTATtgtg                                                                                                                          <  1:651261/143‑1 (MQ=255)                   ggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgtg                                                                                                                          <  1:484113/143‑1 (MQ=255)                   ggTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATtgtg                                                                                                                          <  1:150321/143‑1 (MQ=255)                     ttGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGca                                                                                                                        <  2:529852/143‑1 (MQ=255)                      tGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGcac                                                                                                                       <  1:581718/143‑1 (MQ=255)                           gTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGGGGTAATTTGACCTTAATAAAACGGTCAGATGGGAACAGGAGA‑TTTTATTGTGAACAGaaa                                                                                                                  >  1:465128/1‑143 (MQ=255)                           gTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaa                                                                                                                  >  1:655081/1‑143 (MQ=255)                            tatgaggcGGTCTGAAAACCTGGAGAGTGGTAACGACGATTGTCGGGTTTACCGGATTCTTAATTACCTGGTGCGTGTGGGCGGTAAGTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGaaaa                                                                                                                 <  2:310616/136‑1 (MQ=255)                               cAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGc                                                                                                              >  2:12267/1‑143 (MQ=255)                                 gACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTCATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTCTATTGTGCACAGAAAAGGcaa                                                                                                            >  2:106858/1‑141 (MQ=255)                                      tCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTGAATAAAACCATCCGCTGGCCATCCGACC‑TTTTATTGTACACAGAAAAGGCCAGCCTc                                                                                                       <  1:175613/143‑1 (MQ=255)                                        tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGc                                                                                                     >  2:821776/1‑143 (MQ=255)                                        tGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATCGGCATCGGGCC‑TTGTATTGTGCACAGAAAAGGCCAGCCTCGc                                                                                                     >  1:10695/1‑143 (MQ=255)                                         gAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCgtg                                                                                                    >  2:537957/1‑141 (MQ=255)                                          aaaaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCAGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCTCAGAAAAGGCCAGCCCCGCtt                                                                                                   >  1:476065/1‑143 (MQ=255)                                            aaCCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTga                                                                                                 <  1:900279/143‑1 (MQ=255)                                                 ggACAGTGTTAACGATCATTCTCGGTTTTAGCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTg                                                                                            <  2:559781/143‑1 (MQ=255)                                                   aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGc                                                                                          >  1:700656/1‑143 (MQ=255)                                                   aCAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTAGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGc                                                                                          <  1:84201/143‑1 (MQ=255)                                                     aGTGTGAATGACTTTTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCt                                                                                        <  1:187900/143‑1 (MQ=255)                                                         ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCCCAGAAAAGGCCAGCCTCGCTTGAGACTCGCCTTCCt                                                                                    >  2:268773/1‑143 (MQ=255)                                                         ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGGAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                    >  1:545214/1‑143 (MQ=255)                                                         ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                    <  2:389678/143‑1 (MQ=255)                                                         ttAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                    <  2:389654/143‑1 (MQ=255)                                                         tcagcGACCAGTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGAGGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCt                                                                                    <  2:38882/139‑1 (MQ=255)                                                           aaCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGa                                                                                  >  1:330536/1‑143 (MQ=255)                                                            aCGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAATAAGGTCCGATGGGCATCGTACC‑TTTTATTCTAATCAGAATAGGTCATCGTCTCTTGAGACTGGCCTTTCTGAc                                                                                 <  1:389177/143‑1 (MQ=255)                                                             cGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTAGCTTGAGACTGGCCTTCCTGACa                                                                                >  1:616902/1‑143 (MQ=255)                                                                 cATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGAGAGATg                                                                            >  1:46157/1‑143 (MQ=255)                                                                     ctcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctta                                                                        >  1:412435/1‑143 (MQ=255)                                                                     ctcGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGctta                                                                        <  1:330862/143‑1 (MQ=255)                                                                       gcgTTTTACCGGATTCTACATTCCCCATTACCTCTAGGCGGTAGTGTGAGCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGcttact                                                                      <  1:762939/141‑1 (MQ=255)                                                                                     tCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTCACTTACTCGCGGAACac                                                        >  2:446829/1‑142 (MQ=255)                                                                                       ttAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAgcg                                                      <  2:805600/143‑1 (MQ=255)                                                                                         aaTTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCt                                                    <  1:41289/143‑1 (MQ=255)                                                                                          agtcgcTGGTGCGTATGGGCGGTAAGTTGACCTTAAGAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCtt                                                   <  2:594964/138‑1 (MQ=255)                                                                                           ttACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTc                                                  <  2:670779/143‑1 (MQ=255)                                                                                                  gTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                           >  1:298626/1‑143 (MQ=255)                                                                                                  gTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                           <  1:437534/143‑1 (MQ=255)                                                                                                  gTGCGAATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTc                                           <  1:121143/143‑1 (MQ=255)                                                                                                      gTATGGGCGGTAATTTGACCTTAATAAAAAGGGCCGATGGGCATCGGACC‑TTTTAGTGTGCACAGAAAGGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGcc                                       >  1:696692/1‑143 (MQ=255)                                                                                                               gtaaactGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCATCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGttt                             >  1:214830/7‑143 (MQ=255)                                                                                                               gTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGACATTCAGCCCC‑TGCGttt                             >  2:301164/1‑143 (MQ=255)                                                                                                                  atgtGACCGTAATAAGAAGGTCCGATCGGAGTAGGAAC‑CTATATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCa                          <  2:888623/140‑1 (MQ=255)                                                                                                                   tttGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAg                         <  2:720526/143‑1 (MQ=255)                                                                                                                     tGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGa                       <  2:465479/143‑1 (MQ=255)                                                                                                                         cTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTTGCGTTTCCAGGAttt                    >  2:347561/1‑143 (MQ=255)                                                                                                                          ttACTAAAAAGGTCCGGTGGGCAGCGGCCG‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTcg                  <  2:30940/143‑1 (MQ=255)                                                                                                                           taataaAAATGTCCGATTCTCATCGTACC‑TTTTATTGTGCACAGAAAAGTCCAGCCTCTCTTGAGACTTTCCTTTCTCACAGATGCTTACTTACTCGCGTAACAGCGCTTCTATATTCAGCCCC‑TGCGTTTGCAGGATTTcgc                 <  2:465128/143‑1 (MQ=255)                                                                                                                            aataaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTcgcg                >  1:720521/1‑143 (MQ=255)                                                                                                                              taaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGCACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCa              >  1:52579/1‑143 (MQ=255)                                                                                                                              taaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTGTGCAGGATTTCGCGCa              >  1:570660/1‑143 (MQ=255)                                                                                                                              taaAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCGTTCTGACCGATGCTTACTTCCTCGCGGGACAGCGCTTCGATATTCCGCCCC‑TACGTTTGCAGGCGTTCgcgcc              >  1:263379/1‑142 (MQ=255)                                                                                                                               aaaaaGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCaa             >  1:94507/1‑143 (MQ=255)                                                                                                                                aaaaGGGCCGATGGGCAGCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCaaa            <  1:596975/143‑1 (MQ=255)                                                                                                                                    ggTCCGATGCCCATCGGACCTTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACCGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCAAACgg         >  1:95790/1‑143 (MQ=255)                                                                                                                                     gTCCGATGCCCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCCGGATTTCGCGCAACCGgcg       >  2:852942/1‑143 (MQ=255)                                                                                                                                      tCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCAACCGgcgc      >  1:559311/1‑143 (MQ=255)                                                                                                                                       ccGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGC‑TACTTACTCGCGGAACAGCGCTTCGATATTCAGCACC‑TGCGGTTGCAGGATTTCGCGCAAACGGCGCAg    >  1:711751/1‑143 (MQ=255)                                                                                                                                        cGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCCTACCCACTCGCGGAACCGCGCTTCGATATTCCGCCCC‑TGCGTTTGCAACACCTCGCCCCAACGCCGCac    >  2:101279/1‑142 (MQ=255)                                                                                                                                          aTGGGCATCGGCCC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc  <  1:271298/143‑1 (MQ=255)                                                                                                                                          aTGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc  <  1:271170/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                              TTGTGACCAAATATCTGGGGTTGTCGGTAGCAGACGGTCTGAAAACCTGGACAGTGTTAACGACCATTCTCGGTTTTACCGGATTCTTAATTACCTGGTGCGTATGGGCGGTAATTTGACCTTAATAAAAAGGTCCGATGGGCATCGGACC‑TTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGC  >  NC_000913/2866347‑2866626

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