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

CGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCG  >  NC_000913/2326305‑2326595                                                                                                                                                  |                                                                                                                                                  cgcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGa                                                                                                                                              <  1:197410/151‑1 (MQ=255)  gcgTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGACAGGCGCAGAGCATTGTGGCGGGGGGGATgg                                                                                                                                               >  1:2689/1‑149 (MQ=255)      cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                                          <  2:505494/150‑1 (MQ=255)      cAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAat                                                                                                                                          <  2:448519/150‑1 (MQ=255)        ggCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATg                                                                                                                                       <  2:397021/151‑1 (MQ=255)          cACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATg                                                                                                                                       >  2:529648/1‑149 (MQ=255)              gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGtt                                                                                                                                   <  2:833537/149‑1 (MQ=255)              gTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAggt                                                                                                                                   >  1:502563/1‑147 (MQ=255)               ttAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGAAAAAATGAGttt                                                                                                                                  >  1:1090436/1‑149 (MQ=255)                  aaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATTTGTTTTGccc                                                                                                                             >  1:25411/1‑151 (MQ=255)                  aaaaGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAAAATGAGGTTAGccc                                                                                                                             >  1:538342/1‑151 (MQ=255)                            tGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGGTTAACCCCCtacttac                                                                                                                    >  1:828144/1‑150 (MQ=255)                               cAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGGGTTTTGACCCCtaatt                                                                                                                      >  2:533190/1‑145 (MQ=255)                                aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGAt                                                                                                               <  2:207391/151‑1 (MQ=255)                                aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGAt                                                                                                               <  2:852330/151‑1 (MQ=255)                                aGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAAAATGAGGTTTGCCCCCTACTTACTc                                                                                                                  >  2:552266/1‑148 (MQ=255)                                 gAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATg                                                                                                              <  1:1130387/151‑1 (MQ=255)                                   aaCGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGc                                                                                                             >  2:367950/1‑150 (MQ=255)                                     cGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaa                                                                                                          <  2:538370/151‑1 (MQ=255)                                     cGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaa                                                                                                          <  1:775335/151‑1 (MQ=255)                                      ggTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCaaaa                                                                                                         <  2:527412/151‑1 (MQ=255)                                            cGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGTGGGAAATATGAGTTTAGCCCCCCACTTACTCGATGGAAAAGCa                                                                                                      >  1:460217/1‑148 (MQ=255)                                             ggATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGGCCCCTACTTACTCGATGCAAAAGCACGcc                                                                                                  >  2:1198034/1‑150 (MQ=255)                                             ggATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCCACTTACTCGATGCAAAAACACg                                                                                                    >  2:1085452/1‑149 (MQ=255)                                                 tCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTg                                                                                               <  2:212272/150‑1 (MQ=255)                                                  cACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTgg                                                                                              <  1:1017530/150‑1 (MQ=255)                                                  cACGGTCAATAAAGTATGTGGTGCGGGTCTTAAAAGTGGGGCGCTTGCCGCCCAGGGCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAAGATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGt                                                                                             <  1:651298/151‑1 (MQ=255)                                                    cGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAAAACGCTCTgggt                                                                                            >  2:832574/1‑148 (MQ=255)                                                      gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATc                                                                                         <  2:130052/151‑1 (MQ=255)                                                      gTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCAACCTCTGGtt                                                                                            >  2:1282887/1‑148 (MQ=255)                                                          aTAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCt                                                                                     <  2:958851/151‑1 (MQ=255)                                                              aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGga                                                                                 <  2:701392/151‑1 (MQ=255)                                                              aGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGga                                                                                 <  2:141524/151‑1 (MQ=255)                                                                tATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGgag                                                                                <  1:631043/150‑1 (MQ=255)                                                                 aTGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGGCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGAc                                                                              <  2:6476/151‑1 (MQ=255)                                                                   gtgGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACgg                                                                            <  1:408061/151‑1 (MQ=255)                                                                      gTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGGTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACgg                                                                            >  2:283195/1‑148 (MQ=255)                                                                          gggTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAAAGAGTTTAGCCCCCCACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAgg                                                                       >  1:374782/1‑149 (MQ=255)                                                                                tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTAt                                                                  <  2:208264/148‑1 (MQ=255)                                                                                tAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTTGAGACGGGCAGGGTTAt                                                                  >  1:1029268/1‑148 (MQ=255)                                                                                 aaaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGAAAGGTTTATTACg                                                              >  2:1104933/1‑151 (MQ=255)                                                                                  aaaGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTTCTCCATGCAAAAGCACCCTCTGGGTTTTGTCTTTGAGACGGACAGGTTTTTGaag                                                              >  1:1288481/1‑148 (MQ=255)                                                                                    aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTaa                                                           <  1:22963/151‑1 (MQ=255)                                                                                    aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGTTGGAAAATATTAGGTTAGGCCCCCACTTATTCCGTGCGAAAGCACGCTCTGGGTAACGTTTTTGAGACGGAGAGGGTTATGa                                                                >  2:442521/1‑146 (MQ=255)                                                                                    aGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGAAGGGAAAGATAAGGTTTACCCCCCACCTACTCGAATCCAAAACACGCTCTGGTTATCGGCTTTGGGAGGGGCAGGGTTATTACgga                                                            >  1:171648/1‑148 (MQ=255)                                                                                         ggCGCTTGCCGCCCAGGCCGTTCAGGTAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCTCTGCTTACTCGATGCAAAAGCACGCTCTGGTTATCGTATTGGATACGGACAGGTTTATGACGTAATCCTg                                                      <  1:1217204/151‑1 (MQ=255)                                                                                         ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGTCAGGTTTATGACGTAATCCTg                                                      <  2:1255190/151‑1 (MQ=255)                                                                                         ggCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCt                                                       <  2:861904/150‑1 (MQ=255)                                                                                          gcgcTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTgc                                                     <  2:994348/151‑1 (MQ=255)                                                                                               tGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGCAAACTGAGTTTAAAACCCCCCCTACTCGCTGACAAAGAACGCTCTGGGTATCGTCTTTGAGACGGGAAGGGTTTAGACGTAATACTGcgcg                                                  >  1:728681/1‑149 (MQ=255)                                                                                                gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATg                                               <  1:763188/151‑1 (MQ=255)                                                                                                gccgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAAAATGAGGTTAACCCCCCACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGAGACGGACAGGGTTATGACGTAATCCTGcgcg                                                  >  2:1071456/1‑148 (MQ=255)                                                                                                  cgccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGc                                             <  1:275496/151‑1 (MQ=255)                                                                                                   gccCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGGCTGAAAATATTAGTTTAGCCCCCTACTTACTGGAGGCCAAAACACGCCCTGGTTATCGTCTTGGGGCGCGGCATGGTTTTGACGCAAACCCGCGCGATGGcc                                            >  1:523832/1‑151 (MQ=255)                                                                                                    cccAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCt                                           <  2:544335/151‑1 (MQ=255)                                                                                                      cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTTGTTATCGTCTTGGGGACGGACAGGGTTATTACGTAATCCTGCGCGATGGCCTg                                          >  2:175593/1‑150 (MQ=255)                                                                                                      cAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGTTATCGTCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGATGGcc                                            >  2:342473/1‑148 (MQ=255)                                                                                                         gCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAtgt                                      <  2:1153107/151‑1 (MQ=255)                                                                                                            aTTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACACGCTCTGGGTATCGTCTTGGAGACGGAAAGGGTTATTACGTAATCCTGCGCGATGGCCCGAtg                                       >  2:197573/1‑147 (MQ=255)                                                                                                              tCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAACGCGGTCTGGTTATCGTCTTGGAGAGGGAGAGGGTTATGTTGTAATACTGTGCGATGGCCTtgtgtgg                                    >  1:860937/1‑147 (MQ=255)                                                                                                               caggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGGAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGcc                                 >  2:1238741/1‑150 (MQ=255)                                                                                                                aggcaggTCAGGCGCAGAGCATTGTGGCGGGGGGGATGGAAAATATGAGTTTAGCCCCCTACTTAATCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGGGACGGACAGGTTTTTGACGCAATCCTGCGCGATGGCCTGGTGGGGGcc                                 >  2:451387/1‑149 (MQ=255)                                                                                                                         aGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTa                      <  2:553820/151‑1 (MQ=255)                                                                                                                               agagCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCTCCCATGGTTATCata                 <  1:1282846/150‑1 (MQ=255)                                                                                                                                  gCATTGTGGCGGGGGGGATGGAAAAAATGAGTTTAGCCCCCTCCTTACTCGATGCAAAAGCACGCTCTGGGTATCGTCTTGGGGACGGACAGGGTTATGACGTAATCCTGCGCGACGGCCTGATGTTCGCCCCCCATGGTTATCATATggg             >  1:884726/1‑151 (MQ=255)                                                                                                                                          gCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGAAGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccg     <  1:1290929/151‑1 (MQ=255)                                                                                                                                            ggggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGAGGTGCGCCACCCATGGTTATCATATGGGGATTAccgcc   >  1:603465/1‑151 (MQ=255)                                                                                                                                             gggggTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGAGGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTAccgccg  <  1:763939/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CGCGTCAGGCACTGTTAAAAAGCGGGCTGGCAGAAACGGTGTGCGGATTCACGGTCAATAAAGTATGTGGTTCGGGTCTTAAAAGTGTGGCGCTTGCCGCCCAGGCCATTCAGGCAGGTCAGGCGCAGAGCATTGTGGCGGGGGGTATGGAAAATATGAGTTTAGCCCCCTACTTACTCGATGCAAAAGCACGCTCTGGTTATCGTCTTGGAGACGGACAGGTTTATGACGTAATCCTGCGCGATGGCCTGATGTGCGCCACCCATGGTTATCATATGGGGATTACCGCCG  >  NC_000913/2326305‑2326595

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