Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,791	0	A	C	41.0%	176.0 / 16.6	95	V161G (GTG→GGG)	basS	sensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (50/6);  new base C (0/39);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 1.14e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAT  >  NC_000913/4332642‑4332936                                                                                                                                                         |                                                                                                                                                  cATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CAc                                                                                                                                                  >  2:879332/1‑151 (MQ=255)    gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CAccc                                                                                                                                                <  1:732173/150‑1 (MQ=255)       taTAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGc                                                                                                                                            >  1:90285/2‑151 (MQ=255)       gATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGc                                                                                                                                            >  2:470147/1‑151 (MQ=255)        aTAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGcc                                                                                                                                           >  1:124880/1‑151 (MQ=255)        aTAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGcc                                                                                                                                           >  2:415115/1‑151 (MQ=255)         tAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGcc                                                                                                                                           >  2:443753/1‑150 (MQ=255)          aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGcc                                                                                                                                           >  2:867679/1‑149 (MQ=255)           aTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGt                                                                                                                                        >  2:1125982/1‑151 (MQ=255)            ttACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTg                                                                                                                                       >  1:438933/1‑151 (MQ=255)             tACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTgg                                                                                                                                      >  2:880074/1‑151 (MQ=255)                 ggAAAAAAATGCCTGTCCGGCCCGCGCCAGTGGCAGCAGCGGGGAGACGCTCCCCATCATCTGATCAAGCCGTGCCACAAACGGTGCTACATCAATGGGATGCGTTTTCCCCAGCAGTCCA‑AGATGCAAC‑CC‑CCCCCCCGCCAGTGGCg                                                                                                                                    <  1:783657/149‑1 (MQ=255)                  gaagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTc                                                                                                                                 >  1:1399785/1‑151 (MQ=255)                  gaagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTc                                                                                                                                 >  2:1140968/1‑151 (MQ=255)                    agaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGc                                                                                                                               >  2:1290500/1‑151 (MQ=255)                      aaaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGc                                                                                                                               >  2:300807/1‑149 (MQ=255)                       aaaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGt                                                                                                                            >  2:206555/1‑151 (MQ=255)                        aattccTGTCCGGCCCGCGCCGGTGGCAGCGGCGGGGGGCCCCTCCCATTCATTTGATCAAGCCGGCCCACTAACGGTGTTACATCAATGTGATGCTTTTTCGCCGGCAGTCCC‑AGAGGCAAA‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                           <  2:749030/146‑1 (MQ=255)                          tGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                         >  2:1391178/1‑151 (MQ=255)                          tGACTGTCCGGCACGCGCAAGTTGCAGCAGCTGGGCGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                         >  1:1429/1‑151 (MQ=255)                          tGACTGTCCGCCCCGCGCCAGTTGCAGCAGCGGGGAGCCCCTCTCCATCATCTGATCAAGCCGTCCCACTACCGGTGCTACATAAATGTGTTGCTTTTTCGCCAGCAGTCCC‑AGATGAAAA‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                          <  2:1448598/150‑1 (MQ=255)                            aCTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTg                                                                                                                       >  1:1214313/1‑151 (MQ=255)                                 ccGGCCCGCGCCAGTTGCAGCGGCGGGGAGCCGCTCCCCATCATCTAACCAAGCGGTCCCACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCAGTCCC‑AGATGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                         <  2:325962/144‑1 (MQ=255)                                 ccGGCCCGCGCAAGTTGCAGCGGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAGCGGTGCTACATCAATGTGATGCTTTTCCGCCAGAAGTTCC‑AGATGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                                      <  2:1501160/147‑1 (MQ=255)                                 ccGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGcgcg                                                                                                                   >  2:1192196/1‑150 (MQ=255)                                  cGGCACGCGCCATTTGCGGCGGCGGGGAGCCGCTCTCCACCATCTGATAAAGCCGTCCCACTAACGGTGCTACATCAATGTGATGCTTTTCCGCCGGCAGTTCA‑AGATGAAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGAc                                                                                                                 <  2:1314242/151‑1 (MQ=255)                                   gccccgcgCCAGTGGCAGCGGGTGGGAGACCCTCCCCATCTTCTGATCAAGCCGTGCCACTATGGGTGCTACATCAATGTGTTGCGTTTTCGCCAGCTGTTCC‑AGAGGAAAC‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTgcgc                                                                                                                    <  2:1088207/143‑1 (MQ=255)                                     cACGCCCCAGTGGCGGCGGGGGGGAGCCGCTCCCCACCATCTGATAAGGCGGTCCCACAAACGGTGCTACACAAATGTGATGTTTTTTCCCCGGAAGTTCA‑AGATGCAAC‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGAc                                                                                                                 <  2:733615/148‑1 (MQ=255)                                          gcaattggcagcaGTGGGGAGCCCCTCCCAATCTTCTGACAAAGCCGTCCCACAAAGGGTCCTACACAAATGTGTGGTTTTTCCGCCAGCAGTCCA‑AGTTGAAAC‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCa                                                                                                             <  2:1501406/140‑1 (MQ=255)                                               tggaggcagcGGGGAGCCGCTCCCCATCTTTTGATCAAGCCGTGCCACTAGCGTTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                                       <  2:1303784/143‑1 (MQ=255)                                                    gcaggggggggcccctctCCTTCATCTGATAAAGCCGTCCCACAAACGTTGCTACTTCAATGTGTGGCTTTTCCGCCAGAAGTCCG‑AGA‑GCAAACCG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                    <  2:1051721/132‑1 (MQ=255)                                                          gggAGCCCCTCTCCATCTTCTGATAAAGCCGTGCCACAACCGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTTCC‑AGAGGAAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGtt                                                                                         <  1:185444/151‑1 (MQ=255)                                                          gggAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGtt                                                                                         >  2:58216/1‑151 (MQ=255)                                                               aCGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATcc                                                                                     >  1:1176265/1‑150 (MQ=255)                                                                cgttctcCACCTTTTGATCAAGCCGTCCCACTAGCGGTGCAACATCAATGTGTTGCTTTTCCGCCAGCAGTTCC‑AGATGAAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                   <  1:688903/148‑1 (MQ=255)                                                                cgccctccactttcTGAACAGCCCGTGCAACTAACGGTCCTACAACAATGGGAGGCTTTTTCGCCAGCATTCCC‑AGAGGCAAC‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCa                                                                                    <  1:1392936/138‑1 (MQ=255)                                                                ccctctCCATCCTCTGATCAAGCCGTGCCACAAACGGTGCAAAATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGAAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                                       <  1:1385106/145‑1 (MQ=255)                                                                ccctccCCATCATCTGATCAAGCCGCCCCACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCGTTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                   <  1:1431917/149‑1 (MQ=255)                                                                ccctccCCATCATCTAACCAAGCCGTGCCCCTAAGGGTGCTCCACCAAGGTGATGCTTTTCCCCCAGCATTCCC‑AGATGCAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAg                                                                                   <  1:139425/149‑1 (MQ=255)                                                                 gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                  >  2:1367226/1‑151 (MQ=255)                                                                 gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                  >  2:1167945/1‑151 (MQ=255)                                                                   tctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGt                                                                                >  1:882278/1‑151 (MQ=255)                                                                    ccccATCCTCGGAACAAGCCGCCCCACGAACGGGGCCACATAAATGTGTTGCTTTTTCCCCAGCATTTCA‑AGACGAAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCa                                                                                    <  1:1326289/144‑1 (MQ=255)                                                                        atcatcTAATCAGGCCGCGCCACAAGGGGGGCTACATCAATGTGATGCTTTTCCGCCAGCAGTCCA‑AGAGACAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                            <  1:1290264/150‑1 (MQ=255)                                                                         tcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGt                                                                          >  2:1353892/1‑151 (MQ=255)                                                                           atcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                         >  2:537749/1‑150 (MQ=255)                                                                           atcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                         >  1:164590/1‑150 (MQ=255)                                                                              taatcagcccGTCCCACTACCGGTGCTACATCAATGTGTTGTTTTTTCCCCAGCAGTCCC‑AGAGGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCa                                                                        <  1:272507/140‑1 (MQ=255)                                                                              tGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                     >  2:436811/1‑151 (MQ=255)                                                                               gatcagccctgcccACTACCGGTGCTACATAAATGTGATGGTTTTTCCCGAGC‑GTTCCGAGAGGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                            <  1:928241/131‑1 (MQ=255)                                                                               gATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                    >  2:493659/1‑151 (MQ=255)                                                                                   accccGTCCCACAACCGGGGCTAAATCATTGTGAGGTGTTTCCCCCAGCAGTCCC‑AGAGGCAAA‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAg                                                                       <  1:749532/141‑1 (MQ=255)                                                                                       cGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                             >  1:875082/1‑150 (MQ=255)                                                                                        ggcccACAACCGGTGTTACATCAATGTGTGGGGTTTCCGCCAGCTGTTCC‑AGATGCAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAg                                                           <  1:1375625/148‑1 (MQ=255)                                                                                        gccccACTAACGGGGCTACATAAATGTGATGCTTTTTCGCCAGCATTTCC‑AGAGGAAAG‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCa                                                            <  2:1170105/147‑1 (MQ=255)                                                                                           ccACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                             >  2:1207748/1‑146 (MQ=255)                                                                                           ccACTAACGGTGCTACATAAATGTGATGCTTTTCCGCCAGCATTTCC‑AGATGAAAA‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                        <  2:229243/151‑1 (MQ=255)                                                                                           ccACTAACGGGGTTACAAACATGCGATGTTTTTACGCAAGC‑GTCGC‑AGGAGGAAA‑CGCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCata                                                          <  1:1161419/149‑3 (MQ=255)                                                                                               tAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                    >  1:886189/1‑151 (MQ=255)                                                                                               tAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                    >  1:1282279/1‑151 (MQ=255)                                                                                                aagggTGCTACATAATTGTGATGCGTTTTCGCCAGCAGTCCC‑AGATCCAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGAt                                                         <  1:1328255/142‑1 (MQ=255)                                                                                                aaCGGTCATACACAAATGGGTTGCGTTTTCCCCAGCAGTCCA‑ATATGAAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGtt                                                     <  2:385997/149‑1 (MQ=255)                                                                                                aaCGGGGCTACATAAATGTGATGCGTTTTCGCCAGCGGTTCC‑AGAGGCAAC‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTa                                                   <  2:981014/151‑1 (MQ=255)                                                                                                   ggTGCTACATCATTGTGATGCTTTTCCGCCAGCATTTCA‑AGATGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTaa                                                  <  2:1384427/149‑1 (MQ=255)                                                                                                     tGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGCCTGACCAGATCGTTTAACg                                                >  2:1465784/1‑149 (MQ=255)                                                                                                         aCATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaa                                           >  2:331211/1‑150 (MQ=255)                                                                                                           aTCATTGTGAGGTGTTTTCGCGAGCTGCCCA‑AGCGGCAAG‑CC‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCATCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                               <  1:914302/144‑1 (MQ=255)                                                                                                           aTCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGTCGACTGACCAGATCGTTTAACGCCGAAacc                                        >  1:522427/1‑151 (MQ=255)                                                                                                              aaTGTGATGCTTTTTCCCCAGCAGTCCC‑AGATGAAAC‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                     <  2:1227687/151‑1 (MQ=255)                                                                                                                tgtgATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                                   >  1:609890/1‑151 (MQ=255)                                                                                                                   gttGCTTTTCCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTc                                <  2:592892/149‑1 (MQ=255)                                                                                                                    aTGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGGCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                               >  2:774951/1‑151 (MQ=255)                                                                                                                     tGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCg                               >  2:827950/1‑150 (MQ=255)                                                                                                                      gCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACACCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAt                             >  2:15824/1‑151 (MQ=255)                                                                                                                      gCGTTTCCGCCAGAAGTTCC‑AGATGAAAC‑CC‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacca                                       <  2:226903/141‑1 (MQ=255)                                                                                                                          ttccgcgagcgtttCC‑AGTTGAAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                          <  2:1430326/138‑1 (MQ=255)                                                                                                                           ttCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                         >  1:909042/1‑150 (MQ=255)                                                                                                                              gCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                     >  2:1361064/1‑151 (MQ=255)                                                                                                                                  gcagTCCC‑AGAGACAAA‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgc                 <  1:1142376/151‑1 (MQ=255)                                                                                                                                   cagTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcg                >  2:907139/1‑151 (MQ=255)                                                                                                                                     gTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCt              >  1:878525/1‑151 (MQ=255)                                                                                                                                      tccc‑AGATGCAAA‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                     <  2:1582208/141‑1 (MQ=255)                                                                                                                                       tCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt            >  1:1319485/1‑151 (MQ=255)                                                                                                                                       tCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt            >  2:1310907/1‑151 (MQ=255)                                                                                                                                         caagataaaaa‑CG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg           <  2:134278/143‑1 (MQ=255)                                                                                                                                         c‑AGATGCAAC‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt            <  1:536445/149‑1 (MQ=255)                                                                                                                                          aatgagaaaa‑CC‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTgg                  <  2:982677/136‑1 (MQ=255)                                                                                                                                           aGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGa          >  1:970218/1‑150 (MQ=255)                                                                                                                                            gataaaaa‑CC‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg             <  1:1101007/141‑1 (MQ=255)                                                                                                                                                aaaa‑CG‑CACCCCCGACAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGGGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCCACTGACCCGATCGTTTAACGCCGAAACCACCGCTTCGCTTTCGAGGGCGGCGCTGTGAATGGc     >  1:461235/2‑150 (MQ=255)                                                                                                                                                  aaccG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCaa   <  1:415085/148‑1 (MQ=255)                                                                                                                                                  aa‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGCGGGTGGCGCTGTGAATGGCa    >  1:1122115/1‑149 (MQ=255)                                                                                                                                                  aa‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt  >  1:1500187/1‑151 (MQ=255)                                                                                                                                                   accG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt  <  2:1232304/149‑1 (MQ=255)                                                                                                                                                     cG‑CCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt  <  1:781364/149‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                  CATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑AGATGCAAA‑CG‑CACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAT  >  NC_000913/4332642‑4332936

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