Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	40.9%	122.5 / 16.6	66	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (8/31);  new base G (27/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366377‑366675                                                                                                                                                       |                                                                                                                                                       aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGggt                                                                                                                                                        <  1:11807/150‑1 (MQ=255)  aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGggtgg                                                                                                                                                      <  1:33855/151‑1 (MQ=255)       tGGGGGGCCTAATCAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGCCCAACGCGCGGGGAGAGGC‑GGATTACGGATTGGGCGGCCGggtggactt                                                                                                                                                  >  2:589573/1‑146 (MQ=255)        ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGGGGATTTGGGGGCAGGGGGGgtttt                                                                                                                                                 >  2:337500/1‑150 (MQ=255)          ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGATTTGGGGCCCGGGGGGGttttctt                                                                                                                                              >  1:64089/1‑151 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                                            <  1:136540/150‑1 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                                                                                                            <  2:108053/150‑1 (MQ=255)              ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                                                                                                          <  1:661375/151‑1 (MQ=255)                    gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTg                                                                                                                                    >  2:44868/1‑151 (MQ=255)                     agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                                   <  1:379147/151‑1 (MQ=255)                           cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTTTTGGGCGCCCGGGGGGGTTTTCTTTTCACCAGGgagaacgg                                                                                                                             >  1:370149/1‑147 (MQ=255)                             aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑ggg                                                                                                                            <  2:596550/149‑1 (MQ=255)                             aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCa                                                                                                                          <  2:222616/151‑1 (MQ=255)                              aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCACGGAGAC‑ggg                                                                                                                            >  1:380143/1‑148 (MQ=255)                                tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGGGGGTTTTCTTTTCACCAGTGAGAC‑GGGCaa                                                                                                                         >  1:727533/1‑149 (MQ=255)                                    aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTg                                                                                                                   <  1:560547/151‑1 (MQ=255)                                        aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTg                                                                                                               <  2:223717/151‑1 (MQ=255)                                         ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTCACCAGGGAGAC‑GGGCAACAGCTGAtt                                                                                                                >  1:445653/1‑149 (MQ=255)                                          tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGGGTGTTGGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGTGAC‑GGGGAACAAGTGATTGc                                                                                                              >  1:643991/1‑150 (MQ=255)                                           gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGccc                                                                                                            <  2:734352/151‑1 (MQ=255)                                                 cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCAcc                                                                                                      <  2:418555/151‑1 (MQ=255)                                                  gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCTTGTGTGTTTGGCGCCGGGGGGGGGTTTCTTTTTTCCAGTGAGAC‑GGGGAACAGCTGCTTGGCCTTTAc                                                                                                       >  2:320267/1‑149 (MQ=255)                                                  gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGC‑GGGAGAGGCGGGTTTCGGGTTTGGGGGCCGGGGGGGGTTTTTTTTTC‑CCCGGGAGACGGGGCACAACCTGATTGCCCTTccc                                                                                                       >  2:332767/1‑147 (MQ=255)                                                   cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGAC‑GGGCAACAGCTGATTGGCCTTCCCCGc                                                                                                    >  1:739985/1‑151 (MQ=255)                                                    tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGCGCCAGGGGGGTTTTTTTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACcccc                                                                                                   >  1:778913/1‑151 (MQ=255)                                                     cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                  <  1:95768/151‑1 (MQ=255)                                                       cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTGGGATTTGGCGCCCGGGGGGGTTTTCTTTTTACCCGTGAGAC‑GGGCGAAAGCGGATTTCCCCTTCCCGCCt                                                                                                  >  1:604106/1‑149 (MQ=255)                                                       cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGCGCCAGGGGGGTTTTTTTTTTTACCAGGGAGAC‑GGGCAACAGCTGATTGCCCCTCACCGcc                                                                                                   >  2:230354/1‑148 (MQ=255)                                                         gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGCGGATTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGCGAGAC‑GGGGAACAGCCGATTGCCCtttcccccc                                                                                                   >  2:379472/1‑139 (MQ=255)                                                             gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTga                                                                                          <  1:346172/151‑1 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga                                                                                        <  2:958806/150‑1 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga                                                                                        <  2:330798/150‑1 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGGGGTTTTTCTTTTTACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCCCCCGGccccggg                                                                                         >  1:809335/1‑145 (MQ=255)                                                                   cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGCGCCCGGGGGGTTTTTTTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                                     >  1:428259/1‑150 (MQ=255)                                                                     gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                                   <  2:397537/150‑1 (MQ=255)                                                                     gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGCGTATTGGGCGCCCGGGGGGGTTTTCTTTTTACCCAGGAGAC‑GGGCAACAACCGCTTGCCCTTCACCCCCCGGCCCCGAGAGAGtt                                                                                   >  2:565478/1‑150 (MQ=255)                                                                      tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGGGCCCGGGGGGGTTTTCTTTTCACCAGTGAGAC‑GGGCAACAACTGATTGCCCTTCACCCCCCGGCCCCgaga                                                                                        >  2:974243/1‑144 (MQ=255)                                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgc                                                                                 >  1:700227/1‑149 (MQ=255)                                                                          gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  2:633049/151‑1 (MQ=255)                                                                             aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                          <  1:141449/151‑1 (MQ=255)                                                                                    gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGAC‑GGGCAAAAGCTGATTGCCCCTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCAc                                                                   >  1:662794/1‑151 (MQ=255)                                                                                      gCCGGCTGCATTAATGAAGCGGCCAGCGCCCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCTAGGTTTGTTTTTCTTTTCACCATTGAGAC‑GGGCAACAGCTGAGTGCCCTTCACCGCCTGGCCCTGGGAGAGTTGCAGCGAGCGTTCCACg                                                                  <  2:887977/150‑1 (MQ=255)                                                                                       ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                                <  1:67150/151‑1 (MQ=255)                                                                                        cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGCGCCAGGGGGGTTTTTTTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAGCGGGCCAc                                                                   >  1:547296/1‑147 (MQ=255)                                                                                        cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTg                                                               <  1:900719/151‑1 (MQ=255)                                                                                        cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTg                                                               <  1:776535/151‑1 (MQ=255)                                                                                         aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGCGCCAGGGGGGGTTTTCTTTTTACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGGGAGAGTTGCAGCAAGCGGGCCACGCTgg                                                              >  2:98693/1‑151 (MQ=255)                                                                                             gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAC‑GGGGAACAGCTGATTTTCCCTCCCAGCCTGGCCCCGGGAGAGGTGCAGCAAGCGGGCCCCGCCGGttt                                                           >  1:754468/1‑150 (MQ=255)                                                                                             gCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTGCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCGTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGt                                                             <  2:136353/148‑1 (MQ=255)                                                                                                ttAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAC‑GGGGAAAAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGGTGCAGCAAACGGGCCCCGCTGGTTTGccc                                                       >  1:230010/1‑151 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTGGGCGCCAGGGGGGTTTTTCTTTTTACCCGGGAGAC‑GGGCAACAACTGATTGCCCTTCACCCCCCGGCCCCGGGAGAGGTGCAGCAAACGGGGCACGCCGGGTTTccccc                                                     >  2:1000358/1‑147 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagc                                                   <  1:575421/150‑1 (MQ=255)                                                                                                      aaTCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTGGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGGTTGCCCccgc                                                   >  1:304956/1‑149 (MQ=255)                                                                                                              aaCGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCCGGGTGGTTTTTCTTTTCACCCGTGAGAC‑GGGCAACAGCTGATTGCCCTTCCCCGCCTGGCCCTGAGAGAGGTGCAGCAAGCGGGCCACGCCGGTTTGCCCCCGCAGGCGaaaa                                          >  2:250561/1‑150 (MQ=255)                                                                                                                                ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt                         <  2:695777/150‑1 (MQ=255)                                                                                                                                    tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                    >  2:767861/1‑151 (MQ=255)                                                                                                                                    tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                    >  1:856772/1‑151 (MQ=255)                                                                                                                                    tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                    >  1:477359/1‑151 (MQ=255)                                                                                                                                    tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                     <  2:162868/150‑1 (MQ=255)                                                                                                                                    tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                     <  1:938580/150‑1 (MQ=255)                                                                                                                                     gCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                   >  1:23746/1‑151 (MQ=255)                                                                                                                                        tATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                >  1:127990/1‑151 (MQ=255)                                                                                                                                            gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggggt             >  2:402221/1‑148 (MQ=255)                                                                                                                                            gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat             <  2:210990/150‑1 (MQ=255)                                                                                                                                             ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat           <  2:328863/151‑1 (MQ=255)                                                                                                                                                 ccAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa       <  2:790834/151‑1 (MQ=255)                                                                                                                                                  cAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         >  2:298366/1‑148 (MQ=255)                                                                                                                                                  cAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt      <  2:380137/151‑1 (MQ=255)                                                                                                                                                   aGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg     <  2:445645/151‑1 (MQ=255)                                                                                                                                                     ggtcgtTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa       >  1:1721/1‑147 (MQ=255)                                                                                                                                                       tggtTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAAATGAGc  >  2:842896/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366377‑366675

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