Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	44.8%	107.1 / 17.5	67	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (29/8);  new base C (0/30);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 5.24e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  NC_000913/4537960‑4538258                                                                                                                                                           |                                                                                                                                                       aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCAc                                                                                                                                                       >  2:290891/1‑151 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑t                                                                                                                                                    <  2:250983/151‑1 (MQ=255)    cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑tt                                                                                                                                                   <  1:892122/151‑1 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTa                                                                                                                                                 <  1:744164/151‑1 (MQ=255)        aaaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTAtt                                                                                                                                               >  1:42790/2‑151 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTc                                                                                                                                              >  2:519211/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTca                                                                                                                                             <  2:541602/151‑1 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTca                                                                                                                                             <  2:613671/151‑1 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTc                                                                                                                                              >  2:493523/1‑150 (MQ=255)           aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAGGATAATTCCTTTTTAAATTATTA‑CC‑GGGGAAACCAAGTCCAATTAGGCCATCGGTTCGCAATCTTGGTTTGG‑CTTCGCC‑TTAATTAGGCCAGG‑TTTATTAACC‑TTAAt                                                                                                                                                >  1:523027/1‑147 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTcaca                                                                                                                                           <  1:905785/151‑1 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATATGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAAccg                                                                                                                                       >  1:406976/1‑151 (MQ=255)                 cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCATTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAAccgc                                                                                                                                      <  1:655690/151‑1 (MQ=255)                  aGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAAccgcc                                                                                                                                     >  2:903004/1‑151 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCATTCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACcgccgc                                                                                                                                   >  2:1087195/1‑150 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCAcc                                                                                                                                >  1:845035/1‑150 (MQ=255)                          cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGc                                                                                                                             >  1:1109521/1‑151 (MQ=255)                           aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGcc                                                                                                                            >  1:803740/1‑151 (MQ=255)                             cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGt                                                                                                                          >  2:1227959/1‑151 (MQ=255)                               ggTGATGAGAGGGGAGCAAGCTTATTCCATTTTAAATTATTC‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGT‑CTTCGCC‑ATTAATAAGCCAGT‑TTTTATCCCC‑TTTATTCACAACCGCCGCACCAGCCGTTc                                                                                                                        <  2:662224/151‑1 (MQ=255)                                     gagaCGGGAGCAAGTTTTTTCCTTTTAAAATTATTA‑CG‑GGTAAAATAAAGTTAAAATCCGGCACCCTTTCGCATTCCGGTTTTGG‑CTTCGCC‑ATAAATAAGCCAGT‑TTTTACCCCT‑TTTATTCAAAACCGCCGCACCAGCCGTTCCg                                                                                                                      <  2:791587/147‑1 (MQ=255)                                       gaggggAGAAGGCTT‑TTCCCATTTAAATTTTTACCC‑GGTAAAATCAAGTTAAAATGCGCCGCCCGTCCACATCCCTGGTTGGC‑TTTCCCC‑ATAAATAACCGAGT‑TTTTACCCCT‑TTTATCAACAACCGCCGCACCAGCCGTTCCGTACCAg                                                                                                                <  2:8442/148‑1 (MQ=255)                                           ggAGCAAGTTTTTTCCATTTAAAATTATTC‑CG‑GGTAAAAACAAGTTAAAATACGCCATCCGTTCGAAACCCGGTTTTGG‑CTTCGCC‑ATTAA‑AAGCCCGGTTTTTACCCCC‑TTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCa                                                                                                                 <  2:807979/146‑1 (MQ=255)                                            gAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCg                                                                                                            >  2:15318/1‑150 (MQ=255)                                               cAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCCGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        >  2:137332/1‑151 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  2:432400/1‑151 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                       >  1:957557/1‑151 (MQ=255)                                                 atgttttttccttttAAATTTTTC‑CCGGGTAAAATCAATTTAAAATCCGCCATCC‑TTCCAAATTCTTGTTTGGTCTTCCCC‑ATTAATAGCCCAGT‑TTTTACCCCT‑TTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                        <  1:1145646/139‑1 (MQ=255)                                                 aGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGc                                                                                                      >  2:699355/1‑151 (MQ=255)                                                    ttATCCCATTTTAAATTATTA‑CG‑GGTAAAATAAAGTTCAAATACGCCATCCGTTCGCAACCCTGTTTTGG‑CTCCGCC‑TTTAATAGGCCGGT‑TTTTATCCCC‑TTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                                    <  2:1020563/150‑1 (MQ=255)                                                         ccATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                               >  2:1213261/1‑150 (MQ=255)                                                            tttaAAATTATTC‑CC‑GGTAAAATAAAGTTAAAATACGGCATCCGTCCGAATTCCTGGTTTGG‑CTCCGCA‑ATTAAAAAGCCAGG‑TTTTACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGAACCAGGGCGATTCGCCAGCGTaa                                                                                              <  2:803743/148‑1 (MQ=255)                                                                aaaTTATTC‑CG‑GGTAAAATAAATTTAAAATACGGCATCCGTTCGCAACCCTGGTTTGG‑TTTCCCA‑ATTAAAAACCCAGT‑TTTTACCCCC‑TTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  2:1220809/151‑1 (MQ=255)                                                                aaaTTATTA‑CG‑GGTAAAACCAATTACAAATCCGGCTCCCGTCCGCAACCCTGGTTTGG‑CTCCGCC‑ATAAAAAGGCCAGT‑TTTTACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  1:6543/151‑1 (MQ=255)                                                                aaaTTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       >  2:416429/1‑151 (MQ=255)                                                                 aaTTATTC‑CG‑GGTGAAATAAAGTAAAAACCCGGCTTCTGTCCCCACCCCTGTTTTGT‑CTCCGCA‑ATTAATAGGCAAGT‑TTTTCCCCCT‑TTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                      <  1:1182363/151‑1 (MQ=255)                                                                  aatttttcccgggTAAAATAAATTTAAAATACGGCTTCCTTCCCCAATCCTGTTTTGT‑CTCCCCC‑ATAAATAA‑CCCGG‑TTTTTTCCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                       <  1:590574/140‑1 (MQ=255)                                                                       ttccc‑GGAGAAATAAATTTCAAATACGGCATCCGTCCGCAATCCTGGTTTGG‑CTCCGCC‑ATTAATAAGCGAGT‑TTTTATCACC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                    <  2:152603/145‑1 (MQ=255)                                                                       ta‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                                >  2:462774/1‑151 (MQ=255)                                                                             ggTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                              >  2:758033/1‑149 (MQ=255)                                                                                aaaacaaaGTTAAAATACGGCTCCCGTCCGCAATCCTGTTTTGG‑CTCCGCA‑TTTAATAAGCCGGT‑TTTTACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                           <  2:840817/143‑1 (MQ=255)                                                                                 aaaTCAAGTTAAAATCCGGCACCCGTTCGCAACCCGGTTTGGG‑CTTCGCC‑TTTAATAAGCCAGT‑TTTTATCCCC‑TTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                        <  1:1032364/151‑1 (MQ=255)                                                                                   aTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       >  2:456603/1‑150 (MQ=255)                                                                                   aTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                                       >  2:504302/1‑150 (MQ=255)                                                                                      aaGTTAAAATGCGGCATCTTTCCCCACCCCTGGTTTGG‑CTTCGCC‑ATAAAAAACCCGGT‑TTTTACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                      <  2:556298/148‑1 (MQ=255)                                                                                       aGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                    >  1:891079/1‑149 (MQ=255)                                                                                             aaTACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                             >  2:46586/1‑150 (MQ=255)                                                                                                       ccGTCCACAACCCTGGTTTGG‑CTCCGCT‑ATTAAAAAGCAAGG‑TTTTACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  <  2:633247/151‑1 (MQ=255)                                                                                                          tccGCATTCCTGGTTTGT‑CTCCGCT‑ATTAATAGGCCGGT‑TTTTACCCCT‑TTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                  <  1:1125389/146‑1 (MQ=255)                                                                                                               aaTCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                          >  1:1105061/1‑151 (MQ=255)                                                                                                                aTCTTGTTTTGG‑CTCCCCA‑ATAAAAAACCCTGT‑TTTTACCCCT‑TTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                               <  1:159622/145‑1 (MQ=255)                                                                                                                     ggTTTGG‑CTTCGCC‑ATTAATAAGCCAGT‑TTTTACCCCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                     <  2:933868/150‑1 (MQ=255)                                                                                                                      gTTTTG‑CTTCCCCAATAAAAAAGCCGGG‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  1:1057808/148‑1 (MQ=255)                                                                                                                       tttGG‑CTCCCCA‑ATAAAAAAGCAAGT‑TTTTATCCCC‑TTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataa                                        <  1:929810/145‑1 (MQ=255)                                                                                                                       tttGC‑CTTCCCC‑ATTAAAAAGCCAGT‑TTTTACCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                       <  2:1244028/146‑1 (MQ=255)                                                                                                                          gg‑CTTCCCA‑ATTAAAAACCCAGT‑TTTTATCCCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  2:1110938/149‑1 (MQ=255)                                                                                                                             cTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                             >  2:1109396/1‑151 (MQ=255)                                                                                                                              ttCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            >  2:322900/1‑151 (MQ=255)                                                                                                                                cccc‑ATAAAAAAGCGGGT‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                                 <  1:978688/142‑1 (MQ=255)                                                                                                                                cGCC‑AAAAATAGGCCGGT‑TTTTATCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTACTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                          <  1:755136/151‑1 (MQ=255)                                                                                                                                 cccaataaataaGCCAGT‑TTTTATCCCT‑TTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                            <  1:699352/142‑1 (MQ=255)                                                                                                                                 ccc‑AATAAAAAGCCCGGTTTTTACCCCT‑TTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt                          <  1:637797/150‑1 (MQ=255)                                                                                                                                  cc‑ATAAATAA‑CCAGGTTTTTATCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         <  1:112125/150‑1 (MQ=255)                                                                                                                                   caataaaaaaGCCGGT‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                         <  2:74367/142‑1 (MQ=255)                                                                                                                                         aaaaacccggtttttACCCCT‑TTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                           <  1:225454/132‑1 (MQ=255)                                                                                                                                          taaGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACGACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                  >  1:65966/1‑150 (MQ=255)                                                                                                                                           aaGCCGGT‑TTTTATCCCC‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa                 <  2:72486/150‑1 (MQ=255)                                                                                                                                                           aCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtttt  >  1:559801/1‑150 (MQ=255)                                                                                                                                                           |                                                                                                                                                       ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTA‑CC‑GGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCC‑ATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATTTT  >  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: ‑