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

GGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAG  >  NC_000913/504498‑504786                                                                                                                                                |                                                                                                                                                     gggTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGt                                                                                                                                            <  2:615384/151‑1 (MQ=255) gggTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGt                                                                                                                                            <  1:779122/151‑1 (MQ=255)    ttATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTg                                                                                                                                          >  2:471001/1‑150 (MQ=255)    ttATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTg                                                                                                                                          >  1:240579/1‑150 (MQ=255)    ttATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTg                                                                                                                                          >  2:802981/1‑150 (MQ=255)            ccAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCcg                                                                                                                                 >  1:376776/1‑151 (MQ=255)               aTCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCgcgc                                                                                                                              >  2:618333/1‑151 (MQ=255)                tCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGcgcc                                                                                                                             >  2:961068/1‑151 (MQ=255)                    cATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGcgcccgcc                                                                                                                         >  2:108711/1‑151 (MQ=255)                     aTAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGcgcccgcc                                                                                                                         >  1:226574/1‑150 (MQ=255)                                 gaAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgtt                                                                                                            >  1:695378/1‑151 (MQ=255)                                  aaaaTTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGgttg                                                                                                           >  1:306157/1‑151 (MQ=255)                                      ttCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTc                                                                                                        >  1:590774/1‑150 (MQ=255)                                       tCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCAc                                                                                                      >  2:840074/1‑151 (MQ=255)                                       tCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCAc                                                                                                      >  2:855795/1‑151 (MQ=255)                                           actgaAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAt                                                                                                  >  1:874779/1‑151 (MQ=255)                                           acggaAGCAGCGAATAAATCGCAAGAAGCGGCTATTGGTTCATGTCGTTCAGCAGATGTGAGGGGGTAATAGCACCAAAAACACCAAA‑CGATGTGCGGGTT‑TTGG‑TCGTGGACCCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAt                                                                                                  <  2:1043705/148‑1 (MQ=255)                                            cggaAGCGGCGGATAAATCGCAAGAATCAGCGTTTGGATCATGTCGTCCGGCAGATGTGGGGGGCTGAAAGCACCTAAAATCCCAAA‑CGATGTGGGGTCT‑TTGG‑TCGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCAtt                                                                                                 <  2:188111/149‑1 (MQ=255)                                               agggggggaaaaaaCCGCCAAAATCAGCGTTTGGATCAGGTCTTTCGGCAGTTGTGGGGGGTTGATAGC‑CCAAAAAATCCAAA‑CGGTTGTGGGGCCTTTGG‑TCGTTGCCGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTg                                                                                                <  2:404105/138‑1 (MQ=255)                                               aaccaggGGATAAACCGCAAAAATCAGCGTTTGGATCATGTCGTTAAGAAGATGTGAGAGGCTGATACCCCTAAAAATCCCAAA‑CGAGGTTCGGTCT‑TTGG‑TCTTTGCCGCCGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCAtt                                                                                                 <  1:553581/145‑1 (MQ=255)                                                    gcGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAgg                                                                                         >  2:672342/1‑151 (MQ=255)                                                    gcGGAAAAACCGCAAGAATCAGCGATTGGATCATGTCGTCCAGCAGATGTGAGGGGCTGATAGCCCCTAAAATCCCAAA‑CGATGTTGGGGCT‑TTGT‑TCGTGGCCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCa                                                                                             <  1:413472/147‑1 (MQ=255)                                                        aaaaaCCGCCAAAATCAGCGATTGAATCATGTCGTCCAGCAGATGGGAGAGGCTGATACCACAAAAAAACCAAAA‑CGTTGTGGGGGCT‑TTGG‑TCGTGGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGaa                                                                                       <  1:711310/147‑1 (MQ=255)                                                          aaaCCGCGAGAATCAGCGATTGGATCATGTCGTCCAGCAGAGGGGAGGGGCTGAAACCACAAAAAATCCAAAA‑CGTTTTTGGGGCT‑TGGT‑TGGTGGACGCCGCCGCACCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTc                                                                                   <  1:868966/151‑1 (MQ=255)                                                           aaTCGCAAAAACCAGCGATTGGATCATGTCGTTCAGCAGTGGGGAGGGGCTAAAAGCACCAAAAATCCAAAA‑CGTTTTGGGGTCT‑TTGG‑GCATTAGCCCGCCCCCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCa                                                                                  <  1:443250/151‑1 (MQ=255)                                                            aTCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCAc                                                                                 >  1:717395/1‑151 (MQ=255)                                                                  aaaaTCAGGGTTGGGATCTTGTCGTGCAGAAGGTGGGGGGGGCTGATACCACCAAAAATACAAAA‑CGTTTTTGGGGTT‑TGGT‑TCTTGGACGCGGCCCCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACttttt                                                                            <  2:870041/148‑1 (MQ=255)                                                                  aGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTc                                                                           >  1:778761/1‑151 (MQ=255)                                                                  aGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTc                                                                           >  2:84383/1‑151 (MQ=255)                                                                  aGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCAAAAATCCCAAA‑CGATGTTGGGGCT‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACtttt                                                                             <  1:742458/149‑1 (MQ=255)                                                                  aGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCAAAAATCCAAAA‑CGATGTTCGGGCT‑TTGG‑TGGTGGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACtttt                                                                             <  1:599661/149‑1 (MQ=255)                                                                        agtgattgattCATGTCGTTCGGCAGGTGGGGGGGGCTGATACCCCAAAAAAACCCAAA‑CGGTGTGGGGGTT‑TTGT‑TCGTTGCCGCGCCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAggg                                                                       <  2:594006/146‑1 (MQ=255)                                                                           gATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCg                                                                  >  2:233243/1‑151 (MQ=255)                                                                             ttGGATCATGTCGTTCAGCAGATGGGAGAGGCTAATAGCCCCAAAAATCCAAAA‑CGTTTTTCGGGCT‑TTGG‑TTTTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCg                                                                  <  2:902531/149‑1 (MQ=255)                                                                             ttGGATCATGTCGTTCAGCAGAGGTGAGAGGCTGATAGCACCAAAAATCCCAAA‑CGATGTTGGGGCT‑TTGG‑TCTTTGACGCAGCCGCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGc                                                                   <  1:897670/148‑1 (MQ=255)                                                                               ggATCATGTCGTTCAGCAGATGTGAGAGGCTGAGAGCCCCTAAAATACCAAA‑CGATGTTGGGGCC‑TTGG‑TCGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGAt                                                                <  1:698880/149‑1 (MQ=255)                                                                               gaataatttCGTTCAGAAGATGGGAGGGGCTGAACGCACCAAAAATACCAAG‑CG‑TGTTTCGGTC‑TTGGTTTTTTGCCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATg                                                               <  2:663647/142‑1 (MQ=255)                                                                                  tCATTTCGTTCAGCAGATGTGAGGGGTGGATACCCCTAAAAACACAAAA‑CGGTTTTCGGCCT‑TTGT‑TCGTTGCCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaa                                                            <  2:394381/150‑1 (MQ=255)                                                                                      gtggtgcagcagTTGGGGGAGGATAAAACCACAAAAAATACAAAC‑CGATGTGGGGGCT‑TGGG‑TGGTTGACGCACCGGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaa                                                            <  1:372059/140‑1 (MQ=255)                                                                                          ttCAGCAGAGGTGAGGGGCTGATAGCACCAAAAATACCAAA‑CGATGTTCGGGCT‑TTGG‑TCGTGGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaaa                                                           <  1:743812/143‑1 (MQ=255)                                                                                          tacagcagGTGTGAGAGGCTTATACCCCCAAAAACCCAAAA‑CTATGTTGGGCCT‑TTTG‑TCTTTGACGCACCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATg                                                      <  2:923785/146‑1 (MQ=255)                                                                                               cagATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAtt                                              >  2:870151/1‑151 (MQ=255)                                                                                                 gATGTGAGTGGCTGATACCCCTAAAAATACCAAA‑GGTTGTGGGGCCC‑TTGG‑TCGCTGACGCCGCGGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                                                 <  1:130220/146‑1 (MQ=255)                                                                                                    gtggggggCTGATAGCCCAAAAAATACAAAA‑CGTTGTTGGGCCT‑TTGG‑TCGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt                                             <  1:725394/141‑1 (MQ=255)                                                                                                    gggagaGGCTGATAGCACCTAAAATACCAAA‑CGATGTGGGGGCC‑TTGG‑TCGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                                                 <  1:1006379/141‑1 (MQ=255)                                                                                                    gggagaGGCTAACAGCACCT‑AATCACCAAATCGATGGGGGGGCT‑TTGA‑TAGTTGACACAGCGACCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTGAtt                                              <  2:203684/144‑1 (MQ=255)                                                                                                      gagaGGCTGATAGCACCAAAAATACCAAA‑CGATGTTGGGGCT‑TTGG‑TCGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa                                         <  1:86593/149‑1 (MQ=255)                                                                                                        gagGCGGAACGCACTAAAAATCCCAAC‑CGTTGTTCGGGCT‑TTGG‑TGGTTGACGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatta                                      <  1:362203/150‑1 (MQ=255)                                                                                                           gCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc                                  >  2:854630/1‑151 (MQ=255)                                                                                                            cTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAcc                                 >  2:644154/1‑151 (MQ=255)                                                                                                               aTAGCACCAAAAATACCAAA‑CGATGTTGGGGCT‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAcc                                 <  1:672360/148‑1 (MQ=255)                                                                                                                 aGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCACCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACAAGACt                            >  2:503909/1‑151 (MQ=255)                                                                                                                      cTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACa                        >  1:391724/1‑150 (MQ=255)                                                                                                                        aaaaTCCCAAA‑CGATGTTGGGGCT‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACAt                       <  1:79227/149‑1 (MQ=255)                                                                                                                           aTACCAAA‑CGATGTTCAGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTg                  >  1:56411/1‑151 (MQ=255)                                                                                                                             aacAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGACACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTg                  >  2:1062596/3‑149 (MQ=255)                                                                                                                              ccAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTAt               >  1:700536/1‑151 (MQ=255)                                                                                                                                aaa‑CGTTTTTCGGGCG‑TTGG‑TCTTTGCCGCACCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGt                 <  1:991719/147‑1 (MQ=255)                                                                                                                                aaa‑CGTTGGTGGGGTT‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGc             <  2:1018861/151‑1 (MQ=255)                                                                                                                                aaa‑CGATGTGGGGGCT‑TTGT‑TCGTTGCCGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGc             <  2:231630/151‑1 (MQ=255)                                                                                                                                 aaagggtgtgggggCT‑TGGG‑TGGTGGCCGCGGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGt                 <  2:201219/136‑1 (MQ=255)                                                                                                                                  a‑CGTTGTTCGGGCC‑TTGG‑TCGTGGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgt           <  1:837190/151‑1 (MQ=255)                                                                                                                                  a‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGcgt           >  1:506150/1‑151 (MQ=255)                                                                                                                                         tgcGGGCT‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCt      <  1:1018679/148‑1 (MQ=255)                                                                                                                                            gggCT‑TTGG‑TCGTGGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTgaag  <  1:318951/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     GGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGATAAATCGCCAGAATCAGCGATTGGATCATGTCGTTCAGCAGATGTGAGAGGCTGATAGCACCTAAAATACCAAA‑CGATGTTCGGGCC‑TTGG‑TCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAG  >  NC_000913/504498‑504786

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