Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	504,640	0	C	T	44.4%	57.6 / 31.5	64	K19K (AAG→AAA)	fsr	putative fosmidomycin efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (24/11);  new base T (0/28);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 4.12e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAA  >  NC_000913/504496‑504785                                                                                                                                                  |                                                                                                                                                   gAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑Tc                                                                                                                                             >  2:187817/1‑151 (MQ=255)        aTCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGAc                                                                                                                                       <  2:881420/150‑1 (MQ=255)        aTCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACg                                                                                                                                      >  2:348865/1‑151 (MQ=255)        aTCATGCCAATCGGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGACGTTCGGGC‑CTTGG‑TCGTTGAc                                                                                                                                       <  1:922761/150‑1 (MQ=255)             gCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGcc                                                                                                                                 <  2:461139/151‑1 (MQ=255)               cAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCcgc                                                                                                                               <  2:161221/151‑1 (MQ=255)                aaTCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCcgcg                                                                                                                              >  1:670175/1‑151 (MQ=255)                        tAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCcac                                                                                                                      >  2:69655/1‑151 (MQ=255)                            tGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGc                                                                                                                  >  2:203278/1‑151 (MQ=255)                            tGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGc                                                                                                                  >  2:480195/1‑151 (MQ=255)                             gTCAGAGAAAATTCTGACTGAAGCAGC‑GATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTg                                                                                                                >  2:612590/1‑151 (MQ=255)                                           tgattgaGGCAGCGAATAAACCGCCAGAATCAGCGTTTGGATCATGTCGTTCAGCAGAGGTGAGGGGTTAATACCACCAAAAATCCCAACCGAGGTTGGGGC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACt                                                                                                    <  1:622234/150‑1 (MQ=255)                                            gactgaAGCAGCGAAAAAATCGCAAAAATCAGCGTTTGGATCATGTCTTTCAGAAGATGTGAGGGCCTGATAGCACCAAAAATCCCAAACGATGTGGGGGC‑TTGGG‑TCGTGGCCGCACCCGCGCCCGCCACAGGCTGGGGTTGTTCACt                                                                                                    <  1:955007/149‑1 (MQ=255)                                              ctgaAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAtt                                                                                                >  2:843781/1‑151 (MQ=255)                                                      gggaaTAAATCGCCAAAATCAGCGATTGGATCTTGTCGTTCAGCAGATGTGAGGGGCTGATAGCACCAAAAATCCCAAACGTTGTGGGGGC‑TTTGG‑TCGTGGCCGCACCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAtt                                                                                                <  1:657833/139‑1 (MQ=255)                                                      gcGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAgg                                                                                        >  2:824656/1‑151 (MQ=255)                                                          aaaaaTCGCAAGAATAAGCGATTGGATCATGTCTTTCAGCAGATGGGAGGGGCTGAAACCACAAAAAATCCCAAACGAGGTTGGGGC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGa                                                                                       <  1:354215/146‑1 (MQ=255)                                                             aaccGCCAGAATCGGCGATTGGATCATGTCGTTCAGCAGATGGGAGAGGCTGATAGCACCAAAAAACCCAACCGATGTTCGGGC‑TTTGG‑TCTTTGCCGCAGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCa                                                                                 <  1:726385/148‑1 (MQ=255)                                                               ccGCCAGAATCGGCGTTTGGATCATGTCTTCCAGAAGATGTGAGGGGCTGATCGCACCAAAAATCCCAAACGTTGTGCGGGC‑TTTGG‑TCGTTGACCCAGCGGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACt                                                                               <  1:699364/150‑1 (MQ=255)                                                                  caaGAATCGGAGATTGGATCATGTCTTTCAGGAGATGGGAGGGGCTGATACCACCAAAAACCCCAAACTAGGTGCGGGT‑GTTGG‑TCTTGGACGCAGCGGCACCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAc                                                                                <  1:308411/145‑1 (MQ=255)                                                                          aGCGATGGTATCTTGTCTTTCAGCAGATGGGAGAGGCTGAAAGCACATAAAATCCCAAACGATGTGCGGGC‑CTGGG‑TCTTGGCCGCGCCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGt                                                                     <  2:659614/150‑1 (MQ=255)                                                                            cGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGc                                                                  >  1:203892/1‑151 (MQ=255)                                                                             gatggtatccttcagTTCAGCAAAGGTGGGAGGCGGATAGCACTTAAAATACCAAAGGTTTTCGGGGT‑TTTGG‑GTGTTGACGCACCCCCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGc                                                                  <  2:1010678/136‑1 (MQ=255)                                                                               ttGGTTAATGTCGTCCAGCAGATGGGAGAGGCTGATACCACCTAAAATCCCAAACGGGGTGCGGGC‑TTTGT‑TCTTTAAGGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAgg                                                                       <  1:410813/143‑1 (MQ=255)                                                                                tGGATCATGTCTTGCAGCAGGGGGGGGAGCTTAAAAGCCCCTAAAATCCCAAACGATGTCGGGGC‑CTTGG‑CCTTTAGCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGt                                                                     <  1:411960/144‑1 (MQ=255)                                                                                    tCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaa                                                           >  1:805592/1‑150 (MQ=255)                                                                                      attttgTTCGGCAAATGTGAGAGGCTGATACCACATAAAATCCCAATGGATTTTGGGGC‑CTTGG‑TCGTTGACGCCCCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaaa                                                          <  2:228461/144‑1 (MQ=255)                                                                                      aTGTCGTTCAGCAGAGGGGAGGGGCTGATAGCACCTAAAATCCCAAACGGTGTTGGGGC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAg                                                         <  1:360620/150‑1 (MQ=255)                                                                                       tGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGa                                                        >  2:727865/1‑150 (MQ=255)                                                                                              cagcagATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                                                >  1:983907/1‑151 (MQ=255)                                                                                              cagcagATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                                                >  1:256951/1‑151 (MQ=255)                                                                                                gcagATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt                                              >  2:889698/1‑151 (MQ=255)                                                                                                  agttgggaggGGTTGATAGCCCAAAAAATCCAAACCGATGTTCGGGC‑TTTGG‑TCGTTGACGCAGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAtt                                             <  1:344415/144‑1 (MQ=255)                                                                                                  agATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt                                            >  1:329881/1‑151 (MQ=255)                                                                                                  agATGTGAGAGGCTAAACCCCCAAAAAAAACAAAAGGATTTTGGGGC‑TTTGG‑TCGTGGACCCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATc                                                  <  1:954123/145‑1 (MQ=255)                                                                                                   gATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtg                                           >  2:745851/1‑151 (MQ=255)                                                                                                    aagggagaGGCTGATACCACCAAAAATACCAAACG‑GGGTCGGGCGCTGGG‑CTGTGGCCGGACCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                                                <  2:264485/141‑1 (MQ=255)                                                                                                     tgtgAGAGGCTGATAGCACAAAAAATACAAAACGATGTTCGGGC‑TTTGG‑TCTTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtg                                           <  2:557305/149‑1 (MQ=255)                                                                                                      gtgAGGGGCTGATAGCACAAAAAATACAAAACGATGTTCGGGC‑TTTGG‑TCGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtgt                                          <  2:561925/149‑1 (MQ=255)                                                                                                      gtgAGGGGCTGATACCCCAAAAAATCCCAAACGATGTTGGGGC‑TTTGG‑TCGTGGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATg                                                     <  1:753206/138‑1 (MQ=255)                                                                                                        gagaGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑TTTGG‑TCGTTGCCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat                                       <  2:882871/150‑1 (MQ=255)                                                                                                          gagGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattat                                    >  1:212002/1‑151 (MQ=255)                                                                                                            ggCTAAAAGCACAAAAAATACAAAACGATGTTGGGGC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta                                  <  2:782255/151‑1 (MQ=255)                                                                                                              ttGATAGCACCAAAAACACCAAACGATTTTGGGGT‑TTTGT‑TCTTTGACGCAGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt                                   <  2:250740/147‑1 (MQ=255)                                                                                                              cTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAcc                                >  1:62002/1‑151 (MQ=255)                                                                                                               tGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAcc                                >  1:867732/1‑150 (MQ=255)                                                                                                                 ataccgcaaaaaaCACAAATCGATGTTGGGGC‑TTTGG‑TCTTGGACGCACCGGCGTCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                                                <  2:655503/123‑1 (MQ=255)                                                                                                                 aTAGCACCAAAAAAACCAAACGATGTTGGGGC‑TTTGG‑TCGTTGCCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCa                               <  2:201308/149‑1 (MQ=255)                                                                                                                 aTACCACCAAAAATCCCAAACGATGTTGGGGC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAcc                                <  1:787337/148‑1 (MQ=255)                                                                                                                      aCCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAAc                        >  1:1011415/1‑151 (MQ=255)                                                                                                                       ccTAAAACACCAACGGATGTCCGGGC‑TTTGG‑TCGTTGACGCGGCGGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACa                       <  2:968916/151‑1 (MQ=255)                                                                                                                                 caaaccggtgtCGGGG‑CTTGGCTCTTTGACCAACCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTg                 <  1:605071/138‑1 (MQ=255)                                                                                                                                  aaaCGATGTTCGGGC‑TTTGG‑TGGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATg             <  2:326676/150‑1 (MQ=255)                                                                                                                                  aaaCGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGc            >  2:452988/1‑151 (MQ=255)                                                                                                                                    aCGATGTTGGGGC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgt          <  1:187795/151‑1 (MQ=255)                                                                                                                                       aTGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgtcgt       >  2:251630/1‑151 (MQ=255)                                                                                                                                          ttggggC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCt     <  2:549982/147‑1 (MQ=255)                                                                                                                                          ttCGGCC‑CTTGG‑TCTTTGCCGCGGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgt          <  2:173106/145‑1 (MQ=255)                                                                                                                                             gggC‑TTTGT‑TCTTTGACCCACCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgtcgt       <  1:728357/145‑1 (MQ=255)                                                                                                                                             gggC‑TTTGG‑TCTTGGCCGCACCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTgaa  <  2:620571/150‑1 (MQ=255)                                                                                                                                             gggC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCgtcgtc      <  1:435887/146‑1 (MQ=255)                                                                                                                                             gggC‑TTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTga   <  2:670246/149‑1 (MQ=255)                                                                                                                                             gggC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgt          <  2:262810/142‑1 (MQ=255)                                                                                                                                             gggC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTga   >  1:535092/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAACGATGTTCGGGC‑CTTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAA  >  NC_000913/504496‑504785

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