Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,189,573	A→C	5.9%	*1408Y (TAA→TAC)	rpoC →	RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,189,573	0	A	C	5.9%	174.7 / 2.1	69	*1408Y (TAA→TAC)	rpoC	RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base C (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

CCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAAC  >  NC_000913/4189431‑4189713                                                                                                                                               |                                                                                                                                              ccGTTTACGCGTACCACCAGGTTCGTATGCGTCGCCGTGCTGCGGGTGATGCTCCGGTTGCACTTTAGGTGTCTTCAGAAGACGCATCTGCCAGCCTGGCAGAACTGGTGAACGTAGGTCTGGGCGGTTCTGATAACGAGTaa                                                                                                                                               <  1:62175/143‑1 (MQ=255) ccGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCCGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATACCGAGTaa                                                                                                                                               >  1:292712/1‑143 (MQ=255)   ggTGACGGGTACAACCAGGAACGGAGGCGGCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCAGCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATc                                                                                                                                             <  1:144556/143‑1 (MQ=255)    ggtccgcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCg                                                                                                                                            <  2:486392/139‑1 (MQ=255)     ttACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGt                                                                                                                                           >  1:596732/1‑143 (MQ=255)       aCGCGTAGCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGGAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTa                                                                                                                                         <  2:384516/143‑1 (MQ=255)         gcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAAt                                                                                                                                       <  1:96885/143‑1 (MQ=255)              ccaccgGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa                                                                                                                                  <  2:602823/143‑1 (MQ=255)                  caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGCTAACGAGTAATCGTTACTCCGCaaacc                                                                                                                              >  1:358622/1‑141 (MQ=255)                  caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATa                                                                                                                              >  2:836635/1‑143 (MQ=255)                      aTCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGt                                                                                                                          <  1:794583/143‑1 (MQ=255)                             gCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAcc                                                                                                                   >  2:864548/1‑143 (MQ=255)                               gTCGCCGTGCTGCGGCGGACGCTCCGGCTGCAGCGCAGGTGACTGCAGAAGACGCAGCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCg                                                                                                                 <  2:505040/143‑1 (MQ=255)                                  gCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGCAGACGCATCCGCACGGAGAGTATATCTGCTGAACGCAGGTCTGGGCGGTTCTGCCAACGCGCAATCGGCAATCCGCACAGAACGTACAAACCCGcgt                                                                                                              >  2:826825/1‑141 (MQ=255)                                           cGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑tt                                                                                                    >  2:369192/1‑143 (MQ=255)                                             ggTGAAGTTCCGGCTGCACCGCAGGTGTTTGTGGGTGACCAAACCACCAGCCTGGCAGAACTGCTGAACGCAGGTCGGGGCGGGTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑tttt                                                                                                  <  1:371697/143‑1 (MQ=255)                                             ggTGAAGCTCCGGCTGCACGGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑tttt                                                                                                  <  2:310094/143‑1 (MQ=255)                                                gAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑ttttttt                                                                                               >  1:871842/1‑143 (MQ=255)                                                  aGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGACCTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTGATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTAt                                                                                             >  1:286600/1‑143 (MQ=255)                                                   gCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGGTAACGAGTAATGGTTAATCCGCAAATAATGTAAAAAACAGATTCAGCGGG‑TTGTTTGATg                                                                                            >  2:439494/1‑143 (MQ=255)                                                   gCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATg                                                                                            <  1:784094/143‑1 (MQ=255)                                                    cTCCGGCTGTACCGCAGGTGACTGCAGAAGACGCAGCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATgg                                                                                           <  1:876774/143‑1 (MQ=255)                                                    cTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATgg                                                                                           <  2:602432/143‑1 (MQ=255)                                                      ccGTCTGCAGCGCAGGTGACGGCAGAAGACGCATCTGCCAGCCAGGCAGAACTGCGGAACGCAGGTCTGGGCGGTTCTGATAACGAGTATTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATgggg                                                                                         <  1:235260/143‑1 (MQ=255)                                                      ccGGCTGCGCCGGAGGTCACTGCAGAAGCCGCCGCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAGCGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATgggg                                                                                         <  2:811212/143‑1 (MQ=255)                                                      ccGGCTGCACCGCAGGTGACTGCAGAAGACGCAGCGGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAAACGCTTCGGCGGG‑TTTTTTTATgggg                                                                                         <  2:784142/143‑1 (MQ=255)                                                             cACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGCTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑CCTTTTTATGGGGGGAGtat                                                                                  >  2:882258/1‑141 (MQ=255)                                                               ccGCAGGTGACTGCAGAAGACGCATCTGCCAGTCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAg                                                                                <  1:247995/143‑1 (MQ=255)                                                               ccGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATGACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGtgttg                                                                                >  2:164611/1‑139 (MQ=255)                                                               ccGCAGGGGACTGCGGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCGTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAg                                                                                <  1:738919/143‑1 (MQ=255)                                                                cGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTCAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTCAgg                                                                               >  1:297301/1‑143 (MQ=255)                                                                      tGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAaga                                                                         <  1:864548/143‑1 (MQ=255)                                                                       gACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCAGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTATAGGGAAagag                                                                        >  2:298559/1‑143 (MQ=255)                                                                            cAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAttt                                                                   <  1:667640/143‑1 (MQ=255)                                                                            cAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTGTTACGGGGGGCGTGTAGGGCAAGAGCAttt                                                                   >  2:220817/1‑143 (MQ=255)                                                                             agaagaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTg                                                                  >  2:646299/1‑143 (MQ=255)                                                                              gaagaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGt                                                                 <  1:119995/143‑1 (MQ=255)                                                                                 gaCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGCTTCTTTTGTGGGGCGGGCTCAGGGAAAGCCCATCTGTCa                                                               >  2:674230/1‑143 (MQ=255)                                                                                   cGCATCTGCCAGTCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGGAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGaa                                                            <  1:365025/143‑1 (MQ=255)                                                                                    gCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGTTTAATAGCATTTGTCAGact                                                           >  1:839220/1‑141 (MQ=255)                                                                                         tGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGGTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTCGTCCGAATGTTTa                                                      >  2:858704/1‑143 (MQ=255)                                                                                          gtcAGCGTGGCCGGAGTGCTGAACGGAGGTCTGGGTGGTGCTGATAACGAGTAATCGTTAATCCGCAAATATTGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGGTTAGAGAAAGAGCATTTGTCAGAATATTTaa                                                     <  2:679055/141‑1 (MQ=255)                                                                                          gcaaGCCTGGCAGAACTGCTGAACGCAGGTCGGGGCGGTTCTGATAACGAGGAATCGTTAATCCGCAAATAACGGAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                     <  1:729144/140‑1 (MQ=255)                                                                                          gCCGGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCGTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                     <  2:516558/143‑1 (MQ=255)                                                                                          gCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATGCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTaa                                                     <  1:447838/143‑1 (MQ=255)                                                                                             aGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑GTTTTTTATGGGGGGAGTTTAGGGAACGAGCATCTGTCAGAATATTTAAGGa                                                  >  2:335497/1‑143 (MQ=255)                                                                                                  ggCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGAAGAGAATAGGGAAAGAGCATTTGTAAGAAAAATTAAAGAAGTTc                                             >  1:368842/1‑143 (MQ=255)                                                                                                   gCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCt                                            <  2:873757/143‑1 (MQ=255)                                                                                                     aGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGa                                          >  1:573591/1‑143 (MQ=255)                                                                                                      gAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCTGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGaa                                         <  1:787084/143‑1 (MQ=255)                                                                                                           gctgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTc                                    >  2:21037/1‑143 (MQ=255)                                                                                                             tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGGGTTTAGGGAAAGAGCAGTTGGCAGAATATCTAAGGAATTTCTGGATAATCAt                                  >  2:270599/1‑143 (MQ=255)                                                                                                              gAACGAAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATa                                 >  2:741920/1‑143 (MQ=255)                                                                                                               aaCGCAGGTCTGGGCGGTTCTGATGACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATaa                                <  1:398390/143‑1 (MQ=255)                                                                                                                    aGGTCTGGCCGGTTCTGATAACGAGTAACCGTTAATCCGCAAATAACGCAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAAt                           <  1:164611/143‑1 (MQ=255)                                                                                                                        cTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGtcg                       >  1:726122/1‑141 (MQ=255)                                                                                                                          gggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTagag                     >  1:775799/1‑143 (MQ=255)                                                                                                                           ggCGGTTCTGATAACGAGTAATCGTTAATCCGCAAAGAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAGTTGTGAGAATATTGAAGGCATTTCTGAATACTCATAATCACTGTAgaga                    >  2:609901/1‑143 (MQ=255)                                                                                                                            gCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGAt                   <  1:836635/143‑1 (MQ=255)                                                                                                                                ttCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACt               >  2:99301/1‑143 (MQ=255)                                                                                                                                 tttGATAACGAGTATTCGTTAATCCGCAAAGAACGTAAAGTCCCGCGTCGGCGGG‑TTTATTTATTGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTa              <  2:535956/141‑1 (MQ=255)                                                                                                                                  ctcattaCGAGTGATGGTGAATCCGCACATAACGTAAAAACGCGCTTCGGCTCG‑TTTTTTTCTGGTGGGTATAGAGGTAAAGGGCACTTGTCACAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTaa             <  2:315573/137‑1 (MQ=255)                                                                                                                                    attAACGATTACTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCGTAATCAAGTTAGAGCTGACTAata           >  2:49370/3‑143 (MQ=255)                                                                                                                                     ttAACGATTACTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGttt                                                                                  <  1:179307/70‑1 (MQ=37)                                                                                                                                     ttAACGATTACTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGttt                                                                                  >  2:179307/2‑71 (MQ=37)                                                                                                                                      tAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAA‑GAAATTCTGAATACTCATAATCAATGTAGAGATGACTAATATcc        >  2:373623/1‑143 (MQ=255)                                                                                                                                       aaCGAGTACTCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATcc        <  2:847092/143‑1 (MQ=255)                                                                                                                                          gAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGACTACTCATAATCATTGTTTAGTTGACTTGTATCTTGa     >  1:106311/1‑143 (MQ=255)                                                                                                                                          gAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTATAAATGACTAATATCCtaa     >  2:146344/1‑141 (MQ=255)                                                                                                                                             tAATCGTTAATCCGGAAATAACGTTAAAACCCGCTTTGGGGGGTTTTTTTTGTGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGaaa   <  1:674230/143‑1 (MQ=255)                                                                                                                                             tAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACAAATATCCTGAAac  >  2:749067/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                              CCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAGACGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAGAATATTTAAGGAATTTCTGAATACTCATAATCAATGTAGAGATGACTAATATCCTGAAAC  >  NC_000913/4189431‑4189713

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