Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,326,450	0	T	G	52.4%	81.0 / 103.7	63	G114G (GGT→GGG)	atoB	acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/28);  new base G (33/0);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 6.91e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-12

CGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAA  >  NC_000913/2326302‑2326599                                                                                                                                                     |                                                                                                                                                       cGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTAt                                                                                                                                                      <  2:744946/151‑1 (MQ=255)  ggCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATg                                                                                                                                                     <  1:602191/151‑1 (MQ=255)   gcgcGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATgg                                                                                                                                                    <  2:755934/151‑1 (MQ=255)   gcgcGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                                    >  2:972553/1‑151 (MQ=255)       gTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGaaa                                                                                                                                                 >  2:721383/1‑150 (MQ=255)        tCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaa                                                                                                                                                >  1:651839/1‑150 (MQ=255)          aGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaa                                                                                                                                                >  2:83352/1‑148 (MQ=255)            gCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGgaaaaa                                                                                                                                              >  2:776788/1‑146 (MQ=255)               cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGaa                                                                                                                                                  <  1:479077/141‑1 (MQ=255)               cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                                    >  2:479053/1‑139 (MQ=255)               cTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGtt                                                                                                                                        >  1:343531/1‑151 (MQ=255)                        aGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTTTGGGAAATATGAGTTTTGtcccc                                                                                                                                >  1:196490/1‑150 (MQ=255)                        aGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGccccc                                                                                                                                >  1:1017103/1‑150 (MQ=255)                         gCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCta                                                                                                                              <  1:366796/151‑1 (MQ=255)                         gCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATTGGAAATATGAGATTTGccccc                                                                                                                                >  2:617778/1‑149 (MQ=255)                           gggCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCtact                                                                                                                            <  2:369832/151‑1 (MQ=255)                                 gCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCg                                                                                                                      <  1:971694/151‑1 (MQ=255)                                   aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATAGGTTTAGCCCCCcacttact                                                                                                                        >  2:16372/1‑147 (MQ=255)                                      aaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCCACTTAATCGGTGc                                                                                                                  >  2:742243/1‑150 (MQ=255)                                          gtgtGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAg                                                                                                             <  2:6733/151‑1 (MQ=255)                                           tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGTAAATATGAGGTTAGACCCCCACGTTCTCgccgccaaa                                                                                                              >  2:102050/1‑141 (MQ=255)                                           tgtgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCaaaa                                                                                                              >  2:608152/1‑149 (MQ=255)                                             tgCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCa                                                                                                           >  1:346201/1‑150 (MQ=255)                                              gCGGATTCACTGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAc                                                                                                          <  2:556486/150‑1 (MQ=255)                                              gCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAc                                                                                                          <  2:758492/150‑1 (MQ=255)                                               cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGc                                                                                                        <  2:140189/151‑1 (MQ=255)                                                 gATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACGctc                                                                                                      >  1:144926/1‑151 (MQ=255)                                                  aTTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAACACGctct                                                                                                     >  2:149118/1‑151 (MQ=255)                                                   ttCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTTGCCCCCCACTTTCTCGATGCAAAAGCACGccct                                                                                                     >  1:869526/1‑150 (MQ=255)                                                     cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGt                                                                                                  <  2:901949/151‑1 (MQ=255)                                                     cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGTGTTTAGCCCCCCACTTACTCGATGGAAAAGGACGCTCTgg                                                                                                   >  1:247302/1‑150 (MQ=255)                                                      aCGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGtt                                                                                                 <  1:895416/151‑1 (MQ=255)                                                         gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATc                                                                                              <  1:848364/151‑1 (MQ=255)                                                           cAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTTGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCg                                                                                             >  1:778733/1‑150 (MQ=255)                                                                  gTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGgag                                                                                    <  2:339000/151‑1 (MQ=255)                                                                          ttCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGTGTTTTGTCCCCTCCTTACTCGATGCAAAAGCACGCTCTGGTTTTTGTC‑TTGGAGACgg                                                                                >  1:108160/1‑147 (MQ=255)                                                                            cGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAgg                                                                           <  2:129262/150‑1 (MQ=255)                                                                                tCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGGTATCGTC‑TTGGAGACGGACAGGgtt                                                                        >  1:622133/1‑149 (MQ=255)                                                                                   tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATTAGGTTTGCCCCCTACTTACTCCATGCAAAAACACGCTCCGGGTTTTGTCTTTGGGAACGGACGGGtt                                                                         >  2:880882/1‑146 (MQ=255)                                                                                    aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCCACTTACTCGATGCAAAAGAACGCCCTGGTTATCGTC‑TTGGAGACCGACAGGGTTATGACg                                                                  >  2:184555/1‑151 (MQ=255)                                                                                     aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGt                                                                 <  2:1014540/151‑1 (MQ=255)                                                                                           tgGCGCTTGCCGACCAGGCCATTCAGGTAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCTGGCTGCTTACGCGATGCAAAAGCACGCTCTGGTTAACGTC‑TTGGAGACGGACAGGTTTATGACGTAATcc                                                            <  2:475675/150‑1 (MQ=255)                                                                                             gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGGTTAGCCCCCCACTTACTCGCTGCAAAAGCACGCTCTGGTTATCGTC‑TTTGAGACGGACAGGCTTATGACGTAATcc                                                            >  2:158817/1‑148 (MQ=255)                                                                                               gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTCCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                       <  2:435251/151‑1 (MQ=255)                                                                                               gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTgcgc                                                       <  1:936607/151‑1 (MQ=255)                                                                                               gcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTgcg                                                        >  1:291530/1‑150 (MQ=255)                                                                                                 ttGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCCCCTTACTCGATGGAAAAGCACGCCCTGGGTTTTGTC‑TTGGGAACGGACAGGGTTTTGACGTAATCCTGcgcg                                                      >  2:31417/1‑150 (MQ=255)                                                                                                  tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGGAAAAACACGCTCTGGGTTTTGTC‑TTGGAGACGGACAGGGTTATAACGTAATACTGcgcg                                                      >  1:994205/1‑149 (MQ=255)                                                                                                   gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATg                                                   <  2:984493/151‑1 (MQ=255)                                                                                                    ccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAAACCTGCGCGATgg                                                  >  1:26275/1‑151 (MQ=255)                                                                                                      gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGcc                                                <  1:949638/151‑1 (MQ=255)                                                                                                       cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGcc                                                >  1:479960/1‑150 (MQ=255)                                                                                                              cATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGTGGGGATGGGAAATATAAGGTTAACCCCCCACTTTCTCGATGCAAAAACACGCTCTGGTTATCGTC‑TTTGGGACGGACAGGTTTATGACGTAAACCTGcgcggagggct                                               >  1:602281/1‑137 (MQ=255)                                                                                                                ttCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGGGTTTAGCCCCCCACTTAATCGATGCAAAAGGACGCTCTTGTTATCGTC‑TTTGAGACGGACAGGGTTATGACGTAATACTGCGCGATGGGCt                                               >  2:174121/1‑142 (MQ=255)                                                                                                                  caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCa                                    <  2:229445/151‑1 (MQ=255)                                                                                                                   aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAGCCCCCCACTTACTCGAAGCCAAAGGACGCTCTGGGTATCCTC‑TTTGGGGCGGGAAGGTTTATGACGCAAACCTGTCCCGTGGCCTGGCGTGCGGCa                                    >  2:316473/1‑150 (MQ=255)                                                                                                                         gTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATgg                             <  1:184565/151‑1 (MQ=255)                                                                                                                              gcgcAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATc                        <  2:473214/151‑1 (MQ=255)                                                                                                                                   gagCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAAAACGCTCTGGTTATCTTC‑TTGGGGACGGACAGGGTTTTGACGTAATCCTTCGCGGTGGGCTGATGTGCGGCACCCAAGGTTTTCata                     >  2:322929/1‑149 (MQ=255)                                                                                                                                          gtgGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGAtt             <  2:583373/150‑1 (MQ=255)                                                                                                                                               ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc       <  2:247289/151‑1 (MQ=255)                                                                                                                                               ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc       >  1:942732/1‑151 (MQ=255)                                                                                                                                                    gTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGaaaa  <  2:920660/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CGGCGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTC‑TTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCGAAAA  >  NC_000913/2326302‑2326599

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