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

ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253                                                                                                                                                      |                                                                                                                                                 aTCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCa                                                                                                                                                  <  2:1059227/150‑1 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                               >  2:193526/1‑151 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                               >  1:973308/1‑151 (MQ=255)   ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                               >  2:696368/1‑151 (MQ=255)      cGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                             >  2:879833/1‑150 (MQ=255)       gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACAGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                           >  1:88007/1‑151 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                          <  2:415982/151‑1 (MQ=255)        cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                          >  2:832239/1‑151 (MQ=255)         aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                         >  2:1072359/1‑151 (MQ=255)          aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                         >  2:364796/1‑150 (MQ=255)           aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTcac                                                                                                                                       >  2:1034349/1‑151 (MQ=255)                ccAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccg                                                                                                                                  <  1:216939/151‑1 (MQ=255)                 cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgc                                                                                                                                 <  2:263381/151‑1 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:317284/1‑150 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:133735/1‑150 (MQ=255)                     tACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                              >  2:381090/1‑150 (MQ=255)                        gCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                           >  1:867250/1‑150 (MQ=255)                         ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                           >  2:441158/1‑149 (MQ=255)                                       gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                             >  1:122544/1‑149 (MQ=255)                                       gaCGGGAGCAACTTTATTCCTTTTTAAATTATTCCCGGTAAAATCAATTTCAAATCCGGCATCCGTTCGCAACCCTGTTTTGGCTTCCCCATTAATAACCCAGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCa                                                                                                            <  2:110302/150‑1 (MQ=255)                                        aCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                          >  2:62619/1‑151 (MQ=255)                                        aCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                          >  1:326936/1‑151 (MQ=255)                                          gggAGCAAGTTTATTCCATTTTAAATTATTACCGGTAAAATCAATTTAAAATACGGCATCCTTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGTTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAc                                                                                                              <  2:900586/145‑1 (MQ=255)                                          gggAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                        >  2:207827/1‑151 (MQ=255)                                              gCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                    >  1:951125/1‑151 (MQ=255)                                                aactttATTCCTTTTAAACTTATTCCCGGGAAATTCAAGTTAAAATACGGAACCCGTTCGAAACCCGGGTTTGCCTGCCCCTTTAAAAAGCCAGTTTTTCACCCCTTTTTTCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCg                                                                                                  <  2:122459/147‑1 (MQ=255)                                                aaGCTTATTCCATTTTAATTTATTCCCGGAGAAATCAGGTTCAAATACGCCCTCGGTGCAAAATCCTGTTTTGGCTCCCCCATTAAAAAGCCTGTTTTTATCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                                   <  2:524449/150‑1 (MQ=255)                                                   cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  2:369437/1‑150 (MQ=255)                                                     tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                             >  2:955951/1‑151 (MQ=255)                                                         ccATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa                                                                                          >  2:1128464/1‑150 (MQ=255)                                                               tAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt                                                                                   >  2:719525/1‑151 (MQ=255)                                                                 aaTTATTACCGGTGAAATCAAGTTCAAATACGGCACCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                 >  2:751038/1‑151 (MQ=255)                                                                   ttattcCCGGTGAAACAAATTACAAATCCGGCTTCCGTTCGCACCCCGGGTTTGGCTCCGCAATAAAAAACCCGGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                                <  2:384808/150‑1 (MQ=255)                                                                       taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                           >  2:2838/1‑151 (MQ=255)                                                                        aCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGtt                                                                          >  1:1110720/1‑151 (MQ=255)                                                                         ccGGTAAAATCAAGTTAAAATACGGCATCCGTCCGAAATCCTGTTTTGGCTTCGCCATAAAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                           <  2:828333/149‑1 (MQ=255)                                                                           ggTAAAATCAATTTCAAATACGGCATCCTTTCACATTCCTGTTTTGGTTTCGCAATAAATGACCGAGTTTTTATCACTTTTATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg                                                                       <  2:698421/151‑1 (MQ=255)                                                                           ggTAAAATCAATTTAAAATACGCCATCCTTTCGCAACCCTGGTTTGGCTCCGCAATTAATAACCCGGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                         <  1:157562/149‑1 (MQ=255)                                                                             tGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                     >  1:749382/1‑151 (MQ=255)                                                                               aaataaatTTAAAATCGGCCACCCGTTCGAAACCCTGTTTTGGTTCCGCCATTAAAAAGCCGGGTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                   <  2:804401/151‑1 (MQ=255)                                                                               aaaTAAAGTTAAAATCCGCCATCCTTCCGCAATCCTGTTTTGGCTCCGCCATTAAAAAGCCAGTTTTTATCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                   <  2:788847/151‑1 (MQ=255)                                                                                    aattcaaaaTCCGGCACCCGTTCGAAATCCTGGTTGGCCTTCCCAATAAATAACCCAGTTTTTATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa                                                              <  1:713612/145‑1 (MQ=255)                                                                                    aattacAAATCCGGCACCCTTTCGAAATCCTGGTTTGCCTCCGCAATTAATAAGCAGGTTTTTACCCCCTTTTCTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                <  1:409129/144‑1 (MQ=255)                                                                                    aaGTTCAATTACGGCCTCCGTTCGAAATCCGGTTTTGGCTTCCCCATAAATAAGCCGGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                <  1:541473/149‑1 (MQ=255)                                                                                         cAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                          >  1:1009943/1‑150 (MQ=255)                                                                                          aaagccGGCACCCGTCCCCAACCCTGGTTGGGCTCCCCCATAAATAAGCCGGTTTTTATCCCCTTTTTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                             <  1:441993/141‑1 (MQ=255)                                                                                          aaaTAGGGCACCCTTCCCCAACCCTGTTTTGGCTCCGCCATAAAAACCCCAGTTTTTACCCCTTTTATCCCAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                           <  1:442538/148‑1 (MQ=255)                                                                                          aaaTACGGCATCCGTTCGCAATCCTGTTTTGGCTTCGCCATTAATAAGCCAGGTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                         <  2:877295/150‑1 (MQ=255)                                                                                            aTACGGCACCCTTTCGAAACCCGGTTTGGGTTTAGCCATTAAAAGGCCGGTTTTTATCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                      <  2:483171/151‑1 (MQ=255)                                                                                               cGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                    >  1:200719/1‑150 (MQ=255)                                                                                                        ttCGCATCCCTGGTTTGTTTCCCCCATTAATAAGCCGGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                           <  2:290981/150‑1 (MQ=255)                                                                                                             aaTCCTGTTTTGTCTTCGCCATTAATACGCCAGTTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaa                                          <  2:70694/146‑1 (MQ=255)                                                                                                                ctggttttGGCTTCCCCATTAATAAGCCGGTTTTTACCCCCTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                         <  2:555740/139‑1 (MQ=255)                                                                                                                 cTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                 >  1:335989/1‑151 (MQ=255)                                                                                                                    ttttGGCTCCCCCATAAAAAACCCTGTTTTTATCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                        <  1:333890/140‑1 (MQ=255)                                                                                                                        ggCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                          >  2:4506/1‑151 (MQ=255)                                                                                                                           ttCCCCATAAATAACCCGGTTTTTATCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                         <  2:479866/149‑1 (MQ=255)                                                                                                                            tcgcttttAAAACGCCGGTTTTTACCCCCTTCTTCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                          <  2:1051242/141‑1 (MQ=255)                                                                                                                               ccATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                   >  1:287641/1‑151 (MQ=255)                                                                                                                                 aTTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  >  1:757689/1‑150 (MQ=255)                                                                                                                                  ttAAAAACCCAGTTTTTATCCCCTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                   <  2:477539/148‑1 (MQ=255)                                                                                                                                    aataaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               >  2:1099631/1‑150 (MQ=255)                                                                                                                                    aataaCCGAGTTTTTACCCCCTTTATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGt                   <  1:279518/146‑1 (MQ=255)                                                                                                                                    aaaatcccggttttaATCCCGTTTATGCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  <  2:1107356/136‑1 (MQ=255)                                                                                                                                    aaaaacccggtttttACCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa              <  2:636902/140‑1 (MQ=255)                                                                                                                                    aaaaaGCCTGTTTTTCCCCCCTTTCTCCAAAACCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        <  1:1072224/138‑1 (MQ=255)                                                                                                                                    aaaaaGCAAGTTTTTATCCCTTTTATCCCAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                 <  2:973430/145‑1 (MQ=255)                                                                                                                                          ccggtttttCTCCCTTTTACTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa        <  1:1023200/146‑1 (MQ=255)                                                                                                                                          ccAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt          <  1:107376/149‑1 (MQ=255)                                                                                                                                          ccAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAg           <  1:1034220/148‑1 (MQ=255)                                                                                                                                          ccAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         >  2:525780/1‑150 (MQ=255)                                                                                                                                            aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc      >  1:415053/1‑151 (MQ=255)                                                                                                                                            aGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGc      >  1:235809/1‑151 (MQ=255)                                                                                                                                               ttttACCCCCTTTCTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa            <  1:345813/142‑1 (MQ=255)                                                                                                                                                tttATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  2:870951/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                 ATCCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537960‑4538253

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