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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253                                                                                                                                                        |                                                                                                                                                 aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCAc                                                                                                                                                 <  1:977774/151‑1 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCAcc                                                                                                                                                <  2:388220/150‑1 (MQ=255)     ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCttt                                                                                                                                             <  1:483354/151‑1 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTa                                                                                                                                            <  2:857921/151‑1 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTAt                                                                                                                                           >  2:756669/1‑151 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTAtt                                                                                                                                          <  1:943745/151‑1 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTc                                                                                                                                         >  1:397529/1‑150 (MQ=255)            aCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTcaca                                                                                                                                      <  1:692176/151‑1 (MQ=255)             ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTcac                                                                                                                                       <  1:377552/149‑1 (MQ=255)               gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAAcc                                                                                                                                   <  1:889031/151‑1 (MQ=255)                      aCGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCAc                                                                                                                            >  2:743874/1‑151 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                           >  2:576353/1‑150 (MQ=255)                         ccactgggTGAGGAGGGGGAAGAAGGCTTTTTCAATTTTAATTTATTCCGGGAGAAACCAAGTTCAAATCCGGCATCCGTTCGCAATCCTGGTTGGGCTTCGCCATA‑AATAAGCCAGT‑TTTTCTCCCTTTTATCCACAACCGCCGCAcc                                                                                                                           <  1:232945/143‑1 (MQ=255)                                       gaCGGGAGCAGGTTTATTCCTTTTTAATTTTTTACCGTTGAAATCAAGTTAAAATACGGCACCCTTTCGAACTCTTGTTTGGGCTTCGCAATT‑AAAAAGCCAGG‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                             <  2:78760/149‑1 (MQ=255)                                          gggAGAAAGCTTATTCCATTTAAAATTATTCCCGGAAAAATCAATTTAAAATCCGCCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATA‑AAAAAGCAAGG‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                          <  2:202989/149‑1 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTCCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAACCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                   <  2:289113/150‑1 (MQ=255)                                                aaGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                  >  2:988518/1‑151 (MQ=255)                                                   cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  1:519455/1‑150 (MQ=255)                                                   cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  2:23895/1‑150 (MQ=255)                                                    ttATTCCTTTTTAAATTTTTCCCGGTAAAACCAATTTCAAATACGGCCCCCGTCCGCAATCCTGTTTGGGCTCCCCCATT‑AATAAGCCAGT‑TTTTATCACTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                               <  1:11099/150‑1 (MQ=255)                                                    ttATTCCATTTAAAATTATTCCCGGTAAAATAAATTTAAAATACGGCATCCGTTCGCAACCCGGTTTTGGCTTCGCCATA‑AATAAGCCAGG‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                               <  2:716495/150‑1 (MQ=255)                                                       ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                           >  2:332087/1‑151 (MQ=255)                                                       ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                           >  2:523728/1‑151 (MQ=255)                                                       ttCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑CTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTTGCCAGCg                                                                                            >  1:955857/1‑150 (MQ=255)                                                       ttCCATTTTAAATAATCACCGGAGAAAAAAGGTCCAAATCCGGCACCCGTCCGAAATCCTGGTTGGGCTCCTCAATA‑AAAACGCCAGC‑TTTAACCCCCTTATTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                              <  1:816784/148‑1 (MQ=255)                                                          cttttaAAATTATTCCCGATAAATTAAAGTTAAAATACGCCGCCCGTCCGCAATCCGGTTTTGGCTTCGCCATA‑AATAAGCCAGG‑TTTTATCCCCTTTATCAAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                          <  2:194665/147‑1 (MQ=255)                                                            tttaAAATTTTTCCGGGTAAAATCAAGTTAAAATACGCCTCCCTTTCCCAACCCTGGTTGGGCTCCCCCATA‑AAAAAGCCAGT‑TTTTACCCCCTTTTTCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                             <  1:23480/144‑1 (MQ=255)                                                                   ttattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                               >  2:246421/1‑151 (MQ=255)                                                                       taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                           >  2:692549/1‑151 (MQ=255)                                                                        aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                          >  2:26066/1‑151 (MQ=255)                                                                         cgggTAAAATAAATTTAAAATCCGCCCCCTGTCCGCATCCCTGTTTTGTCTTC‑CCATTAAAAAACCCGGG‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtgtg                                                                            <  1:244546/146‑1 (MQ=255)                                                                         ccGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACACCATTGCTGTGTTg                                                                         >  2:835549/1‑151 (MQ=255)                                                                         ccGGTGAAATCAAGTTAAAATCCGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGT‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                          <  2:310438/150‑1 (MQ=255)                                                                                 aTCAAGTTCAAATCCGGCATCCGTTCGCAATCCTGGTTTGGCTTCCCCATT‑AATAAGCCAGG‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 <  2:292973/151‑1 (MQ=255)                                                                                  tCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 >  2:95895/1‑150 (MQ=255)                                                                                    aatttaaaaTAGGGCTTCCGTTCGAAACCCTGGTTTGGTTCCCCAATT‑AATAACCAAGT‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 <  2:872690/142‑1 (MQ=255)                                                                                      gTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                               >  1:116997/1‑148 (MQ=255)                                                                                       ttaaaaTCCGGCACCCGTTCGCAATCCTGGTTTGCTTTCGCCATT‑AAAAAGCCAGG‑TTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                            <  2:1024402/147‑1 (MQ=255)                                                                                               cGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    <  2:998460/150‑1 (MQ=255)                                                                                                ggCACCCGTTCGCATTCCGGTTTTGGCTTCCCCATT‑AATAAGCCAGT‑TTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                  <  2:807302/151‑1 (MQ=255)                                                                                                ggCACCCGTTCGCAACCCTGTTTTGGCTTCGCCATT‑AATAACCCAGG‑TTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAAc                                                  <  2:206949/151‑1 (MQ=255)                                                                                                   aTCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                >  2:1002660/1‑150 (MQ=255)                                                                                                     cctttCGCAACCCTGGTTTGGCTCCGCTTTA‑AAAACCCCAGT‑TTTTCTCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                <  1:202794/145‑1 (MQ=255)                                                                                                       gTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTg                                            >  2:57063/1‑150 (MQ=255)                                                                                                         tCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                           >  2:373390/1‑149 (MQ=255)                                                                                                          ccaatcccTGGTTTGTCTTCGCAATA‑AAAAACCCGGT‑TTTTACCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                         <  1:272853/144‑1 (MQ=255)                                                                                                          cGCAACCGTGGTTTGGTTGCGCTATT‑AATAAGCCAGT‑TTTTACCCCTTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                         <  1:858526/150‑1 (MQ=255)                                                                                                             aacccTGTTTTGGCTTCGCCATA‑AAAAACCCGGG‑TTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                          <  2:833450/143‑1 (MQ=255)                                                                                                             aaTCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                     >  2:958591/1‑151 (MQ=255)                                                                                                                     tttGGCTTC‑CCATTAAATAACCCGGT‑TTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                <  2:626832/148‑1 (MQ=255)                                                                                                                     tttGGCTCCCCCATT‑AAAAAGCCAGT‑TTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                               <  2:528673/149‑1 (MQ=255)                                                                                                                     tttGGC‑TCCCCCTT‑AAAAAACCAGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                              <  1:20215/150‑1 (MQ=255)                                                                                                                     tttGCCTTCCCCTTA‑AAAAAGCCAGT‑TTTTACCCCCTTTATTCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                      <  1:637557/142‑1 (MQ=255)                                                                                                                     tttGCCTCCGCCATA‑AAAAACCCAGT‑TTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                 <  1:86904/147‑1 (MQ=255)                                                                                                                     tttGCCTCCGCAATA‑AATAACCCAGG‑TTTTACCCCCTTTTTTCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                <  1:162765/148‑1 (MQ=255)                                                                                                                        ggCTTCGCCATT‑AAAAAGCGAGT‑TTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                          <  1:650/151‑1 (MQ=255)                                                                                                                           ttCCCCATT‑AATAAGCCGGG‑TTTTTTCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                         <  2:973653/149‑1 (MQ=255)                                                                                                                             ccccATA‑AAAACGCGAGT‑TTTTATCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        <  1:54062/146‑1 (MQ=255)                                                                                                                              gCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                    >  2:645309/1‑151 (MQ=255)                                                                                                                                ctgaaaatagGCCGGG‑TTTTATCACGTATATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAAACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                    <  2:510857/145‑1 (MQ=255)                                                                                                                                 aTT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  >  2:846086/1‑150 (MQ=255)                                                                                                                                     aataaGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  1:101131/1‑150 (MQ=255)                                                                                                                                         aacccggtttttATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg           <  1:1033963/143‑1 (MQ=255)                                                                                                                                           ccAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         >  1:487499/1‑150 (MQ=255)                                                                                                                                            cAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         >  1:694712/1‑149 (MQ=255)                                                                                                                                               gtttttATCCCGTTAATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAg       <  2:669215/147‑1 (MQ=255)                                                                                                                                                 ttttATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc      <  2:715055/148‑1 (MQ=255)                                                                                                                                                   ttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  1:415884/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                 ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATT‑AATAAGCCAGG‑TTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253

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