Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,101	0	G	T	44.4%	118.0 / 42.5	72	T222N (ACC→AAC)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base G (33/7);  new base T (0/32);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 6.82e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537952‑4538249                                                                                                                                                          |                                                                                                                                                       cATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAgg                                                                                                                                                        <  1:353631/150‑1 (MQ=255)  aTTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGtt                                                                                                                                                      >  1:416783/1‑151 (MQ=255)  aTTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGtt                                                                                                                                                      >  1:267083/1‑151 (MQ=255)   ttGCTTATCCCCGTAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CTGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGttt                                                                                                                                                     >  2:443412/1‑151 (MQ=255)     gCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGtttt                                                                                                                                                    <  1:425617/150‑1 (MQ=255)        tATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATc                                                                                                                                                 >  2:275688/1‑150 (MQ=255)         aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCa                                                                                                                                                <  2:626175/150‑1 (MQ=255)                 aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑c                                                                                                                                     >  1:104935/1‑151 (MQ=255)                  aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑c                                                                                                                                     >  1:170267/1‑150 (MQ=255)                        ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAAccg                                                                                                                              >  1:530760/1‑151 (MQ=255)                             tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGCTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACtgccgc                                                                                                                          >  1:154171/1‑150 (MQ=255)                                 ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTACATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCAcc                                                                                                                       >  2:397076/1‑149 (MQ=255)                                 ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCAcc                                                                                                                       >  2:332419/1‑149 (MQ=255)                                  cATCTGTTAATGAGACGGAAGCAACCTTATCCCTTTTAAAATTATTCCCGGGAAAATCAATTTAAATTA‑GGGCAC‑CCGTTCGCAATCCGGGTTGGC‑CTCCCCCATTAAAAAG‑CCAGTTTTAATCCCTTTT‑ATT‑CACAACCGCCGCACCAGc                                                                                                                    <  2:273976/151‑1 (MQ=255)                                  cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCCCACCAGc                                                                                                                    >  2:79584/1‑151 (MQ=255)                                        gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTcc                                                                                                              >  1:253223/1‑151 (MQ=255)                                             gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                          >  2:199588/1‑150 (MQ=255)                                             gagaCGGGAGAAACCTTTTTCATTTTAAAATTAATACCGGTAAATTAAAGTAAAAATA‑CGGTAT‑CTGTTCACAATCCTGTTTTGC‑TTCCGCCATAAAAAAG‑CCGGTTTTTTCCTCTTTT‑TTT‑AACACCCGCCGCACCAGCCGTTCCGTa                                                                                                           <  2:477427/149‑1 (MQ=255)                                                  gggAGCAAGTTTTTTCCATTTAAAATTATTCCCGGTAAAATAAATTTAAAATC‑CGGCTC‑CCTTTCGCAACCCTGTTTTGG‑TTTCCCCTTAAATAAG‑CCAGGTTTTATCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAggg                                                                                                     <  1:731851/150‑1 (MQ=255)                                                   ggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                   >  2:643006/1‑151 (MQ=255)                                                   ggAGAAAGTTTATTCCATTTAAATTTATTACCGGAAAAATAAATTTAAAATA‑CGGCAC‑CCGTTCGCAATCCTGGTTTGG‑CTTCCCCATTAATAAG‑CCGGTTTTTATCCCTTTT‑ATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                   <  1:830727/151‑1 (MQ=255)                                                    gAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                   >  1:395354/1‑150 (MQ=255)                                                         aGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                             >  1:879893/1‑151 (MQ=255)                                                          gttttttCCATTTTAAATTTTTCCCGATAAAAT‑AAATTAAAACACCGACAC‑CCGTTCGCATTCGTGTTTTGG‑TTGACCCATTAAT‑CG‑CCCGTTTTTTGCCCCTTTAATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                              <  2:284335/144‑1 (MQ=255)                                                          gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  1:585335/1‑151 (MQ=255)                                                          gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  2:427224/1‑151 (MQ=255)                                                            ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                            >  1:1003299/1‑149 (MQ=255)                                                              aTTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CCTCAT‑CCGTTAGCAATTCTGGCTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGATCCGTACCAGGGCGATTCGCCAGCg                                                                                        >  2:146664/1‑151 (MQ=255)                                                                tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                       >  2:926221/1‑150 (MQ=255)                                                                  cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAAc                                                                                    >  2:291703/1‑151 (MQ=255)                                                                     tttAAAATTTTACCGGTTAAAATAAGTTAAAAAC‑CGGCCTGCCTTCCGAAATCTGGGTTGGT‑CTTCGCCATAAAT‑AG‑GCAGGTTTTGCCCCCTTT‑TTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                   <  2:277875/150‑1 (MQ=255)                                                                        aaaTTATTCCCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCCCCATAAATAAG‑CCAGTTTTTATCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                              <  1:992296/151‑1 (MQ=255)                                                                        aaaTTATTCCCGGTGAAATAAATTTAAAATA‑CGCCAC‑CCTTCCGCAATCCTGTTTTGT‑CTTCCCCATTAATAAG‑CAAGTTTTTATCCCTTTT‑ATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                 <  1:610123/148‑1 (MQ=255)                                                                        aaaTTATTCCCGGTAAAATCAATTTAAAATA‑CGGCTC‑CCGTCCCAAACCCTGTTTGGG‑CTTCCCCTTAAAAAAC‑CCAGTTTTTACCCCTTTT‑TTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                                <  1:964841/149‑1 (MQ=255)                                                                                aCCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                      >  1:505174/1‑151 (MQ=255)                                                                                 cctggaaaaTCAATTTAAAATA‑CGCCA‑‑CCCTTCCCAATCCTTGTTTTG‑GTTCCGCAATAAAAAGCCCAGTTTTTACCCCCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                      <  2:168560/144‑1 (MQ=255)                                                                                 ccGGTAAAATAAAGTTCAAATC‑CGGCCC‑CCGTCCGCAACCCTGTTTGGG‑CTCCCCCATTAATAAC‑CCGGTTTTTACCCCCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                     <  1:144441/151‑1 (MQ=255)                                                                                 ccGGTAAAACAAAGTTAAAATC‑CGCCAC‑CCGTCCGAATTCCGGTTTTGG‑CTCCCCATTTAAAAAG‑CCGGTTTTTATCCCTTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                       <  1:336688/149‑1 (MQ=255)                                                                                    gaaaaacaaattacAAATC‑CGCCTT‑CCGTCCCC‑ATCCCGGTTTGGCCTCCCCCATAAAAAAC‑CCGGTTTTTCTCCCTTTT‑ACC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                       <  2:70658/133‑1 (MQ=255)                                                                                       aaataaatTTCAAATA‑CGGCAC‑CCTTTCGCAACCCTGTTTTGG‑CTTCCCCATAAATAAG‑CCAGGTTTTATCCCCTTT‑ATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                 <  1:808844/149‑1 (MQ=255)                                                                                       aaaTCAATTTAAATTC‑GGGCAC‑CCTTTCGAAATCTTGTTTTGT‑TTTCGCCTTAAAAAAG‑CCAGTTTTTATCCCTTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                  <  1:947050/148‑1 (MQ=255)                                                                                       aaaTCAAGTTAAAATA‑GGCCAC‑CCGTCCCCAACCTTGTTTGGG‑CTTCCCTATAAAAAAG‑CCAGTTTTTCTCCCTTTA‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa                                                                    <  2:757566/146‑1 (MQ=255)                                                                                                  aaaaA‑CGGCAC‑CCGTTCGCAATCCTGTTTTGT‑CTCCCCCATTAATAAG‑CCGGTTTTTATCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                       <  1:682053/148‑1 (MQ=255)                                                                                                   aaTA‑GGGCAT‑CTGTGCGCATTGCGTGTTGGT‑CTGCCCTATAAATAAG‑CCTGTTTTTACCCCTTTT‑CTT‑AAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaaca                                                   <  2:782446/151‑1 (MQ=255)                                                                                                      a‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                 >  2:743077/1‑150 (MQ=255)                                                                                                         ggCAT‑CCGTTCGCAACCCTGTTTTGG‑CTTCGCCTTAAAAAAC‑CCAGTTTTTATCCCTTTT‑ATT‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                              <  1:66274/151‑1 (MQ=255)                                                                                                            cacccGTTCACACTCGTTTTTTGG‑TTTCCCCATTAATAAG‑CCAGTTTTTATCCCTTTT‑ATT‑CCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                               <  1:717083/145‑1 (MQ=255)                                                                                                               cctttcccatcccTGGTTTGG‑CTTCCCCATTAATAAC‑CCGGTTTTTATCCCCTTT‑ATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                            <  2:416194/137‑1 (MQ=255)                                                                                                                  ttCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       >  1:681749/1‑150 (MQ=255)                                                                                                                    cGCCATCCTTGTTTGT‑CTTCCCCTTAAAAAAC‑AGAGTTTTTACCTCTTTT‑CTT‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                     <  2:509154/150‑1 (MQ=255)                                                                                                                    cGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                       >  1:451184/1‑148 (MQ=255)                                                                                                                    cGCAATCCGGTTTTGG‑CTTCCCCATTAAAAAG‑CCAGTTTTTATCCCCTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                    <  1:157828/151‑1 (MQ=255)                                                                                                                       aaTCCTGTTTTGG‑TTTCCCCATTAAAAAG‑CCGGTTTTTATCCCCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 <  1:789457/151‑1 (MQ=255)                                                                                                                        aTCCTGTTTTGC‑CCTCCCTATTAATACC‑CCAGTTTTTATCCCTTTA‑ATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                  <  2:569578/149‑1 (MQ=255)                                                                                                                          ccTGGTTTGG‑TTTCGCCATAAAAAAG‑CCGGTTTTTATCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                 <  1:963886/148‑1 (MQ=255)                                                                                                                                ttggcctccGCCATTAAAAGG‑CCAGTTTTTACCCCTTTT‑TTC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                          <  2:759409/142‑1 (MQ=255)                                                                                                                                ttggccTTCCCCTTAAATAAC‑CCGGTTTTTACCCCTTTT‑ATA‑CAAATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                               <  2:254765/140‑1 (MQ=255)                                                                                                                                ttGG‑CTTCGCCATAAATAAG‑CCGGTTTTTACCCCCTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                         <  1:351009/150‑1 (MQ=255)                                                                                                                                ttGG‑CTCCCCCATTAATAAC‑CCAGTTTTTATCCCTTTT‑ATC‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                         <  2:64677/150‑1 (MQ=255)                                                                                                                                ttGC‑CTTCCCCATAAATAAG‑CCATGTTTTACCCCCTTT‑ATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                         <  2:928270/150‑1 (MQ=255)                                                                                                                                 tGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                       >  1:526658/1‑151 (MQ=255)                                                                                                                                        cGCCATAAAAAAC‑CCGGTTTTTACCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                  <  1:643276/150‑1 (MQ=255)                                                                                                                                          ccATTAATAGG‑CAGGGTTTAATCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                   <  2:253054/147‑1 (MQ=255)                                                                                                                                             ttAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa             >  1:271354/1‑150 (MQ=255)                                                                                                                                             ttAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa             >  1:748641/1‑150 (MQ=255)                                                                                                                                              ttaaaagcccAGTTTTTATCCCCTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa             <  1:376267/142‑1 (MQ=255)                                                                                                                                               aataaC‑CCAGTTTTTATCACTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa           <  1:500195/150‑1 (MQ=255)                                                                                                                                                 taaG‑CAAGTTTTTATCCCTTTT‑ATT‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt         <  1:290076/150‑1 (MQ=255)                                                                                                                                                  aaG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGAGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt         >  1:384331/1‑149 (MQ=255)                                                                                                                                                   aG‑CCAGGTTTTATCACCTCT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt      >  1:605172/1‑151 (MQ=255)                                                                                                                                                      ccAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg     >  2:701927/1‑150 (MQ=255)                                                                                                                                                       cAGTTTTTATCCCTTTT‑ATC‑CACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg     <  2:151824/149‑1 (MQ=255)                                                                                                                                                        aGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc  >  2:178537/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       CATTGCTTATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGC  >  NC_000913/4537952‑4538249

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