Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,311,650	A→G	F367S (TTC→TCC)	ompC ←	outer membrane porin C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,311,650	0	A	G	95.5%	146.7 / ‑6.2	44	F367S (TTC→TCC)	ompC	outer membrane porin C
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (23/19);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGG  >  NC_000913/2311418‑2311892                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   gCCTTCAAGAACGGTAACTTCAGTAGCTTTGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATGAGGACTGTTAAACCAGAccc                                                                                                                                                                                                                                   <  2:452147/250‑1 (MQ=255)      cAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGGTCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGGTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGAGAT‑CAAGCGAGAGTAGGACTGGTAAACCAGACCCAGAGc                                                                                                                                                                                                                              <  2:55125/250‑1 (MQ=255)        aGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAAAGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCCGACCCAGAGCTAc                                                                                                                                                                                                                           >  2:90989/1‑251 (MQ=255)          aaCGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAAGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAAACCCAGAGCTACg                                                                                                                                                                                                                          >  2:9744/1‑250 (MQ=255)           aCGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGAt                                                                                                                                                                                                                        <  1:328760/251‑1 (MQ=255)             ggTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGGAAAGCTTATTTCCATACGCGATTATCCTCATGCGAACGGTCGCAAGAGGACACCAAAAAACTGAGTTTTTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTTTTCGCTAT‑CAATCGAGATTAGGACTGGGAAACCAGACCCAGAGCTACGATGt                                                                                                                                                                                                                      >  2:342309/1‑251 (MQ=255)              gTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCTCGCGGGTAAAGCTTAGTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGGACGCTGAAAACAATGAAAAAAGTGCCCGCAGGCCCTTTGTTCGAGAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGtt                                                                                                                                                                                                                     <  1:250639/251‑1 (MQ=255)              gTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGtt                                                                                                                                                                                                                     <  1:177925/251‑1 (MQ=255)                         aGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCCAAAAAGACGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCCTGCGAACGGTCGCAAGAGTACACCAAAAAAATGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGGTCT‑CAATCGAGATTAGGACCTGCAAACCAGACCCAGAGCTACGATGTTCTCAGGGtttg                                                                                                                                                                                                          >  2:159120/1‑249 (MQ=255)                              gggTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCCGCAGTGCCACGCGGGGAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGcc                                                                                                                                                                                                      >  2:106970/1‑250 (MQ=255)                                     ggCGGGTTGGTGTTTTTGATCGCAACCAACAAAGAAGAGCGCGTCAGCCAGCCGTGCCTCGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAATCTGAGTTTGTCCGCTGAAAACAATGAAAAAAGTGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGGAAACCAGACCCAGATCTACGATGGTATCTGTGTTGATGCCAGCGTCAc                                                                                                                                                                                              <  2:346689/251‑1 (MQ=255)                                          gTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCGGTGTTCGATAT‑CAATCGAGAGTAGGACTGGTAAACCAGACCCATAGCCACGATGTTATCAGTGGTGATGCCAGCGTCACGAGTg                                                                                                                                                                                         <  1:171421/251‑1 (MQ=255)                                                      gATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                                                                                                                                                             >  1:198801/1‑251 (MQ=255)                                                        tCGCAACCAACAAAGAAGGGCGCGGCCGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCACGAGTACACCACAAAACTGAGTTTGTCCGCTGAAAACACTGACAAAAGGGACCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTAAGATGTTATCAGTGTTGATGCCAGCGTAACGAGTGAACTGGTTgtcgtc                                                                                                                                                                           >  2:328503/1‑251 (MQ=255)                                                                                                cgggggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAATCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCTCTTTGTTCGATAT‑CAATCGGGATTAGGACTGGTAAACCAGACCCAGTGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAGCTGGTTGTCGTTCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttt                                                                                                                                   <  1:469478/248‑1 (MQ=255)                                                                                                 gcgGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGtttt                                                                                                                                  <  1:1945/251‑1 (MQ=255)                                                                                                   gggTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttg                                                                                                                                <  1:9744/251‑1 (MQ=255)                                                                                                           cTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAagt                                                                                                                         >  1:365241/1‑250 (MQ=255)                                                                                                                  ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAgtaggt                                                                                                                  >  2:145325/1‑250 (MQ=255)                                                                                                                      aCGTGATTATCCTCATGCGAACGGCCGCCAGAGTACACCACAAACCTGAGTTTGCACGCTGAAAACAATGAAAAACCGGCACGCAGGCCCTTTGTTCGATAT‑CAATCGCGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGCGTGACCTGGTTGTCGTCCAGCCGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCAc                                                                                                             >  2:306403/1‑251 (MQ=255)                                                                                                                        gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACca                                                                                                           <  2:329720/251‑1 (MQ=255)                                                                                                                                    aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCAAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACatat                                                                                               >  2:329233/1‑251 (MQ=255)                                                                                                                                            ggTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGat                                                                                       <  1:392296/251‑1 (MQ=255)                                                                                                                                               cGCAAGAGTACACCAAAAAACTGAGTTGGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGGCCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGGTGATTTTGTAGTCAACGTAGGTGGACATGTTGTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGatat                                                                                     <  1:76324/250‑1 (MQ=255)                                                                                                                                                gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCt                                                                                   <  2:410997/251‑1 (MQ=255)                                                                                                                                                 cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGGCCAGCAGGTTGATTTTGTAATCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCtt                                                                                  <  1:419926/251‑1 (MQ=255)                                                                                                                                                  aaGAGTACACCAACAACCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCCCGAGTGAACTGGTTGTCGTCCCGCAGGTTGATTTTGTAGTCAACGTAGGGGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCACCATATTTCAGGATATCTtc                                                                                 >  2:446238/1‑251 (MQ=255)                                                                                                                                                     agTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTtcgtc                                                                              >  2:479625/1‑251 (MQ=255)                                                                                                                                                             aaaaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATAGCTTCGTCGTCGTAGc                                                                      <  1:328503/251‑1 (MQ=255)                                                                                                                                                                 aaCTGAGTTTGTACGCTGAAGACAATAAAACAAGGGACAGCAGGCCCTTTGTTAGATATACAA‑CGAGCTTAGAACTGGTAAACCAGACACAGAGAAACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGCACTGGTTGTCGTCCAGCCGGTTGATTTTGAAGTCAACGTAGGTGGGCATGTTTTTGTTGCAGTAGAAGGTAGCACCAACATCACCATATTTCAGGATATCTTCGTCGTCGTAGCCACg                                                                  >  2:313495/1‑251 (MQ=255)                                                                                                                                                                    tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTTAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAcc                                                               >  1:437271/1‑251 (MQ=255)                                                                                                                                                                           gTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGttt                                                        >  2:223323/1‑251 (MQ=255)                                                                                                                                                                           gTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGttt                                                        <  2:218560/251‑1 (MQ=255)                                                                                                                                                                                       cAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGAc                                            <  1:223323/251‑1 (MQ=255)                                                                                                                                                                                       cAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGAc                                            >  2:57825/1‑251 (MQ=255)                                                                                                                                                                                               aaaGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGt                                     >  2:315502/1‑250 (MQ=255)                                                                                                                                                                                                aaGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTa                                    >  2:305084/1‑250 (MQ=255)                                                                                                                                                                                                  gggCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTATTGAAGTAGTCGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGc                                 >  1:435252/1‑251 (MQ=255)                                                                                                                                                                                                        gCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAggg                            >  2:162837/1‑250 (MQ=255)                                                                                                                                                                                                          aggccccttgcccGATAC‑CAATCAAGCTTAGAACCGGTCAACAAGCCCCCGAGCTACGATGCTATCAGTGTTGATGCCAGCGCCAAGCGAGAACTGGTTATCGTCCAGCCGGTTGATTTTGCAGTCCACGACGGTGGAAATGTTTTTGTTGAAGTAGTAGGTAGCACCAAACTACACATATTTCAGGCTCTATTCGTCGTCGTAGCCACGCCCCAGGTTTTTACCTTTCGACTGCAGGTAAGCCAGggacg                         >  2:200292/13‑251 (MQ=255)                                                                                                                                                                                                          aGGCCCTTTGTTCGATAT‑CAATCGAGATTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGCACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCACCGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGAAAGCCAGggacg                         >  2:361577/1‑251 (MQ=255)                                                                                                                                                                                                                tttGTTCGATAT‑CACTTGAGATTAGGACTGGTCAACAAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGACCTGGTTGTCGTCCAGCAGGTTGCTTTTGTAGTCCACGTAGGTGGACATGTTTTTGTTGCCGTAGTAGGTAGCACCAACATCAACATATTTCAGGCTATCTTCGTCGTCGTAGCCCCGCCCCAGGTTTTTAACTTTAGCCTGCAGGTa                                    >  2:212611/1‑234 (MQ=255)                                                                                                                                                                                                                tttGGTCGATAT‑CAATCGAGATTAGGAGTGGTAAACCAGACCCAGAGCTACGATGTGATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCGGGTTTATTTTGTAGTCAACGTAGGGGGACATGTTTGTTTGGAATGAGTAGGTAGCGCCAACATCAACAGATTTCAGGATGGTTTTTTCGTCGAAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTa                                    <  1:212611/234‑1 (MQ=255)                                                                                                                                                                                                                                  gagaTTAGGACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGCAAGCCAGGGACGGACGCAGACCGAAGTCGAACTgg  >  2:200575/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   GCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATAT‑CAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGG  >  NC_000913/2311418‑2311892

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