Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,785	0	A	G	43.1%	97.4 / 20.0	65	F25S (TTC→TCC)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (3/34);  new base G (28/0);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.29e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCC  >  NC_000913/1535637‑1535927                                                                                                                                                     |                                                                                                                                                    gTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGaaa                                                                                                                                                   <  1:1240492/151‑1 (MQ=255)   gACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGgaa                                                                                                                                                   >  1:1246245/1‑149 (MQ=255)    aCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaaca                                                                                                                                                <  2:102127/151‑1 (MQ=255)     ccGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaacaa                                                                                                                                               <  2:43104/151‑1 (MQ=255)         aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaaaaac                                                                                                                                           >  2:194129/1‑147 (MQ=255)         aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaa                                                                                                                                               >  2:1072448/1‑147 (MQ=255)         aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAaaac                                                                                                                                            >  2:806546/1‑150 (MQ=255)           aTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCgg                                                                                                                                         <  2:1177822/151‑1 (MQ=255)             cGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑a                                                                                                                                      <  1:900723/151‑1 (MQ=255)                 aCTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGaaaaaaaa                                                                                                                                             >  1:482579/1‑139 (MQ=255)                  cTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACACCGGC‑AG‑aa                                                                                                                                  >  2:512828/1‑149 (MQ=255)                   ttCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGGC‑Aggaat                                                                                                                                 >  2:1175286/1‑146 (MQ=255)                    tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCAGC‑CG‑aaat                                                                                                                                >  1:125669/1‑147 (MQ=255)                    tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AAtt                                                                                                                                <  2:778259/149‑1 (MQ=255)                    tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AAtt                                                                                                                                <  1:1157851/149‑1 (MQ=255)                       aCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCGCC‑CG‑AATTGGcg                                                                                                                            >  1:1206812/1‑150 (MQ=255)                          gTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGaa                                                                                                                        <  1:518669/151‑1 (MQ=255)                             ggCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAc                                                                                                                      <  2:247623/150‑1 (MQ=255)                             ggCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAc                                                                                                                      <  2:35021/150‑1 (MQ=255)                              gCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑gg                                                                                                                   <  1:644880/151‑1 (MQ=255)                               cGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGC‑CG‑AAATTGCGCGGAAccg                                                                                                                    >  1:697217/1‑148 (MQ=255)                                gTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGGC‑AG‑AATTTGCGCGgaa                                                                                                                       >  2:853812/1‑146 (MQ=255)                                 tAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCGGG‑AATGGCCCGGAAAC‑gg                                                                                                                   >  1:535915/1‑149 (MQ=255)                                  aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGC‑AA‑GATTGGCGCGGAAC‑gggt                                                                                                                 >  2:1117879/1‑147 (MQ=255)                                    aGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAACAAAACAACGGC‑GGCAAAAGCCGCGaa                                                                                                                        >  1:229894/1‑142 (MQ=255)                                          ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTaaa                                                                                                       <  1:474212/151‑1 (MQ=255)                                              gcgcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGc                                                                                                    <  2:314314/150‑1 (MQ=255)                                                      cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAACCGGC‑AG‑AAATTTCGCGAAAC‑GGGTGTGAGGGAAAAGAATCg                                                                                                >  2:95421/1‑146 (MQ=255)                                                       gCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACAACCGC‑AG‑AAATGCCGCGCAAC‑GGGTGTGGGGGAAAAGAAACGACgggg                                                                                          >  2:1249820/1‑151 (MQ=255)                                                          gtcgtcGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGC‑AG‑AAATGCCGCGCAAA‑GGGTGTGGGGGAAAAACAACGACggggg                                                                                         >  2:183674/1‑148 (MQ=255)                                                                     gcagcaGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGttt                                                                           <  1:1033208/151‑1 (MQ=255)                                                                       agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCa                                                                         <  1:995469/151‑1 (MQ=255)                                                                          agAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCGGC‑AG‑GAATGCCGCGAAAA‑GGGTGGGAGGGGAAAAAATCGACGGGGGGCA‑cccccgggtttttt                                                                       >  1:153297/1‑138 (MQ=255)                                                                          agAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTa                                                                     <  2:440998/151‑1 (MQ=255)                                                                                     ttCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCt                                                           <  1:313662/150‑1 (MQ=255)                                                                                       ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGGAAAAACAACAGC‑GG‑AATTTTCGCGGAAA‑GGGTGTGTGGGAAAAGCAACATCCCGGGGCA‑CCTGCCGTTTCC‑TTTACggttttccgt                                                         >  2:648214/1‑141 (MQ=255)                                                                                         gAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGAc                                                      <  2:575190/151‑1 (MQ=255)                                                                                              gCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGC‑CG‑GAATTCCCCGGAAC‑GGGTGTGAGGGAAAAAAATCGACCTGGGACA‑CCCGCGGGTTCT‑TTTAAGTTTTCCTTTGTCTgggg                                                 >  2:163903/1‑149 (MQ=255)                                                                                               ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGC‑AG‑AAATGCCGCCGAAA‑GGGTGGGGGGGAAAAAAAACCACGCGGGACA‑CCCGCGGGTTTA‑TTAAAAGTCTCCtcttgctgggg                                                 >  2:57868/1‑140 (MQ=255)                                                                                                cTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGAC‑AG‑AATAGGCGCGCAAC‑AGTGTTGAGGTGAAAAAAAAGACGTGGGGCA‑ACCGCGGCTTCATTTAACGtattctct                                                          >  1:90678/1‑134 (MQ=255)                                                                                                  aaCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGaa                                               <  1:1056644/149‑1 (MQ=255)                                                                                                       gcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGC‑AG‑AAATGGCGCGGAAA‑GGGTGTGGGGGAAAAAGAACGACGGGGGACA‑CCTCCGGGTTTC‑TTAACGTTTTTCtgttgttgggggaatat                                           >  1:584864/1‑134 (MQ=255)                                                                                                         gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGC‑CG‑AATTTTCGCGGAAA‑GGGTGTGGGGGGAAAAAAACCACGTGGGACA‑CCTGCGGGTTTA‑TTTAAGTTTTCCTCTTACTgggg                                                 >  1:900916/1‑138 (MQ=255)                                                                                                         gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAg                                      <  1:355027/151‑1 (MQ=255)                                                                                                             gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTg                                  <  2:1049866/151‑1 (MQ=255)                                                                                                               tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGG‑AG‑AATTAGCGCGGAACAGGTTTTAGGGGAAAAG‑AACCAGGTGGGAAA‑ACCCCGGTTTTA‑TTTAACTTTTCCCTTGACTAGGGAAATATTATAGCCtccc                                >  2:1019137/1‑148 (MQ=255)                                                                                                                aaGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCa                                <  2:915682/150‑1 (MQ=255)                                                                                                                 aGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAg                              <  1:331247/151‑1 (MQ=255)                                                                                                                 aGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCCGC‑CG‑GAATTCCGCGGAAC‑GGGTGTGAGGGAAAAACAACGACGTGGGGCA‑ACCTCGGGTTTA‑TTTAAGTTCTCCTGTGACTGGAGGACTATTATAgac                                    >  2:354556/1‑143 (MQ=255)                                                                                                                 aGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGCGGAAAAAAAACCGGC‑AG‑AATTCCCGAGGAAA‑GGCGGTGGGGGGAAAACATATGCCCGGTAAACCCTCCGGGTTC‑‑TTCACGTTCTCCCGTGACTGGGGGACAATCAtca                                      >  1:1071178/1‑141 (MQ=255)                                                                                                                  gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAACACCCGC‑CG‑AAATGCCGCGAAAA‑GGGTGTGGGGGAAAAAAATCGACGTGGTACA‑ACTGCGGGTTCT‑TTTACGTTTTCCTCTGACTGGAGAACTATCATAGcc                                    >  2:81230/1‑144 (MQ=255)                                                                                                                  gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGt                             <  2:93562/151‑1 (MQ=255)                                                                                                                    tGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCGGC‑CG‑AAATGCCGCGAAAC‑GGTTGTGGGGTAAAAAAAACGACGGGGGACA‑CCTGCGGGTTTA‑TTTAAATTCTCCTTTGACTGGGGAAATATCATAACCTGCAAGTgg                           >  2:1080162/1‑151 (MQ=255)                                                                                                                     gCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGc                          <  2:1102441/151‑1 (MQ=255)                                                                                                                      cATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAACCCGC‑CG‑AATTGCCGCGAAAC‑GGGTGTGAGGGAAAAACAACGGCGTGGTACA‑CCTGCGGGTTCT‑TTTACGTTCTCCTGTGGCTGGAGAAATAACATAACCTGCAAGTgggc                         >  2:1094832/1‑149 (MQ=255)                                                                                                                      cATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGcc                         <  2:512357/151‑1 (MQ=255)                                                                                                                           gcgcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgag                     <  2:319258/150‑1 (MQ=255)                                                                                                                            cgcTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGC‑CG‑GAATGCCGCGGAAA‑AGTTGTGTGGGAAAAACAAAGACGTGGGACA‑ACCGCGGGGTCA‑TTAAAGTTCTCCTGTGGCTGGAGAAATATAATAACCTGCaaaggggcgggg                     >  1:670147/1‑139 (MQ=255)                                                                                                                             gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGc                  <  1:972242/151‑1 (MQ=255)                                                                                                                                tCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGa                <  1:755044/150‑1 (MQ=255)                                                                                                                                 cGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa               <  1:163899/150‑1 (MQ=255)                                                                                                                                 cGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAg              <  1:654144/151‑1 (MQ=255)                                                                                                                                  gCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg             <  2:1068830/151‑1 (MQ=255)                                                                                                                                     aaGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGc           <  1:336844/150‑1 (MQ=255)                                                                                                                                     aaGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCt          <  1:1072386/151‑1 (MQ=255)                                                                                                                                          tATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATcc      <  1:1221575/150‑1 (MQ=255)                                                                                                                                               gcggGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGcc  >  1:1183722/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGC‑AG‑AATTGCCGCGAAAC‑GGTTGTGAGGTAAAAGCATCGACGTGGTACA‑CCTGCGGTTTCA‑TTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCC  >  NC_000913/1535637‑1535927

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