Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,756,047	T→C	intergenic (+92/‑112)	ssrA → / → intA	tmRNA/integrase IntA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,756,047	0	T	C	98.1%	171.5 / NA	52	intergenic (+92/‑112)	ssrA/intA	tmRNA/integrase IntA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (27/24);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTT  >  NC_000913/2755805‑2756285                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               gcctgGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCcg                                                                                                                                                                                                                                          <  1:229346/249‑1 (MQ=255)          tGAAGCGTTAAAACTTAATCAGGCTAGTTGGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCAAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGa                                                                                                                                                                                                                                <  2:294263/250‑1 (MQ=255)            aaGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGGCCGTCCCCAGCTGGCAAGCGAATGTAAAGACTGACTAAAAATGTAGTACCCAGGATGTAGGAATTTTGGACGCGGGTTCAACTCCCGCCCGCTCCCCCAAAATTCTCCAACGGTGATTACCAGAGTCATCCGATTAAGACCTAAGAGCCCGCACGGGGCAAGCCCCGCGGGCTTTTTTGTGCCCTCAACTCGTACCCCGAAGTCCGAagg                                                                                                                                                                                                                             >  1:684962/1‑250 (MQ=255)                gTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACGGACTAAGCATGGAGTACCGAGGATGTAGGAAGTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGTCTTTTGTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAAc                                                                                                                                                                                                                         <  1:641563/251‑1 (MQ=255)                       cTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGATGGCCAGAGAAAGTAAAGACTGACTAACAATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGGTCCACTCCCGCCCACCCCCCCAAAATTCTCCATCGGGGATTACCAGAGGCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTa                                                                                                                                                                                                                   >  1:475922/1‑250 (MQ=255)                           aTCAGGCTAGTTTGTTAGTGGCGGTTCCGTCCGCAGCTGGCAAGCAAATTTAAAGACTGACTAGGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATcc                                                                                                                                                                                                              <  1:521324/251‑1 (MQ=255)                                 cTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTAACGAGGATGTAGGAATTTCGGACGCGGGGTCAACTCCCCCCCACTCCACCCCAATTCTCCATCGGTGATTACCAGAATCATCCGATGAAGTCCTAAGAGCCCGCCCGGCGCAAACCCCGCGGGCTTTTTTGTGCCCCCAACTCGCCCCGCCACGCCCGCAGAGAACAAATTAAAACCGAACCt                                                                                                                                                                                                        >  1:875265/1‑251 (MQ=255)                                  tAGTTTGTTAGTGGCGTGTCCGTCCCCCGCTGGGAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCCAGGATGTAGGAATTTCGGACCCGGGGTCAACTCCCCCCAGCTCCACCAAAATTCTTCCTCGGTGAATACCAGAGACATCCGATGAAATCCTAAAAGCCCGCCCGGCGCAAGCCCTGCGGGCCTTTTTTTGGCCTCAACTTCTCCCCCGAAGTCCGACGAGAACTTATTAAATCCCAAccct                                                                                                                                                                                                       >  1:695167/1‑249 (MQ=255)                                              ggCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTg                                                                                                                                                                                           <  1:281058/251‑1 (MQ=255)                                                  tgtCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTAACGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGACGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGAt                                                                                                                                                                                         >  1:536098/1‑249 (MQ=255)                                                   gtccgtacGCCGCTGGAAAGCGAATGAAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCAAGCTCCACCAAAATTCTCCCTCGGTGATTACCAGAGTCATCCGCTGAAGTCCTAAGAGCCCGCACGGCGCCAGCCCTGAGGACTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTACTTAAATCCGAACCTTTTAGGCCCATTGCTAgg                                                                                                                                                                                      >  2:455371/1‑251 (MQ=255)                                                      cgtccgAAGCTGGCAATCGAATGTAAAGACTGACTAAGCATTCAGTCCCAAGGATGTAGGATTTTCGGACGAGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGACCCTTTTAGGCCCATTGATAGGccc                                                                                                                                                                                   <  2:684962/251‑1 (MQ=255)                                                           gCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACg                                                                                                                                                                               >  1:70021/1‑250 (MQ=255)                                                                     aGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTAt                                                                                                                                                                    >  1:333967/1‑251 (MQ=255)                                                                       cGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGACTTTCGGACGCGGGTTCAACTCCCGCTAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTa                                                                                                                                                                                                                   >  1:480814/1‑202 (MQ=255)                                                                       cGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCTAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCTCTCAACTCGTCCCGCGAAGTCCGGAGAGAACTAATTa                                                                                                                                                                                                                   <  2:480814/202‑1 (MQ=255)                                                                          aTGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGAAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGttt                                                                                                                                                               <  2:363729/251‑1 (MQ=255)                                                                              aaaGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGt                                                                                                                                                           >  2:56336/1‑251 (MQ=255)                                                                                aGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGACGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGATCTATTGTTTACGTTg                                                                                                                                                         >  2:169939/1‑251 (MQ=255)                                                                                     gactAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCt                                                                                                                                                    >  1:385508/1‑251 (MQ=255)                                                                                           gCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGc                                                                                                                                              >  2:238005/1‑251 (MQ=255)                                                                                              tgttgtaCCGAGGATGTAGGAATTTCGGACGCGGTATCAACTACCGACAGCTCCACCAAAATTCTCCATCGGTGAATAACAGAGTCATCCGATGAAGTCCTAAGAGGCCGCACGGCGCAAGCCCTGCGAGCTTTTTTGTGCAATCAACCCGTCCCGCGACGTCCGAAGCGCACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAAAGCAACGCTCTATTGTTTACGTTGGGCCTAAACGCAgg                                                                                                                                           >  1:614710/5‑251 (MQ=255)                                                                                                          gATGTAGTAATTTGGGACGCGGTTCCACCTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAAGCCGATCCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCAt                                                                                                                               <  2:858514/251‑1 (MQ=255)                                                                                                                      tCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCACGAAAAAc                                                                                                                   >  2:735630/1‑251 (MQ=255)                                                                                                                                  tCAACTCCCGCCAGCTCCACCAAAATTCGCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTGTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTaa                                                                                                        <  1:622597/250‑1 (MQ=255)                                                                                                                                    aaCTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTg                                                                                                     >  1:152920/1‑251 (MQ=255)                                                                                                                                        cccGCCATCTCCACCAAAATTCTCCATCGGTGATTGCCAGAGCCATCCGATGAATGCCTAATAGCCCGCACGGGGCAGGCCCTGCGGGATTTTTTGTGCCCGCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGACCCAACGAAAAGCTCTATTTTTTACTGTGTTCCTAAACGCAGGGAGTCTCCCCATGGCAAGAAAAACCAATCCGTTAACTGAt                                                                                                   <  2:62117/249‑1 (MQ=255)                                                                                                                                           gCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACgg                                                                                               >  1:229220/1‑250 (MQ=255)                                                                                                                                            ccAGCTCCACCAACATTCTCCAACGGCGATTACCCGAGACATCCGGTGAAGTCCTAAGAGCACGCCCGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGCAAAGCTCTATTGTTTACGTTGGGCATAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGaaaagga                                                                                              >  2:879374/1‑244 (MQ=255)                                                                                                                                             cAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACTCAGGGAGACTCCCCCTGGCAAGAAAAACCAAGCCGTTAAATGATACGGaa                                                                                             >  1:654152/1‑250 (MQ=255)                                                                                                                                                  ccaccaAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATGGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCa                                                                                         <  1:56336/249‑1 (MQ=255)                                                                                                                                                   caccaAAATTCTCCCTCGGTGATTACCAGATTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGATCCTTGTAGCCCCATTTATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAg                                                                                      <  1:879374/251‑1 (MQ=255)                                                                                                                                                          aTTCTCCATCGGTGATTCCCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCCCGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCaa                                                                               <  2:160700/251‑1 (MQ=255)                                                                                                                                                                 aTCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTaa                                                                         >  2:269647/1‑250 (MQ=255)                                                                                                                                                                       gATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGcc                                                                  <  1:933156/251‑1 (MQ=255)                                                                                                                                                                        aTTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTg                                                                                                                                                         >  1:6619/1‑163 (MQ=255)                                                                                                                                                                        aTTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTg                                                                                                                                                         <  2:6619/163‑1 (MQ=255)                                                                                                                                                                               gagTCATCCGATGAAGTCATCAGAGTCCGCACGGCGCCAGCACTGCGGGCTTTATTGTGCCCTCAACTCGTCCAGCGAAGTCCGACGAGAACTCATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCAAACGCAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATCCGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTaac                                                           >  2:877067/1‑248 (MQ=255)                                                                                                                                                                                agTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAACGATGCCGATTAcc                                                           >  2:738800/1‑249 (MQ=255)                                                                                                                                                                                   caccctaTGAAGGCCTAAGAGCCCGCACGGCGCAAGCCCGGCGGGCTTTTTTGGGCCCCCAACTCGCCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTACCGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCtt                                                      <  2:70021/245‑1 (MQ=255)                                                                                                                                                                                       cGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTAtg                                                  <  1:23295/251‑1 (MQ=255)                                                                                                                                                                                         aTGAAGTCCTAACAGCCCGCACGGCGCAAGCCCGGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATGAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCAGGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTAtgac                                                <  2:695167/251‑1 (MQ=255)                                                                                                                                                                                           gAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTAtgacgg                                              >  2:570768/1‑251 (MQ=255)                                                                                                                                                                                            aaGTCCTACGACCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTCtgacggg                                             >  1:714502/1‑250 (MQ=255)                                                                                                                                                                                                ccTAAGAGCCCGCCCGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGCCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCAATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTAtgacggtgacg                                         <  1:78619/251‑1 (MQ=255)                                                                                                                                                                                                 cTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGCACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTAtgacggtgccgg                                        >  2:79562/1‑251 (MQ=255)                                                                                                                                                                                                    agagCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCt                                     >  1:296221/1‑251 (MQ=255)                                                                                                                                                                                                     gagCCCGCACGGCGCAAGCCCGGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCtt                                    <  2:296221/251‑1 (MQ=255)                                                                                                                                                                                                            cACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGt                             <  1:269647/251‑1 (MQ=255)                                                                                                                                                                                                             aCGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGACCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGCTAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGtt                            >  1:307156/1‑251 (MQ=255)                                                                                                                                                                                                                      gCCCTGCGGGCTTTTTTGTGCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTc                   <  1:805177/251‑1 (MQ=255)                                                                                                                                                                                                                                       gtgCCCTCAACTCGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGACCCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGGGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCttt  <  1:53619/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               GCCTGGGGTTGAAGCGTTAAAACTTAATCAGGCTAGTTTGTTAGTGGCGTGTCCGTCCGCAGCTGGCAAGCGAATGTAAAGACTGACTAAGCATGTAGTACCGAGGATGTAGGAATTTCGGACGCGGGTTCAACTCCCGCCAGCTCCACCAAAATTCTCCATCGGTGATTACCAGAGTCATCCGATGAAGTCCTAAGAGCCCGCACGGCGCAAGCCCTGCGGGCTTTTTTGTGCCCTCAATTTGTCCCGCGAAGTCCGAAGAGAACTAATTAAATCCGAACCTTTTAGGCCCATTGATAGGCCCAACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAAACCTAAAGATGCCGATTACCAGCTTTATGACGGTGACGGGCTTACTCTGTTAATCAAGTCCAGTGGCAGTAAGCTTT  >  NC_000913/2755805‑2756285

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