Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,511,923	0	A	C	46.8%	134.6 / 10.3	79	V141G (GTT→GGT)	yhfZ	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/7);  new base C (0/37);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 8.38e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCATTG  >  NC_000913/3511774‑3512073                                                                                                                                                             |                                                                                                                                                       aaCGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAac                                                                                                                                                       >  2:1324603/1‑151 (MQ=255)    ggTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAaccac                                                                                                                                                    >  1:95144/1‑151 (MQ=255)      ttATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCg                                                                                                                                                  >  1:1453840/1‑151 (MQ=255)       tATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCg                                                                                                                                                  >  2:251203/1‑150 (MQ=255)           cAGCCCCACGCGCTTCACGTTGGCGGACTCGCC‑TTTCCGGAAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTC‑CCGCCGC‑TAGTCGCGAAACCACCGCCATg                                                                                                                                             <  2:805562/151‑1 (MQ=255)                ccACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCAt                                                                                                                                         >  1:1256457/1‑150 (MQ=255)                  aCGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATac                                                                                                                                       >  1:579085/1‑150 (MQ=255)                    gcgcTTCGCGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCACCGTAGGTGTGCGCCCCCAC‑CCCCAGCGGGAGGCATAAGCC‑TTTTT‑GCGTGAGATACCTTC‑CCGCCGC‑TAGCCGCGAAACCACCGCCATGTCATACACg                                                                                                                                    <  2:61572/151‑1 (MQ=255)                    gcgcTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCACCGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATac                                                                                                                                       <  1:1099123/148‑1 (MQ=255)                        ttCACGTTTGCGGACTCGCC‑TTTACGGAAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGcc                                                                                                                                  <  1:1464157/149‑1 (MQ=255)                          cccGTTTGCGGACTCGCC‑TTTCCGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAC‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCAt                                                                                                                                <  1:68488/147‑1 (MQ=255)                           cactttGGGGACCCCCCTTTTACGGAAAATAAACT‑GGTGCCCCCAACCTTGGGTGTGCGCCCCAAC‑CCCCGGCGCG‑GGGAAAAGCCTTTTTT‑GGGGGAAATACCTTC‑CCCCCGT‑TAGTCGCAAACCCCCCGCCATGTCATACACGCCATTa                                                                                                                              <  2:1612483/148‑1 (MQ=255)                             gTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGc                                                                                                                           >  1:413486/1‑151 (MQ=255)                                  cGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGac                                                                                                                       >  2:1130511/1‑150 (MQ=255)                                     accccccttttCCGGCAAA‑CAACTGGGTGCCCCCCAAGTTGGGTGGGCGCCCCAAC‑CCCCGGCGCG‑GGCAAAAACC‑TTTTTGGCGGGAAATACCTTC‑CCCCCCC‑TGGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACt                                                                                                                    <  2:1523150/143‑1 (MQ=255)                                              ttACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGa                                                                                                            >  1:106397/1‑150 (MQ=255)                                              ttACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGa                                                                                                            >  1:628293/1‑150 (MQ=255)                                              ttACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGa                                                                                                            >  2:782814/1‑150 (MQ=255)                                               tACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAat                                                                                                          >  1:1492748/1‑151 (MQ=255)                                               tACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAat                                                                                                          >  2:1642304/1‑151 (MQ=255)                                                  ggaaaacaaaCG‑GGTGCCCCCCACCGTGGGTGTGCGCCCCCAC‑CCCCAGCGCGAGGCATAACCT‑TTTTT‑GCGTGAAATAACTTC‑CCCCCGC‑TAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAata                                                                                                         <  1:362908/141‑1 (MQ=255)                                                  ggCAAATCAACT‑GGTGCTCGCCACCGTGGGTGTGCGCCCCCAC‑CCCCAGCGCGAGGCAT‑AGCC‑TTTTTGGCGTGAAAAAACTTC‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAata                                                                                                         <  2:464091/149‑1 (MQ=255)                                                  ggCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATc                                                                                                       >  2:1536138/1‑151 (MQ=255)                                                      aaTCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCcgcg                                                                                                   >  1:771328/1‑151 (MQ=255)                                                         cAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGcc                                                                                                 >  2:1399706/1‑150 (MQ=255)                                                         cAACT‑GGTGCCCGCCAACGTAGGTGTGCGCCCCCAC‑CCCCAGCGCGGGGCAAAAGCT‑TTTTG‑GCGTGAGATAACTTC‑CCGCCGC‑TGGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGcc                                                                                                 <  1:945038/150‑1 (MQ=255)                                                          aaCT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCAc                                                                                               >  1:192687/1‑151 (MQ=255)                                                           accgggTGCTCGCCAACGTGGGGGTGCGCCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTG‑GCGTGAGATACCTTT‑CCGCCGC‑TAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCAc                                                                                               <  1:1200325/147‑1 (MQ=255)                                                               ggTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCa                                                                                            <  1:285094/150‑1 (MQ=255)                                                               ggTGCTCGCCAACGTAGGTGTGCGCCCCCAC‑CCCCAGCGCGGGGAATAAGCT‑TTTTT‑GCGTGAGATACCTTC‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAt                                                                                           <  1:686763/151‑1 (MQ=255)                                                                gTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtg                                                                                          >  1:1077452/1‑151 (MQ=255)                                                                   ccccccAACGTGGGTGGGCGCCCCCAC‑CCCCGGCGCGGGGCAAAAGCT‑TTTTG‑GCGGGAAAAAATTTC‑CCGCCAC‑TAGTCGGAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgt                                                                                         <  1:943079/145‑1 (MQ=255)                                                                     cgcaacgttgggTGTGCGCCCCCAC‑CCCCAGCGCGAGGCATACTCT‑TTTTT‑GCTGGAAATACCTTC‑CCGCCGC‑TAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgt                                                                                         <  2:1157916/137‑1 (MQ=255)                                                                     ccccACCGTAGGTGTGCGGCCCCAC‑CCCCAGCGCGGGGCATAACCT‑TTTTG‑GCGTGAGATAACTTC‑CCGCCGC‑TAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTgc                                                                                       <  1:1705701/147‑1 (MQ=255)                                                                       ccAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   >  1:757459/1‑151 (MQ=255)                                                                       caaaCGAAGGTGTGCGCCCCAAA‑CCCCAGCGCGGGGCAAAAGCT‑TTTTG‑GCGGGAGATAACTTC‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTgc                                                                                       <  1:859009/145‑1 (MQ=255)                                                                         aaCGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   >  1:1029244/1‑149 (MQ=255)                                                                          aCGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGt                                                                                >  1:1427376/1‑151 (MQ=255)                                                                               ggggggCGGCCCCAC‑CCCCGGCG‑GGGGGAAAAGCCTTTTTT‑GCGGGAGATACCTTC‑CCGCCGC‑TGGCCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaa                                                                            <  2:770753/145‑1 (MQ=255)                                                                               ggTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaa                                                                            >  1:1335377/1‑150 (MQ=255)                                                                               ggTGTGCGCCCCCAC‑CCCCGGGGCGGGCCAAAACCT‑TTTTG‑GGGTGA‑AAAACTTTCCCCCCGC‑TGGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaaa                                                                           <  2:1205877/151‑1 (MQ=255)                                                                                gtggGCGCCCCCAA‑CCCCAGCGCGGGGCAAAAGCT‑TTTTT‑GCGTGAGATACCTTC‑CCGCCGC‑TGGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaa                                                                              <  2:624218/147‑1 (MQ=255)                                                                                 tgtgCGCCCCCAA‑CCCCAGCGCGAGGCAAAAGCT‑TTTTT‑GCGGGAGATAACTTT‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAg                                                                          <  1:1328549/150‑1 (MQ=255)                                                                                  gggCGGCCCCAC‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAAATACCTTC‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGa                                                                        <  1:1572667/149‑1 (MQ=255)                                                                                  gcgcgcCCCCAC‑CCCCGGCGGGAGGCATAACCT‑TTTTT‑GGGTAAGAAACCTTC‑CCCCCGC‑TAGTCGGAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGa                                                                        <  2:1579958/149‑1 (MQ=255)                                                                                    gCGGCCCCAA‑CCCCAGCGCGAGGCATAAGCT‑TTTTT‑GCGTGAGATAACTTC‑CCGCCGC‑TAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGaa                                                                       <  1:448205/150‑1 (MQ=255)                                                                                      ggCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGc                                                                    >  1:1335904/1‑151 (MQ=255)                                                                                        ccccAC‑CCCCAGCGGGGGGCAAAACCT‑TTTTG‑GCGTGAGATACCTTC‑CCCCCGC‑TAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATg                                                                     <  2:1132958/148‑1 (MQ=255)                                                                                         cccAACCTCCAGCCGGAGGCAAATGCC‑TTTTT‑GCTTGAAATCTCTTT‑CC‑CCGCTTAGGCGCGAACCCCCGGCCATTTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATg                                                                     <  1:655649/148‑1 (MQ=255)                                                                                          ccAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATc                                                                >  2:1648518/1‑151 (MQ=255)                                                                                            aa‑CCCCAGCGCGAGGCAAAAGCT‑TTTTT‑GCGGGAAATAACTTT‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCa                                                               <  1:1310639/150‑1 (MQ=255)                                                                                             accTCCAGCGCGGGGCATAAGCT‑TTTTT‑GCGTGAGATAATTTC‑CCGCCGC‑TAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCaa                                                              <  1:747952/149‑1 (MQ=255)                                                                                             accTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATATCTTT‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCa                                                               <  1:1471105/148‑1 (MQ=255)                                                                                             accTCCAGCCGGAGGAATAAGCT‑TTTTT‑GCGTGAGAAAACTTC‑CCGCCGC‑TAGCGGCAAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCa                                                               <  2:683282/148‑1 (MQ=255)                                                                                               ccccAGCGCGGGGCAAAACCT‑TTTTG‑GCGGAAAATAACTTC‑CCCCCGC‑TGGTCGGAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCaaa                                                             <  2:1426055/148‑1 (MQ=255)                                                                                                    gcgcGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGc                                                       >  1:665413/1‑151 (MQ=255)                                                                                                     cgcgAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCt                                                      >  2:717297/1‑151 (MQ=255)                                                                                                        ggggCATAAGCT‑TTTTT‑GCGTGAGATAATTTC‑CCGCGCC‑TATTCGCAAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTg                                                          <  1:683943/142‑1 (MQ=255)                                                                                                        gAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTc                                                   >  2:342083/1‑151 (MQ=255)                                                                                                            acaaaacctttttg‑GCGGGAAAAAACTTC‑CCCCCGC‑AAGTCGCGAACCCCCCGCCATGTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTggg                                                        <  2:1719382/134‑1 (MQ=255)                                                                                                               aaccc‑TTTTT‑GGGTAAGATAACTTC‑CCGCCGT‑AAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCtt                                                     <  2:697137/139‑1 (MQ=255)                                                                                                               aaGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGc                                            >  1:984725/1‑151 (MQ=255)                                                                                                                acgctttttg‑GCGTGAGATACCTTC‑CCCCCGC‑TAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATcc                                              <  1:104410/144‑1 (MQ=255)                                                                                                                       ttt‑GCGGGAGATACCTTC‑CCGCCGC‑TAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCa                                      <  2:499137/150‑1 (MQ=255)                                                                                                                            cGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCt                                 >  2:946976/1‑151 (MQ=255)                                                                                                                              tGAGATACCTTC‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCtt                                <  1:457454/150‑1 (MQ=255)                                                                                                                                 aaatattttcccgccgC‑TAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCg                              <  2:265117/140‑1 (MQ=255)                                                                                                                                    accTTC‑CCGCCGC‑TAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTc                               <  1:598235/143‑1 (MQ=255)                                                                                                                                       ttt‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACg                      >  1:845238/1‑151 (MQ=255)                                                                                                                                        ttcccgccgC‑TAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACaa                         <  2:807065/145‑1 (MQ=255)                                                                                                                                         t‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg                    >  1:303413/1‑151 (MQ=255)                                                                                                                                           ccgccgC‑TGGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg                    <  2:413466/150‑1 (MQ=255)                                                                                                                                           ccgccgC‑TAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg                    <  1:1415289/150‑1 (MQ=255)                                                                                                                                              ccgc‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtata                <  1:884350/151‑1 (MQ=255)                                                                                                                                               cgc‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtat                 >  1:45125/1‑149 (MQ=255)                                                                                                                                                gc‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagg              >  1:1488959/1‑151 (MQ=255)                                                                                                                                                   tAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggc            >  2:193612/1‑151 (MQ=255)                                                                                                                                                    aGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggc            >  1:1653862/1‑150 (MQ=255)                                                                                                                                                       cgcaAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggcaagg        <  1:547759/151‑1 (MQ=255)                                                                                                                                                         cgAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggcaagggc      >  2:534861/1‑151 (MQ=255)                                                                                                                                                          gAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggcaagggc      >  2:594929/1‑150 (MQ=255)                                                                                                                                                            aaCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCAtt   >  1:541709/1‑151 (MQ=255)                                                                                                                                                            aaCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCAtt   >  2:174828/1‑151 (MQ=255)                                                                                                                                                             accaccGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCATTg  >  2:1578218/1‑151 (MQ=255)                                                                                                                                                             |                                                                                                                                                       AACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCC‑TTTACGGCAAATCAACT‑GGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCC‑TTTTT‑GCGTGAGATAACTTT‑CCGCCGC‑TAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCATTG  >  NC_000913/3511774‑3512073

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