Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,382	0	A	G	45.1%	85.4 / 25.2	71	intergenic (+84/‑63)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base A (7/32);  new base G (32/0);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 9.94e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCG  >  NC_000913/15239‑15529                                                                                                                                                |                                                                                                                                                          ccGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTcc                                                                                                                                                   <  1:160281/151‑1 (MQ=255)     gcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGA‑ATccc                                                                                                                                                  >  1:39964/1‑148 (MQ=255)     gcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTccc                                                                                                                                                  >  2:1012769/1‑148 (MQ=255)     gcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATa                                                                                                                                               <  2:445231/151‑1 (MQ=255)        cAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGGAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCG‑CTccc                                                                                                                                                  >  1:520888/1‑145 (MQ=255)              gCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CCCCCCCAAACG‑CT‑aa                                                                                                                                     >  1:285675/1‑150 (MQ=255)                    ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTaa                                                                                                                                <  2:592181/149‑1 (MQ=255)                    ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTaa                                                                                                                                <  2:502183/149‑1 (MQ=255)                    ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTa                                                                                                                                 <  2:77387/148‑1 (MQ=255)                    ttGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAg                                                                                                                               <  2:118766/150‑1 (MQ=255)                     tGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCACAAACG‑CT‑AACTTAAagg                                                                                                                             >  2:254386/1‑151 (MQ=255)                      gATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAggg                                                                                                                             <  2:856435/150‑1 (MQ=255)                         ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑CGCCCCTCAAAG‑CT‑AACTTTAgg                                                                                                                              >  2:668866/1‑146 (MQ=255)                            gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTgg                                                                                                                       <  1:66979/150‑1 (MQ=255)                            gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGt                                                                                                                      <  2:133978/151‑1 (MQ=255)                             tgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGA‑CTCCCAAAAGCG‑GC‑AAATTTAGGGGTGTgg                                                                                                                       >  1:14511/1‑149 (MQ=255)                                 gaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCCCAAAAG‑CT‑AAATTTAgggg                                                                                                                            >  2:1013187/1‑139 (MQ=255)                                  aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACg                                                                                                                <  2:39966/151‑1 (MQ=255)                                  aagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACg                                                                                                                <  1:182841/151‑1 (MQ=255)                                   agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑c                                                                                                              <  1:570257/151‑1 (MQ=255)                                       ttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACA‑CACCCCAAAGAA‑CT‑CACTTATGGGGTGTGGGATTatg                                                                                                                >  2:279307/1‑144 (MQ=255)                                        tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCCAAGCG‑CC‑CACTTAAAGGGTTTGGGGTTATG‑CCTgc                                                                                                          >  1:922967/1‑149 (MQ=255)                                        tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGa                                                                                                          <  2:92419/150‑1 (MQ=255)                                          tGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CCCCCCTAAGCG‑CT‑CAATTTAGGGGTggggggttttgcctggt                                                                                                          >  1:562993/1‑136 (MQ=255)                                              gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCATAAAAG‑CC‑AACTTTAGGGGTGTGGGGTTATG‑CCTGCTATg                                                                                                      >  1:680411/1‑148 (MQ=255)                                              gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCCTTAAAG‑CC‑AACTTTAGGGGTGTGGGATTAAG‑CCCGATATGG‑tt                                                                                                  >  2:799985/1‑151 (MQ=255)                                                ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CT‑AACTTTAAGGGTGTGGGATTACG‑CCTGATATTA‑ttt                                                                                                 >  2:70093/1‑150 (MQ=255)                                                  tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑ATCCCCCAAACG‑CC‑AAATTTAGGGGTGTGGGATTTCG‑CCCGATATTA‑ttt                                                                                                 >  2:50998/1‑148 (MQ=255)                                                  tGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCCTAAAAG‑GT‑AAATTTAAGGGTGTGGGATTTTG‑GCTGATATGA‑ttt                                                                                                 >  2:85461/1‑148 (MQ=255)                                                   gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACg                                                                                             <  1:90447/151‑1 (MQ=255)                                                    aCCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGt                                                                                            <  1:161803/151‑1 (MQ=255)                                                     cccGCTAACCTTCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAACG‑CT‑AACTTAAAGGGTGTGGGGTTTCG‑CCTGATATTA‑ttt                                                                                                 >  2:659679/1‑145 (MQ=255)                                                      ccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CCCCCCTAAACG‑CT‑AACTTAAAGGGTGTGGGATTTAG‑CCCGATATGA‑TTTAAAGt                                                                                            >  1:48197/1‑149 (MQ=255)                                                      ccGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGc                                                                                          <  2:351006/151‑1 (MQ=255)                                                       cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCTTAAGGC‑CA‑AATTTAAGGG‑TGGGGGATTACG‑cccggctttgttt                                                                                                  >  1:182816/1‑132 (MQ=255)                                                        gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑CTCCCCTAAACG‑CC‑AACCTAAAGGGTGTGGGTTTATG‑CCCCGTATGA‑ttt                                                                                                 >  1:970161/1‑142 (MQ=255)                                                         cTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CTCCCCTAAACG‑CC‑AACTTTAGGGGTGTGGGATTACG‑CCTGATATGA‑TTTAAAGt                                                                                            >  2:237602/1‑146 (MQ=255)                                                             ccTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGA‑CCCCCCTAAGCG‑CC‑AACTTAAAGGGTGTGGGATTAAG‑CCTGATATGA‑ttttaagggcc                                                                                         >  1:467820/1‑137 (MQ=255)                                                                    aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CC‑AACTTTAAGGGTGTGGGATTATG‑CCTGGTATTA‑TTTAAAGTGCCCATGAAtt                                                                                 >  2:413812/1‑146 (MQ=255)                                                                      aaGCCTGTCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑AACCCCTTAACGCCC‑AACTAAAAGG‑TGTGGTATTTCG‑CCTGATATTT‑TTTaatcgggcggtg                                                                                     >  1:682245/1‑130 (MQ=255)                                                                        gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CT‑AACTTAAGGGGTGTGGGATTACG‑CCTGATATGA‑TTTTAAATGCCGATGAAtt                                                                                 >  2:894173/1‑142 (MQ=255)                                                                            gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACA‑CTCCCATAAGCG‑CT‑CACTTAAAGGTTGTGGGATTTTG‑CCCGATATTA‑TTTTAAGTGCCGATGAATTTCTCTCCCGATa                                                                     >  1:273020/1‑150 (MQ=255)                                                                              ccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCATCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATa                                                                     <  1:813608/148‑1 (MQ=255)                                                                               cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                 <  1:80035/151‑1 (MQ=255)                                                                               cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                 <  2:285737/151‑1 (MQ=255)                                                                                gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTgg                                                                <  2:221084/151‑1 (MQ=255)                                                                                gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGCA‑CTCCCATAAGCG‑CT‑AACCTAAGGGGTGTGGGGTTAAG‑CCTGATATTA‑TTTTACGTGCCGCTGAAATTTTCTCCCCATAAATgg                                                                >  1:773469/1‑151 (MQ=255)                                                                                     cAGGCCTGGGTAAAAATAGGGGGCGTTGAAGATATGCGTGCAGCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAg                                                            <  2:889876/150‑1 (MQ=255)                                                                                        gCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATTA‑CTCCCCTAAGAG‑GT‑AACTTTAAGGTTGTGGTTTTACGACCTG‑TATTA‑GTTAATGTGCCGCCGAATTAATCCCATGCTAACTGGGCAgga                                                          >  1:894026/1‑147 (MQ=255)                                                                                         ccTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAACG‑CT‑AACTTTAGGGGTGTGGGATTTCG‑CCTGGTATTA‑TTTTAAGTGCCGGTGAAATACTCTCACGATAAATGGGCAGCaaat                                                       >  2:991667/1‑149 (MQ=255)                                                                                           tGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑ATCCCCTAAAAG‑GTAAAATTAAGGG‑TGTGGGGTTTCG‑CCTGATATGG‑TTTTAAGTGCCGGTGAATTACTCTCTCGATAACTGGTCAGCAAtt                                                       >  1:951817/1‑149 (MQ=255)                                                                                            gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt                                                     <  1:498502/150‑1 (MQ=255)                                                                                              gTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAGGA‑CCCCCACAAGAG‑CT‑TACTTAAGAGGTGTGGGATTAGG‑CCTGTTATGATTTTAAAGT‑CGGATGAATTATTCTCTCCATTACttagtagccattcttg                                                   >  1:648249/1‑135 (MQ=255)                                                                                               tAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCA‑CTCCCCTAAGCG‑CC‑AACTTAAGGGGTGTGGGATTACG‑CCTGGTATGA‑TTTTACGTGGCGATGAAATACTCTCACGATAACTGGGCAGCAATTTTTGcc                                                 >  2:73462/1‑151 (MQ=255)                                                                                                   aaTAGGGGGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCA‑ATCCCCTAAAAG‑CG‑AAATTTAGGGGTTGGGGATTACC‑GCCGATATGA‑TTTTAAGTGCCCATGAATTAATCTCACCGAAAATTGTTTGCAATTCTTGccc                                                >  2:272594/1‑148 (MQ=255)                                                                                                     tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTA‑CTCCCATAAACG‑CT‑AAATTTAGGGGTGGGGGATTAAG‑CCCGATATTA‑TTTTAAGTGCCGATGAAATTATTTCCCGATAACTTGGTAGCAATTTTTGCCCCTA‑tt                                          >  1:952474/1‑151 (MQ=255)                                                                                                       gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTgg                                        <  2:150152/151‑1 (MQ=255)                                                                                                         gTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTa                                      <  2:201183/151‑1 (MQ=255)                                                                                                                gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCg                                <  1:728554/150‑1 (MQ=255)                                                                                                                  aGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCgaag                             <  2:567626/151‑1 (MQ=255)                                                                                                                  aGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCgaag                             <  1:854663/151‑1 (MQ=255)                                                                                                                    atatGCGAGCACCTGTAAAGTGGCGGGGGGGA‑CTCCCCTAAGCG‑CT‑TAATCAATGGGTTTGGTATTTAG‑CCCTATATGA‑TTTTACGTGCCGCTGAATTATTTTTACGATAACTTGTTAGCAATTTCTGCCCATATTTGGTAAcccca                                >  1:128497/1‑146 (MQ=255)                                                                                                                    atatGCGAGCACCTGTAAAGTGGCGGGGGGCA‑CCCCCCTAAACG‑CT‑AAATTAAAGGGTGTGGGATTTCG‑CCCGATATTA‑TTTTACGTGCCGATGAAATTCTCTCACCCTAACTGGTCAGCAATTCTTGCCCCTA‑TTTGTAAAccc                                 >  2:936608/1‑145 (MQ=255)                                                                                                                       tGCGAGCACCTGTAAAGTGGCGGGGGGCA‑ATCCCCTCAACG‑CC‑AAATTTAGGGGTGTGGTATTACG‑CCTGGTATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACCGGGCCGCAATTCTGGCCCCTA‑TTGGTAAGCCCcaagaa                           >  1:896565/1‑148 (MQ=14)                                                                                                                        gCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTg                        <  1:409357/150‑1 (MQ=255)                                                                                                                               ccTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTc                <  1:365679/151‑1 (MQ=255)                                                                                                                                  gTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCg               <  2:273081/149‑1 (MQ=40)                                                                                                                                        tGGCGGGGATCA‑CTCCCATAAGCG‑CG‑AACTTAAGGGGTGTGGTATGACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCACACGAAAACTGGTCAGCAATTCTGGCCCATC‑CTGGTAAGGCCGAAGAACTGGATAATTCGGccagta        >  1:432725/1‑145 (MQ=11)                                                                                                                                         ggCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt      >  2:48349/1‑151 (MQ=35)                                                                                                                                             gggATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg  <  1:263477/151‑1 (MQ=33)                                                                                                                                             gggATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg  <  2:469817/151‑1 (MQ=33)                                                                                                                                                |                                                                                                                                                          CCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCA‑CTCCCATAAGCG‑CT‑AACTTAAGGGTTGTGGTATTACG‑CCTGATATGA‑TTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATA‑TTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCG  >  NC_000913/15239‑15529

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