Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,504,799	0	A	G	43.3%	88.2 / 23.1	60	G361G (GGA→GGG)	insQ	IS609 transposase B
Reads supporting (aligned to +/- strand):  ref base A (4/30);  new base G (26/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAC  >  NC_000913/1504654‑1504948                                                                                                                                                   |                                                                                                                                                         gTTGCCCGTACTGTGGTCATACCGGAAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGAt                                                                                                                                                     <  1:288510/151‑1 (MQ=255)    gcgcGTACTGTGGTCATCCAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGt                                                                                                                                                  <  1:58059/151‑1 (MQ=255)         tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGTTCTTTTCCCTGGGGAGAAGGTgg                                                                                                                                                >  2:837176/1‑147 (MQ=255)               ggTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGGTCTTGCCCGTGGGGAGAAGGGGGAGGCAg                                                                                                                                          >  1:208892/1‑148 (MQ=255)                     aCAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGt                                                                                                                                <  1:788799/151‑1 (MQ=255)                     aCAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCCTTTTTTCCTGTGGGGGGATGGTGCGGTCCAGCCG‑ccc                                                                                                                                  >  2:730647/1‑149 (MQ=255)                        gcaaaaGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGa                                                                                                                             <  2:420256/148‑1 (MQ=255)                        gcaaaaGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGgccg‑c                                                                                                                                    <  1:290436/141‑1 (MQ=255)                         cGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGGTCTTTCCTGTGGGGGGATGGGGCAGTCAGGCCG‑ccc                                                                                                                                  >  2:479927/1‑145 (MQ=255)                          gAAGGAGAACCGCCTTTCTCAAAGTAATTTCAGATGCCTGGTATGTGGATATACAGCGAACGCCGATGTAAGTGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGa                                                                                                                             <  2:182798/149‑1 (MQ=255)                          gAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCCCTCTTTCCTGTGGGGAGATGGGGCAGTCAGGCCG‑ccc                                                                                                                                  >  1:481/1‑144 (MQ=255)                            aaGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGGTCTTGCCTGTGGGGAGAAGGGGGAGTCAGGCCG‑ccccgt                                                                                                                               >  2:650222/1‑142 (MQ=255)                             agagAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGGTCTTTCCTGTGGGGAGAggggggaggcaaggc                                                                                                                                       >  2:678232/1‑125 (MQ=255)                                gaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACACCCGTCCTCCCTGTGGGGGGAGGGGGAAGTCCGGCCGCCCCGTTGaaaacg                                                                                                                        >  1:550635/1‑145 (MQ=255)                                   ccGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑CCGCCGGTCTTTCCTGTGGGGGGATGGGGCAGGCAGGCCG‑CCCCTTTGAGCAGG‑AAacc                                                                                                                 >  2:543659/1‑151 (MQ=255)                                      ccTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGGTCTTTCCTGTGGGGAGAAGGGGGAGGCAGGCCG‑CCCCTTGAAACAGG‑AAACcc‑cc                                                                                                             >  1:1052792/1‑151 (MQ=255)                                      ccTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑cc                                                                                                             <  2:1006914/151‑1 (MQ=255)                                       cTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGGTCTTGCCTGTGGGGAGATGGGGCAGGCAGGGCG‑CCCGTTGGAGCAGG‑AAcccc‑cc                                                                                                             >  2:908747/1‑150 (MQ=255)                                        tGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGCTCTTTCCTGTGGGGGGAAGGTGGAGGCAGGCCC‑CCCCCTGAAACAGG‑AAacca‑cc                                                                                                             >  2:338716/1‑149 (MQ=255)                                             caAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGa                                                                                                      <  2:747802/151‑1 (MQ=255)                                              aaaGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑AAGCCCGTCTTGCCTGTGGGGAGAGGGGGTAGTCCGGCCC‑CCCCCTGAACGATG‑AACCCA‑CCGAAATGAt                                                                                                     >  2:198361/1‑151 (MQ=255)                                               aaGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGAtt                                                                                                    <  2:952440/151‑1 (MQ=255)                                                      ttCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTTCCTGTGGGGAGATGGGGGAGGCAGGCCG‑CCCCTTGAAGCAGG‑AACCCA‑CCCAAAAGATTTAGGCg                                                                                              >  2:1032327/1‑150 (MQ=255)                                                            tGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGCTCTTGCCTGGGGGGGGGAGGGGGAGTCAGGGCG‑CCCGGTGTAGGAGGAAACCCC‑CCGAAA‑GAATTAAGCGACAggc                                                                                        >  1:846066/1‑148 (MQ=255)                                                              ccAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGa                                                                                     <  2:236861/151‑1 (MQ=255)                                                               cAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGaa                                                                                    <  1:1175833/151‑1 (MQ=255)                                                                aGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCAACGCCGTT‑TTTCCCGTGGGGGAAAGGGGGAATGAAGGCG‑CCCGATGTAGCAGG‑AAACCA‑CCCAAAAGAATTTGGCGGCAGGCTgta                                                                                    >  2:351394/1‑148 (MQ=255)                                                                 ggTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACg                                                                                  <  2:555410/151‑1 (MQ=255)                                                                      gtgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTTCCTGTGGGGAGATGGTGGAGGCAGGCCG‑CCCCGTGAAGCAGG‑AACCCA‑CCGAAAAGATTTAGGCGACAGCCTGAACGTAg                                                                               >  2:762528/1‑149 (MQ=255)                                                                      gtgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGGGGGGGGGAGGGGCAGTCAGGGCT‑CCCCTTGAAGCAGG‑AAACCA‑CCGAAAAGATTTAGGGGGCAGCCCGaa                                                                                    >  2:113399/1‑144 (MQ=255)                                                                             taCAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGGTCTTTTCTGTGGGGAGAAGGGGGAGTCAAGGCG‑CCCCCTGGAACAGG‑AAAACACCCGAAAGAATTAAGG‑GACAGCCCGaaattaaacagg                                                                          >  1:255790/1‑137 (MQ=255)                                                                               cAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCCTTCTTTCCTGTGGGGAGATGGGGGAGTCAGGCCG‑CCCCTTGAAACAAG‑AAACCA‑CCGAAAAGAATCAGGCGACAGCCTGAACGGAGCAGGGAGcc                                                                      >  2:213160/1‑149 (MQ=255)                                                                               cAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACCCCCTTCTTTCCTGTGGGGGGATGGGGTTGGCAGGCCG‑CCCCCTGAAATAAG‑AACCCC‑CCCCAATGATTAAGGCGACAGCCCGAAAGTAGCAagg                                                                          >  2:706659/1‑145 (MQ=255)                                                                                      cGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCtt                                                              <  2:947210/150‑1 (MQ=255)                                                                                       gCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCa                                                            <  2:437332/151‑1 (MQ=255)                                                                                        ccGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCCCTCTTGCCTGTGGGGGGGTGGTGGAGTCAGGGCG‑CCCGTTGTAGCAGG‑AAACCA‑CCCCAAAGCTTCAGGCGGCAGCCCGCACGGAGCAGGGGTCCACGccc                                                                >  1:83707/1‑146 (MQ=255)                                                                                          gATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGG‑ACGCCGCTTTTTCCTGTGGGGAGAGGGGGGAGGTCGGGCG‑CCCGTTGCAGCAGG‑AAACCC‑CCGAAAAAATTTAAGCGGCCACCCGAACGTGGCAGGGGGTCCCGTCCTTc                                                             >  1:561394/1‑147 (MQ=255)                                                                                          gATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAggg                                                         <  2:298426/151‑1 (MQ=255)                                                                                           aTGTAAATGGCGCTCGCAACATTTTAGCGGTGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGc                                                        <  2:473199/151‑1 (MQ=255)                                                                                           aTGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGc                                                        <  2:729049/151‑1 (MQ=255)                                                                                               aaaTGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg                                                     <  1:141850/150‑1 (MQ=255)                                                                                                 aTGGCGCTCGCAACATTTTAGCGGCGGGGG‑GCGCCCCTCTTTCCTGTGGGGAGAGGGGGCAGGCAGGCCG‑CCCGCTGTAGCAGG‑GAACCA‑CCCAAAAGATTCAAGCGGCAGACTGGACAGAGCAGGGGTCCACGGCCTTCAcggggggggg                                                  >  1:621899/1‑144 (MQ=255)                                                                                                     cgcTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATg                                              <  1:420785/151‑1 (MQ=255)                                                                                                      gcTCGCAACATTTTAGCGGCGGGGG‑GCGCCGTTTTTGCCTGTGGGGGGGAGGGGGAGTCTGGGCG‑CCCCTTGGAGCCGG‑AAACCA‑CCCAAAAGATTCAGGCGACAGCCTGGAAGGAGCAGGGAGCCACGCCCCTCAGGGCGTgggggggg                                              >  1:874124/1‑147 (MQ=255)                                                                                                            aaCATTTTAGCGGCGGGGC‑ACGCCGGTCTTTCCTGTGGGGAGATGGTGTAGGCAGGCCG‑CCCGCTGTAGCAGG‑AAACCA‑CCCCAAAAATTTAGGCGGCAGCCCGAAAGGAGCAGGGGTTCACCTCCCTCTGGGCGGGGAGGAGGTCAata                                        >  1:434535/1‑150 (MQ=255)                                                                                                             aCATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAatat                                       <  1:508327/150‑1 (MQ=255)                                                                                                              cATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCg                                     <  2:130202/151‑1 (MQ=255)                                                                                                                   tAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGttt                                 <  2:1144484/150‑1 (MQ=255)                                                                                                                     gcggcgGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTg                              <  1:479916/151‑1 (MQ=255)                                                                                                                           gggC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCg                        <  2:1129200/151‑1 (MQ=255)                                                                                                                           gggC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTAAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCg                        <  1:221943/151‑1 (MQ=255)                                                                                                                                    cGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAAtt                >  1:1079257/1‑151 (MQ=255)                                                                                                                                       tCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCa              >  1:257350/1‑150 (MQ=255)                                                                                                                                          tGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCaa             <  2:909423/148‑1 (MQ=255)                                                                                                                                          tGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAaca           <  2:1046997/150‑1 (MQ=255)                                                                                                                                           gCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcacc         <  2:1108386/151‑1 (MQ=255)                                                                                                                                           gCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGCGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCaa             >  1:348038/1‑147 (MQ=255)                                                                                                                                             cTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcaccac       >  1:796643/1‑151 (MQ=255)                                                                                                                                               gtgGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGt     <  1:965631/151‑1 (MQ=255)                                                                                                                                               gtgGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGt     <  1:587837/151‑1 (MQ=255)                                                                                                                                                   agagaTGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAc  >  1:1065418/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                         GTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGC‑ACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCG‑CCCGTTGAAGCAGG‑AACCCA‑CCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAC  >  NC_000913/1504654‑1504948

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