Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	51.4%	80.9 / 124.2	74	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (34/1);  new base C (0/37);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 8.83e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-13

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                          |                                                                                                                                                        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                                                                                                     >  2:427224/1‑151 (MQ=255)                                                  gCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                     >  1:585335/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)                                                                                                                                      ttaaaaacccgttttttATCCCTTTT‑ATT‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                      <  2:530508/137‑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)                                                                                                                                                 gtttttACCCCTTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:177851/149‑1 (MQ=255)                                                                                                                                                  tttttACCCCTTTT‑ATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa          <  1:408867/149‑1 (MQ=255)                                                                                                                                                    tttATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  1:549900/1‑151 (MQ=255)                                                                                                                                                    tttATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg       >  2:257702/1‑151 (MQ=255)                                                                                                                                                        tCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAttt   >  2:978633/1‑151 (MQ=255)                                                                                                                                                         cACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:438264/1‑151 (MQ=255)                                                                                                                                                         cACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:100202/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                        ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATA‑CGGCAT‑CCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAG‑CCAGGTTTTATCACCTTT‑ATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258

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