Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	38.2%	128.1 / 17.9	71	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (37/5);  new base C (0/26);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 3.88e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTT  >  NC_000913/4537962‑4538257                                                                                                                                                         |                                                                                                                                                       ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCt                                                                                                                                                     <  1:1411282/151‑1 (MQ=255)  cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCtt                                                                                                                                                    >  2:1384178/1‑151 (MQ=255)  cccGAAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GT‑TTTTATCACCtt                                                                                                                                                    <  2:315113/151‑1 (MQ=255)     gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTAt                                                                                                                                                 <  2:81320/151‑1 (MQ=255)          aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑cac                                                                                                                                            >  2:184596/1‑150 (MQ=255)            cGCCAGCTACGCCATCTGGTGATGTGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CGTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑c                                                                                                                                        >  2:41932/1‑151 (MQ=255)            cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACaa                                                                                                                                          <  2:1181251/150‑1 (MQ=255)               cAGCTACGACATCTGGTGATGAGACGGGAGCAAGCGTATTCCAGTTTAAATTATTACCGTTGAAATCAAGTTCAAATACGTCA‑TCCTTTCGCAATCCTGA‑TTCGG‑CTTCGCCATTAATAAGCCA‑GG‑GTTTATCTCCTTTAAT‑CACAA‑c                                                                                                                                        <  2:1269226/148‑1 (MQ=255)                    aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCAc                                                                                                                                >  1:785072/1‑151 (MQ=255)                    aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCAc                                                                                                                                >  2:401551/1‑151 (MQ=255)                      gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCAcc                                                                                                                               >  2:204033/1‑150 (MQ=255)                          tCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCg                                                                                                                          >  1:673567/1‑151 (MQ=255)                          tCTGGGGAAGAGACGGGAGAAAGCTTTTTCCATTTTAAATTATTATCGTAGCACTCAAGTTCAAAAAAAGCA‑TCTGTCCGCAATCCTGC‑TTTGG‑CTTCGCCATTAAGAAGCCT‑GG‑TTTCAACACCTTTATT‑CACAA‑CCGCCGCAcc                                                                                                                               >  2:717083/1‑146 (MQ=255)                                     gaCGGGAGCACGCTTACTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAACTACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTAc                                                                                                                  >  2:951146/1‑148 (MQ=255)                                      aCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                              >  1:374894/1‑151 (MQ=255)                                       cGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                             >  2:1020859/1‑151 (MQ=255)                                        gggAGCAATCTTTTTCCTTTTAAAATTATTCCCGGTAAAATAAATTTAAAATACGGCT‑CCCGTTCGCAATCCTGT‑TTGGG‑CTTCCCCATTAATAACCCA‑GT‑TTTTATCCCCTTTATC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                            <  2:1323646/151‑1 (MQ=255)                                            gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                        >  1:1058699/1‑151 (MQ=255)                                            gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                        >  2:868156/1‑151 (MQ=255)                                                 cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                    >  1:355644/1‑150 (MQ=255)                                                  ttttccctttttAAATTATAACGGGTGAAATCAATTTAAAATCGGCCA‑TCGGTTCCAAGTCCGGT‑TTGGG‑CTTCGCTATAAATAACCCA‑GT‑TTTTATCCCCTTTATC‑AACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                                   <  1:745634/142‑1 (MQ=255)                                                     ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                               >  2:637837/1‑151 (MQ=255)                                                     ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                               >  1:504963/1‑151 (MQ=255)                                                        cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                            >  1:909455/1‑151 (MQ=255)                                                        cATTTTAAAATTTTTCCGGGGAAAATAAATTCAAAATCCGCACCCCTTTCCCATTCCTGT‑TTTGG‑CTTCCCCATTAAAAAGCCG‑GT‑TTTTATCCCCTTTTTC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                                 <  1:485341/147‑1 (MQ=255)                                                           tttAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                         >  2:983639/1‑151 (MQ=255)                                                             tAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                       >  1:392895/1‑151 (MQ=255)                                                              aaaTTATTCCCGGTAAAATAAATTTAAAATCCGGCA‑TCCGTCCACAATCTTGT‑TTTGG‑TTTCGCAATAAATAAGCCG‑GT‑TTTTATCCCCTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                      <  1:90075/151‑1 (MQ=255)                                                               attttttcccGGTAAAATCAAGTGCAAATCGGCCC‑GCCGTGCACAA‑CCCGG‑TTTGGCCGCCCCCATTAAGAGCCCG‑GG‑TTTTATCACTTTTATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                       <  1:1128881/141‑1 (MQ=255)                                                                 ttgttcccGGTAAAATGAATTAAAAATCCGCCA‑TCTGTTCGCACTCGTTG‑TTGGG‑C‑TCGCCCATAATGAGGCATTG‑TTTTATC‑CCTTTTTC‑CCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                      <  2:1051785/142‑1 (MQ=255)                                                                  tattaCCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                  >  1:101333/1‑151 (MQ=255)                                                                  tattaCCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                  >  2:1369563/1‑151 (MQ=255)                                                                    ttcccGGTAAAATCAGTTTCAAAACCGCCT‑CCCTTTCACACTCTGGT‑TTTGC‑CTGCCCAATTAATAAGCCT‑GT‑TTTTACCCCTTTTTTC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgt                                                                                 <  2:414580/147‑1 (MQ=255)                                                                     taCCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                               >  1:934905/1‑151 (MQ=255)                                                                       ccGGTAAAATCAATTTCAAATACGGCA‑CCCGTTCGCAATCCTGT‑TTTGG‑CTTCCCCATAAAAAACCCA‑GT‑TTTTATCCCCTTTTTT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                                <  1:195468/148‑1 (MQ=255)                                                                         gggaaaaTCAATTCAAAAAACGAAA‑TCCGTTCGCAA‑CCTGGTTTTGG‑TTTCCCTATAAATAGCCGA‑GG‑TTTTCCCTCTTTAATC‑CCAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                               <  1:471058/143‑1 (MQ=255)                                                                            gAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                         >  1:1026619/1‑150 (MQ=255)                                                                             aaaTCAATTAAAAATAGGCCA‑CCCTTTCGAACCCCTGT‑TTTGG‑‑TTCCCCTTTAAAAAACCA‑GGTTTTTACCCCTTTTATC‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                       <  1:664644/151‑1 (MQ=255)                                                                              aaTCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                      >  1:456654/1‑151 (MQ=255)                                                                                  aatttaaaaTCCGGCA‑TCCGTCCGCAATCCTGT‑TTGGT‑CTTCGCCATAAATAAGCCA‑GT‑TTTTATCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                    <  1:510730/143‑1 (MQ=255)                                                                                  aaGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                   >  2:179590/1‑150 (MQ=255)                                                                                             cGGCA‑TCCTTTCCCACTCGTGG‑TTTGC‑CTCCCCCATTAAAAACCCT‑GT‑TTTTACCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGgaa                                                       <  2:640828/151‑1 (MQ=255)                                                                                                cA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                    >  1:9707/1‑151 (MQ=255)                                                                                                  cacctgttccaaatcTGG‑TTTTG‑TTTCCCCCTTAAAAAGCCA‑GG‑TTTTTACCCCTTTTTTCCACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                      <  1:204083/135‑1 (MQ=255)                                                                                                    ccttccGCAATCTTGT‑TGTGG‑CTCCCCAATAAAAAACCCA‑GT‑TTTTCCCCCTTTTATC‑CTCAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                       <  1:1031905/140‑1 (MQ=255)                                                                                                      gtccGCAACCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCG‑GT‑TTTTACCCCCTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                                <  1:1373014/147‑1 (MQ=255)                                                                                                       ttCGAAATCCTGT‑TTTGC‑TTCCGCAATTAAAAACCAA‑GT‑TTTTACCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                               <  1:1200707/150‑1 (MQ=255)                                                                                                        tCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                               >  2:11358/1‑149 (MQ=255)                                                                                                           cAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                           >  1:1169165/1‑150 (MQ=255)                                                                                                            aaTCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                         >  2:900622/1‑151 (MQ=255)                                                                                                             actcttg‑TTTGGTCTTCGCTATAAAAAAGCCT‑GT‑TTTTACCCCCTTTATT‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                           <  2:887107/143‑1 (MQ=255)                                                                                                               cgtGG‑TTTGG‑CTTCGCCATTACTAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                        >  2:1027038/3‑149 (MQ=255)                                                                                                                   gttttGC‑CTTCCCATTTAATAACCCG‑GT‑TTTTATCCCCTTTATC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                    <  1:401226/148‑1 (MQ=255)                                                                                                                     tttGG‑TTTCGCCATTAAAAAGCCG‑GT‑TTTTACCCCTTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                  <  1:1307888/150‑1 (MQ=255)                                                                                                                       tGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                               >  2:1097647/1‑151 (MQ=255)                                                                                                                       tGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                               >  2:97071/1‑151 (MQ=255)                                                                                                                         gcctgctcaaTTAATAACCAT‑GG‑TTTTATCCCTTTTATT‑CCCAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                               <  1:296190/141‑1 (MQ=255)                                                                                                                           cTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                            >  1:443356/1‑151 (MQ=255)                                                                                                                                ccATAAATAAGCCG‑GT‑TTTTATCCCCTTTATA‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                        <  1:531342/150‑1 (MQ=255)                                                                                                                                   taaataaGCCG‑GT‑TTTTATCCCCTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                     <  1:595829/148‑1 (MQ=255)                                                                                                                                     taaatacacttg‑TTTTTTCCCCTTTTTC‑CACAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                         <  2:1319774/135‑1 (MQ=255)                                                                                                                                       aaacccA‑GT‑TTTTATCCCTTTTATC‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                <  2:355561/147‑1 (MQ=255)                                                                                                                                        agaccA‑GG‑TTTAACCCCCTTTATA‑CAAAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                 <  1:1328194/146‑1 (MQ=255)                                                                                                                                        aaGCCA‑GT‑TTTTATCCCTTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg               <  1:1042884/151‑1 (MQ=255)                                                                                                                                           ccA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg             >  1:936405/1‑150 (MQ=255)                                                                                                                                           ccA‑GG‑TTTTACCCCTTTTATC‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa            <  1:24604/151‑1 (MQ=255)                                                                                                                                                  ttttATCCCTTTTATT‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt        <  1:401642/150‑1 (MQ=255)                                                                                                                                                   tttATCCCCTTTATT‑CACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc          <  2:384602/147‑1 (MQ=255)                                                                                                                                                    ttATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg      >  1:743684/1‑150 (MQ=255)                                                                                                                                                      aTCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt   >  1:1068631/1‑151 (MQ=255)                                                                                                                                                       tCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt  >  1:1103858/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCA‑TCCGTTCGCAATCCTGG‑TTTGG‑CTTCGCCATTAATAAGCCA‑GG‑TTTTATCACCTTTATT‑CACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTT  >  NC_000913/4537962‑4538257

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: ‑