Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,086	0	A	C	41.6%	151.3 / 14.2	78	Y580D (TAC→GAC)	ggt	gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (40/5);  new base C (0/32);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 9.45e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584942‑3585236                                                                                                                                                       |                                                                                                                                                       gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCccg                                                                                                                                                   <  1:527482/150‑1 (MQ=255)  aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCccgcc                                                                                                                                                 >  2:1103674/1‑151 (MQ=255)    aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGt                                                                                                                                               >  1:328541/1‑151 (MQ=255)     aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGtt                                                                                                                                              >  1:401332/1‑151 (MQ=255)      aGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGtt                                                                                                                                              >  2:112082/1‑150 (MQ=255)         cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAat                                                                                                                                          >  1:992352/1‑151 (MQ=255)         cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAat                                                                                                                                          >  1:613641/1‑151 (MQ=255)         cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAat                                                                                                                                          >  2:783196/1‑151 (MQ=255)         cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAat                                                                                                                                          >  2:79186/1‑151 (MQ=255)              gtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATcc                                                                                                                                     >  2:86530/1‑151 (MQ=255)                 aatttctttctCGGGGAACGCAAAACCCGCTGAAAAAATACCTACAGGATAACAAAAAAAAAAAGG‑CTCCCTTCGGCTTGCCC‑TAAAAAAATA‑CCCCTCTCCCCACGAAGAGGCCCC‑CTACCT‑TA‑AGTCCCCCGCCGTTAAATCATCCAcc                                                                                                                                  <  2:228239/148‑1 (MQ=255)                  aaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTTCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                 >  1:465415/1‑151 (MQ=255)                  aaTTCATTCTCGGGGAACGGAAAAACCGCTGAAAGAATATCTAAAGGATTAAAAAAAAAAAAAGG‑CTCCCTTCGGCTTGCCC‑TAAAAAAAGA‑CCCCTCTCCCCACGAAGAGGCCCG‑CTACCC‑TT‑AGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                 <  2:766100/151‑1 (MQ=255)                   attcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                 >  1:686097/1‑150 (MQ=255)                   attcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                 >  2:702389/1‑150 (MQ=255)                    ttcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                 >  1:890975/1‑149 (MQ=255)                        ttcTCGGGGAACGGAAAACCCGCTGGAAAAATATCTAAAGGATTACAAAAAAAAAAAGG‑CTCCCTTCGGCTTCCCC‑TGAAAAAATA‑CCCCTCTCCCCACGAAGAGGCCCG‑CTAACT‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAgcg                                                                                                                             <  2:653375/149‑1 (MQ=255)                        ttcTCGGGAAACGGAAAAGCCGCGGGAAAAATATCTAAAGGATAACAAAAAAAAAAAGG‑CTCCCTTCGGCTTCCCC‑TAAAAAAATA‑GCCCTCTCCCCACGAAGGGGGCCG‑CAACCT‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAg                                                                                                                               <  1:262648/147‑1 (MQ=255)                           gcaggaaaggacAAAGCCGCTGGAATATAACCTACAGGATTA‑GAAAAAAAAACGGTCTTCTTTCGGCTGCCCC‑TGAAAAAAGA‑CCCCTC‑TCCCCCGGAGGGGGCCC‑CCAACCTTT‑AACCCCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                            <  1:340709/136‑1 (MQ=255)                              gggAACGGCAAAACCGGGGAAAAAATATCTAAAGGATTACAAAAAAAAAAAGG‑CTCCTTCCGGTTTCCCC‑TGAAAAAATA‑CCCCTCTCCCCACGAAGGGGCCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                                         <  2:191429/147‑1 (MQ=255)                              gggAACGGAAAACCCGCTGGAAGAATATCAACAGGATTACAAAAAAAAAAAGG‑CTCCCTTCGGCTTGCCC‑TGACAAAATC‑CCCCTCTCCCCACGAGGAGGCCCG‑CTACCC‑TT‑ATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                      <  1:898278/150‑1 (MQ=255)                               ggAAGGGAAAACCCGCGGGAAAAATTTCTACGGGATTACAAAAAAAAAAAGG‑CTCCCTCGGGTTTGCCC‑GAAAAAAATA‑GCCCTCTTCCCACAAAGAGGCCCG‑CTAACT‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGa                                                                                                                    <  2:363492/151‑1 (MQ=255)                                gAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTCACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     >  1:640423/1‑149 (MQ=255)                                   cGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                                >  1:785136/1‑151 (MQ=255)                                        aaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  1:230415/1‑150 (MQ=255)                                            cGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                       >  1:476611/1‑151 (MQ=255)                                              cTGAAAGAATATCTACAGAATTACGAAAAAAAAAAGG‑CTCCCTTCGCCTTGCCC‑TAACAAAATC‑CCCCTCTCCCCACGAAGAGGGCCG‑TTACCC‑TT‑AGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                      <  1:529126/150‑1 (MQ=255)                                               tgaaaaaatcccaaCAGGATTACAAAAAAAAAAGGG‑CCCCTTCCGCCGTGCCC‑GGACAAACTC‑CCCTTCCCCCAAAGAAGAGGGCCG‑TTGTCG‑TG‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                      <  1:675736/136‑1 (MQ=255)                                               tGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     >  1:134729/1‑150 (MQ=255)                                                ggAAAAATATCTACAGGATTACAAAATAAAAAAGC‑TTCCCTTCGGTTTCCCC‑TAAAAAAATA‑CCCCTCTCCCCACGAAGGGGGCCC‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                      <  1:1092376/148‑1 (MQ=255)                                                 gaagaaTATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  >  2:861869/1‑151 (MQ=255)                                                 gaagaaTATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  >  1:996300/1‑151 (MQ=255)                                                  aaaaaTATCTACAGGATAACAAAAAAAAAAGGG‑CTCCCTGCGGTTGCCCC‑TAAAAAAAAA‑CCCCTCTCCCCACGAAGAGGCCCG‑CTACCC‑TA‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  2:808223/147‑1 (MQ=255)                                                     aaTTTCTACAGGATTAAAAAAAAAAAAAGG‑CTCCCTCCGGCTTGCCC‑TGAAAAAATA‑CCCCTCTCCCCACGAAGAGGCCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              <  2:1099061/151‑1 (MQ=255)                                                         ttaaCAGAATGACAAAAAAAAAAAGC‑CTCCCTTCGGCTCCCCC‑AGAAAAAATA‑GCCCTCTCCCCACGAATGGGGCCG‑CAAACC‑TG‑AGTTCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  2:786749/145‑1 (MQ=255)                                                             acggaataacaaaaaaaaAGGG‑CTCCCTTCGGCTT‑CCC‑TGAAAAAAAACCCCCTCTCCCCCCGAAGGGGCCCG‑CTACCT‑TT‑ACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                        <  1:159472/139‑1 (MQ=255)                                                              aGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  2:979035/1‑151 (MQ=255)                                                               gaaTTAAGAAAAAAAAAAGC‑CTACCTTGGCTTTCCCC‑AAAAAAAATC‑CCCCTCTCCCCACGAAG‑GGGCCGCCAAACT‑TT‑ACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     <  1:80501/148‑1 (MQ=255)                                                                 ataacaaaaaaaaaaaGG‑CTCCCTGCGGCTTCCCC‑TGAAAAAAGA‑CCCCTCTCCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                                    <  2:163383/143‑1 (MQ=255)                                                                 aTTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                                   >  2:392277/1‑150 (MQ=255)                                                                 aTTACAAAAAAAAAAAGC‑CTCCCTTCGGCTTGCCC‑TGAAAAAATA‑CCCCTCTCCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                                    <  2:796644/149‑1 (MQ=255)                                                                     aaaaaaaaaaaaAG‑CCTCCTTGCGCTTCCCCTTGAAAAAATA‑GCCCTCTCCCCACGAAGGGGCCCC‑CTAACC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                <  2:652956/148‑1 (MQ=255)                                                                      gAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             >  2:877474/1‑151 (MQ=255)                                                                      aaaaaaaaaaaGG‑CTCCCTTCGCCTTCCCC‑GAAAAAAAAA‑GCCCTCTCCCAACGAGGGGGGCCG‑CTACCC‑TT‑ATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                               <  2:891850/148‑1 (MQ=255)                                                                       aaaaaaaaaaGG‑CTCCCTTCGGCTTCCCC‑TGAAAAAATA‑CCCCTCTCCCCACGAAGGGGCCCG‑CAAACC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            <  2:763437/151‑1 (MQ=255)                                                                        aaTAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                           >  1:594456/1‑151 (MQ=255)                                                                         aTAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCGGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                          <  2:785714/151‑1 (MQ=255)                                                                           aaaaaaGG‑CTCCCTTCGGCTTGCCC‑TGAAAAAATA‑CCCCTCTCCCCAGGAAGAGGGCCG‑CTACCT‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                        <  1:839109/151‑1 (MQ=255)                                                                           aaaaaaGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                        >  2:481289/1‑151 (MQ=255)                                                                            aaaaaGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                       >  2:202625/1‑151 (MQ=255)                                                                                 gg‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  >  1:38339/1‑151 (MQ=255)                                                                                    ctccttccgccTTCCCC‑TAAAAAAATA‑CCCCTCTCCCCACAAAGAGGCCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                                 <  2:717676/140‑1 (MQ=255)                                                                                    cTACCTTCGGCTAGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  1:436446/1‑151 (MQ=255)                                                                                     tACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  1:590501/1‑150 (MQ=255)                                                                                      aCCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                              >  1:349108/1‑151 (MQ=255)                                                                                      aCCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                              >  1:586718/1‑151 (MQ=255)                                                                                          tcggttgcccc‑TAAAAAAATA‑CCCCTCTCCCCCAGAAGGGGCCCG‑CAACCC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           <  2:324070/142‑1 (MQ=255)                                                                                          tCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGcc                                                          >  1:950370/1‑151 (MQ=255)                                                                                             gCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                       >  2:189539/1‑151 (MQ=255)                                                                                               ttcccc‑TAAAAAAATA‑CCCCTCTCCCCACAAGGAGGCCCG‑CAACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                       <  2:335918/146‑1 (MQ=255)                                                                                                  ccc‑TAAAAAAATC‑CCCCTCTCCCCACGAGGAGGGCCC‑CAACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  2:53802/148‑1 (MQ=255)                                                                                                          aaaaTC‑CCCCTCTCCCCACGAAGGGGCCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                           <  1:177523/151‑1 (MQ=255)                                                                                                           aaaacgcccccCTCCCCACAAAG‑GGGCCC‑CCAACC‑TTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCAGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCAc                                            <  2:213191/143‑1 (MQ=255)                                                                                                            aaTA‑GCCCTCTCCCCAGGAGGAGCGCCG‑CTACCC‑TG‑ATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  1:533961/151‑1 (MQ=255)                                                                                                              tA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       >  2:1001201/1‑151 (MQ=255)                                                                                                                 gCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg                                      >  1:259940/1‑150 (MQ=255)                                                                                                                  cccTCTCCCCACGGAGGGCGCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  2:526418/148‑1 (MQ=255)                                                                                                                  cccTCTCCCCACGAAGAGGCCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  2:827449/148‑1 (MQ=255)                                                                                                                        tCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg                              >  2:1078275/1‑151 (MQ=255)                                                                                                                         cccAGGAGGAGGCCCG‑CTACCC‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt                                  <  2:435934/146‑1 (MQ=255)                                                                                                                           cACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                           >  2:425181/1‑151 (MQ=255)                                                                                                                                agagGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                      >  1:390688/1‑151 (MQ=255)                                                                                                                                  ggggcccccaaccc‑TT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCACCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagca                     <  1:877344/139‑1 (MQ=255)                                                                                                                                    ggCCC‑CTAACC‑TT‑ATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                   <  2:434643/150‑1 (MQ=255)                                                                                                                                       cG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                >  1:775598/1‑150 (MQ=255)                                                                                                                                                ccttgctACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGCCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg            <  2:555737/142‑1 (MQ=255)                                                                                                                                                  t‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCg      >  2:519734/1‑151 (MQ=255)                                                                                                                                                    aGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       <  1:113579/149‑1 (MQ=255)                                                                                                                                                       aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  2:218879/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       GATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATA‑GCCCTCTTCCCACGAAGAGGGCCG‑CTAACC‑TT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584942‑3585236

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