Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	42.4%	101.4 / 16.5	61	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/25);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 4.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAT  >  NC_000913/4537961‑4538255                                                                                                                                                        |                                                                                                                                                      tCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCAcc                                                                                                                                                     >  1:180762/1‑151 (MQ=255)  ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCAcc                                                                                                                                                     <  2:714930/150‑1 (MQ=255)    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCttt                                                                                                                                                  >  2:851109/1‑151 (MQ=255)    ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCttt                                                                                                                                                  <  2:44658/151‑1 (MQ=255)     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCttt                                                                                                                                                  >  2:893301/1‑150 (MQ=255)     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCttt                                                                                                                                                  >  2:362056/1‑150 (MQ=255)     cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑a                                                                                                                                                <  2:812807/151‑1 (MQ=255)      gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑At                                                                                                                                               >  2:625105/1‑151 (MQ=255)          aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑ca‑c                                                                                                                                         <  2:947879/151‑1 (MQ=255)          aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑ca‑c                                                                                                                                         >  1:275120/1‑151 (MQ=255)          aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑ca                                                                                                                                           >  2:53620/1‑150 (MQ=255)           aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑ca‑c                                                                                                                                         >  1:504891/1‑150 (MQ=255)              gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAAcc                                                                                                                                     >  2:621921/1‑151 (MQ=255)                cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAAccgc                                                                                                                                   >  1:321301/1‑151 (MQ=255)                    tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACcgccgc                                                                                                                                >  1:466268/1‑150 (MQ=255)                       gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCAcc                                                                                                                             >  2:143708/1‑150 (MQ=255)                        ccATCTGTTGATGAGACGGGGACAAGCTTATTCCATTTTAAATTATTACCGTTGAAATCAATTTAAAATACGGCATCCCTCCGCAACCCTGGTTT‑GGCTTCGCC‑ATAAAAAACCGAGGTTTT‑TTCCCCTTT‑ATT‑CA‑CAACCGCCGCACCAg                                                                                                                           <  2:477777/151‑1 (MQ=255)                        ccATCTGGTGATGAGAGGGGAGCAAGCTTATTCCTTTTTAAATTATTCCCGGTGAAATAAAGTTAAAATACGACACCCTTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCCTTT‑ATT‑CA‑CAACCGCCGCACCa                                                                                                                            <  1:362051/150‑1 (MQ=255)                          aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGcc                                                                                                                         >  2:158994/1‑151 (MQ=255)                                    gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTCCCGGTAAAATCAAGTTCAAATACGGCATCCTTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCCCCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTAc                                                                                                                <  2:854903/150‑1 (MQ=255)                                         gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                          >  2:191658/1‑151 (MQ=255)                                         gggACCAATCTTATTCCTTTTAAATTTATTCCGGGTAAAACCAATTCAAAATACGCCATCCTTTCCAAATCCTGTTTT‑GGCTCCGCC‑ATAAAAAAGCTAGTTTTT‑ATCCCCTTT‑ATT‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCa                                                                                                              <  1:966957/147‑1 (MQ=255)                                                   ttATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                  >  2:307642/1‑149 (MQ=255)                                                    tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                               >  1:7573/1‑151 (MQ=255)                                                         cATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTaa                                                                                           >  1:798816/1‑150 (MQ=255)                                                          aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                         >  2:504111/1‑151 (MQ=255)                                                            tttAAATTATTCCCGGGGAAATAAAGTTCAAATACGGCATCCGTTCGCAATCCTGTTTT‑GGTTCCGCC‑ATTAATAAGCCAGTTTTT‑ACCCCCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc                                                                                        <  2:297044/150‑1 (MQ=255)                                                                aaTTATTGCCGGTGAAACCAGGTTAAAATACAAAATCCGTTCCCAATCTTAGTTG‑GATTTCCCC‑ATAAATAAGCCAGTTTTT‑ATCGCCTTT‑ATT‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                   <  2:436017/151‑1 (MQ=255)                                                                aaTTATTACCGGAAAAATAAAGTCAAAATACGCCCCCCTTTCCCAATCCGGTTTT‑GGCTCCGCC‑ATAAAAAAGCCAGTTTTT‑ATCCCCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                   <  1:533592/151‑1 (MQ=255)                                                                   tattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg                                                                                >  1:527895/1‑151 (MQ=255)                                                                          gggaaaaTAAAGTTCAAATACGGCTTCTGTTCGCAACCTTGTTTT‑GGCTTCCCC‑TTTAAAACGCCAGGTTTT‑ATCCCTTTT‑TTC‑CA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                         <  2:321303/147‑1 (MQ=255)                                                                           gtaaaaTCAAGTTAAAATCCGCCATCCGTTCGCAATCCTGTTTT‑GGCTTCCCC‑ATTAATAAGCCAGTTTTT‑ATCCCCTTT‑ATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg                                                                        <  2:486396/148‑1 (MQ=255)                                                                             gAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                       >  2:444939/1‑150 (MQ=255)                                                                             gAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                       >  2:574668/1‑150 (MQ=255)                                                                              aaataaatTTCAAATACGGCATCCGTTCGCAATCCTGTTTT‑GGCTTCGCC‑ATTAAAAAGCCAGGTTTT‑ATCCCCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                      <  2:749486/150‑1 (MQ=255)                                                                              aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGa                                                                      >  1:970864/1‑150 (MQ=255)                                                                              aaaTCAAGTCCAAATACGGCTCCCGTTCGCAATCCTGGTTT‑GGTTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                       <  2:87196/149‑1 (MQ=255)                                                                                  cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                  >  2:816869/1‑150 (MQ=255)                                                                                   aattcaaaacccGGCACCCGTCCACAATCTTGTTTT‑GGCTCCCCC‑ATTAA‑AAGCCGGGTTTTAACCCCCTTT‑TTC‑CC‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                 <  1:305302/139‑1 (MQ=255)                                                                                      ttCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                             >  1:969106/1‑151 (MQ=255)                                                                                         aaaTCCGGCTTCCGTTGGCATCCCTGGTTT‑GGTTTCGCC‑TTAAATAAGCCGGTTTTT‑ACCCCTTTT‑ATT‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                           <  2:480486/150‑1 (MQ=255)                                                                                         aaaTACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                          >  1:719455/1‑151 (MQ=255)                                                                                               ggCTTCCGCTCCCAATCCTGTTTT‑GGTTTCGCA‑ATTAATAAGCCAGGTTTT‑ATCCCCTTT‑ATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                    <  2:827657/151‑1 (MQ=255)                                                                                                    ccGTTCGCAATCCTGTTTT‑GGCTTCGCC‑ATTAATAAGCCAGTTTTT‑ATCCCCTTT‑ATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                   <  2:708921/147‑1 (MQ=255)                                                                                                       ttccccaaccTTGTTT‑TGCCTCCCCAATTAAAAACCCGGTTTTT‑ACCCCCTTT‑‑TT‑CACCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAt                                                <  1:584002/140‑1 (MQ=255)                                                                                                         cgaaaTCTTGTTTT‑GGCTTCGCA‑ATAAATAAGCCGGTTTTT‑ATCCCTTTT‑ATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                          <  2:180762/148‑1 (MQ=255)                                                                                                            cactcTGGTTTGGGCTCCGCA‑ATCAA‑AACCCCGGTTTT‑CCCCCCTTTAATT‑CA‑CAACCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                         <  1:353809/146‑1 (MQ=255)                                                                                                               ccTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                    >  2:404415/1‑151 (MQ=255)                                                                                                                    ttt‑GGCTTCCCT‑ATAAATAAGCCAGTTTTT‑ACCCCTTTT‑ATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                                <  2:884451/150‑1 (MQ=255)                                                                                                                     ttggccTTCCCC‑TTAAATAAGCCAGGTTTT‑ACCCCCTTT‑ATC‑CA‑AACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                  <  1:191299/143‑1 (MQ=255)                                                                                                                        ggCTTCACA‑ATTAAAAGGCATGGTTTT‑ACCTCTTTT‑CTT‑CA‑AACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                               <  1:24160/132‑1 (MQ=255)                                                                                                                         ggttCCCC‑AATAAAAAACCAAGTTTT‑TACCCCTTT‑TTTCCA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                              <  1:724013/147‑1 (MQ=255)                                                                                                                            tccccaataaataaGCCGGGTTTT‑ACCCCTTTT‑ATT‑CA‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                           <  2:935149/140‑1 (MQ=255)                                                                                                                             cccc‑ATTAATAACCCGGGTTTT‑ATCCCCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                        <  1:773500/148‑1 (MQ=255)                                                                                                                                cattaaaaaacccAGGTTTT‑ATCCCCTTT‑ATT‑CA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                     <  2:357471/140‑1 (MQ=255)                                                                                                                                   ttAATAAGCCGGTTTTT‑ATCCCTTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                   <  2:30110/150‑1 (MQ=255)                                                                                                                                     aaaaaGCCAGGTTTT‑TCCCCCTTT‑ATC‑CA‑AAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                 <  2:566932/147‑1 (MQ=255)                                                                                                                                      ataaGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt               >  2:169914/1‑151 (MQ=255)                                                                                                                                       taaGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt               >  1:887388/1‑150 (MQ=255)                                                                                                                                                   t‑ATCCCTTTT‑ATC‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGa   <  1:133354/150‑1 (MQ=255)                                                                                                                                                       cACCTTT‑TTT‑AC‑CACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt  <  1:754971/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                      TCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTT‑GGCTTCGCC‑ATTAATAAGCCAGGTTTT‑ATCACCTTT‑ATT‑CA‑CAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAT  >  NC_000913/4537961‑4538255

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