Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,639,556	A→T	6.7%	intergenic (+26/+34)	creD → / ← arcA	putative inner membrane protein/DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,556	0	A	T	6.7%	166.9 / 2.3	65	intergenic (+26/+34)	creD/arcA	putative inner membrane protein/DNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base A (37/23);  new base T (2/3);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGT  >  NC_000913/4639414‑4639697                                                                                                                                                |                                                                                                                                               gcgcTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaa                                                                                                                                                >  1:996151/1‑143 (MQ=255)   gcTGGCCGGCATGATGTTTGGGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTGACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaaaa                                                                                                                                              <  1:524562/143‑1 (MQ=255)    cTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAACGAAGTTACAACGGACGATG‑AGTTACGTATCTGGCATTAAGGTCGAAAAATCAAAACGGCTCTAAAAAg                                                                                                                                             >  1:47919/1‑143 (MQ=255)         cGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCg                                                                                                                                        >  1:363212/1‑143 (MQ=255)         cGGCATGATGTTTGTGACCCGTAATATCGACTCAACTGCATTAGCACTGCAGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AAGTACGGATCAGGAAATGAGGGTGAAAAATAAAAACGGAGCTAAAAAGCAcca                                                                                                                                        >  2:789676/1‑142 (MQ=255)          ggCATGATGTTTGGGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑t                                                                                                                                      <  1:669847/143‑1 (MQ=255)            cATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGATTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttt                                                                                                                                    >  1:344539/1‑143 (MQ=255)                  gttttagACCCCTAGTATCTGTTGGTATGCGGTTGCACGGACGATAATGACAGCCAGTAAAATAGTTAGGACGGACGATG‑GGTTACGGATCTGGACATAAGTTTGTTTAATTAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTg                                                                                                                              <  2:395197/137‑1 (MQ=255)                   tttGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCG‑tttttatttt                                                                                                                             <  1:663496/143‑6 (MQ=255)                   tttGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑GTTTTTTTGa                                                                                                                             >  1:233391/1‑143 (MQ=255)                    ttGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAACATGAAAGCCAGTGAAGAAGTTACAAAGGAGGGTT‑ATGTGTGTATATGGTAAGACGGGTGAAAAATAAAAGGGGCGCTATAGAGCGCGG‑TTATGTTTGAc                                                                                                                            >  2:784230/1‑143 (MQ=255)                      gtgACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACgg                                                                                                                          <  2:935499/143‑1 (MQ=255)                        gACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACggtg                                                                                                                        <  1:1104342/143‑1 (MQ=255)                         aCCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACggtgg                                                                                                                       <  1:917326/143‑1 (MQ=255)                            cGTAATATCGACTGGTATGCGTTTTCACTGCCTCATTTCAACCCCAGTATAGAAGTTACTACTGCCCATG‑AGTTACTTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACggtggcac                                                                                                                    >  2:574153/1‑140 (MQ=255)                            cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTCAAGAAGTTGCAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAACAGCGCAG‑GTTTTTTTGACGGTGGTag                                                                                                                    >  2:941125/1‑142 (MQ=255)                            cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACCACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTaa                                                                                                                    >  2:278333/1‑143 (MQ=255)                            cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAAAGGCGCTAAAAAGCGCCG‑TTTTTTTTAAAGGTGGTaa                                                                                                                    >  2:818908/1‑143 (MQ=255)                               aaTATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGc                                                                                                                 <  1:595172/143‑1 (MQ=255)                                gtgtCGACGGGTAGGCGCTTTCAATGCCGAAAATGAGAGCCAGTAAAAAAGTTATAACGGACGATG‑AGTTACGTATCTGGAAGTAGGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGcc                                                                                                                <  2:508904/140‑1 (MQ=255)                                atatCGACTAGTATTCTTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATC‑AGTTACGTATCTTGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGcc                                                                                                                <  1:565081/143‑1 (MQ=255)                                atagCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGCTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGcc                                                                                                                <  1:581008/143‑1 (MQ=255)                                 tatCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGCCGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTCACGCTAGTAAACCCg                                                                                                               >  1:785783/1‑143 (MQ=255)                                  atCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTACAGCCGa                                                                                                              >  2:708425/1‑143 (MQ=255)                                  atCGACTGGGATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTTGTAAAGCCGa                                                                                                              <  1:577436/143‑1 (MQ=255)                                   tCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGac                                                                                                             >  2:513263/1‑142 (MQ=255)                                       ttGGTTTGCTGTTTAACGGCCGAAAATGAAA‑TCAATAAAGGAATTTAAAAAGAGCATGAAGTTAGGTATTTGAAAGTAGGTTTAAAAATTAAAAACGGCGCTAAAAGGCCCGT‑TTTTTTTTGACGGTGGTAAAGCCGATTAAt                                                                                                         <  2:917215/142‑1 (MQ=255)                                        tGGTATGCGTTTGCACTGCCGAAAATGAAAGCCAGGAAAGAAGTTACAAGGGACGATG‑AGTTACGTAGCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATc                                                                                                        <  1:902993/143‑1 (MQ=255)                                              gCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACAGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGGGGTCAAGCCGATTAATCTTCCat                                                                                                  >  1:89279/1‑142 (MQ=255)                                               cGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGCCAACTAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGa                                                                                                 <  1:42671/143‑1 (MQ=255)                                                gTTTTCACTACCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGAt                                                                                                >  1:1006442/1‑143 (MQ=255)                                                     cACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCACATCACCg                                                                                           >  1:817433/1‑143 (MQ=255)                                                      aCTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTAAGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGCCGGTGGTAAAGCCGTTTAATCTTCCACATCActc                                                                                           >  1:557571/1‑141 (MQ=255)                                                        tGCCGAAAATGAAGGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAATAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAg                                                                                        <  2:457743/143‑1 (MQ=255)                                                          ccGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAgca                                                                                      >  1:328429/1‑141 (MQ=255)                                                          ccGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGaa                                                                                      >  1:863803/1‑143 (MQ=255)                                                            gAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGAGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGACGATTAATCTTCCAGATAACCGCAGAAGc                                                                                    >  1:372815/1‑143 (MQ=255)                                                              aaaTGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGATTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGa                                                                                  <  2:344539/143‑1 (MQ=255)                                                                     aGCCAGGAAAGAAGTTAAAACGGACGAGG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCtt                                                                           <  2:996151/143‑1 (MQ=255)                                                                        cAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCAc                                                                        >  2:48526/1‑143 (MQ=255)                                                                             aagaagTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGa                                                                   >  1:226994/1‑143 (MQ=255)                                                                                  gTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGCAcggt                                                              >  2:448788/1‑143 (MQ=255)                                                                                   ttACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAtggtg                                                             <  2:1006442/143‑1 (MQ=255)                                                                                       tagggACGACG‑AGTTCCCGATCTGGAAGTAAGCTTGGAGAATAAAAACGGCGCTAAAAAGCTCCG‑CTTTTTTGGACGGTGGGAAAGCCGATTAATCTTCCAGATAACCGCTGAAGCTTGAACCTTCACCGTGAATGGTGGCga                                                         <  2:395993/140‑1 (MQ=255)                                                                                       aaCGGACGATG‑AGTTACGTATCCGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCga                                                         >  2:236478/1‑143 (MQ=255)                                                                                          gggcGATG‑AGTGACGTGTGTGGGAATAAGGTCGAATGATAAACACGGCTCTGAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCCATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAAGGTGGCgatga                                                      <  2:147910/140‑1 (MQ=255)                                                                                           gACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatgat                                                     >  2:103031/1‑143 (MQ=255)                                                                                                    gatACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGCTTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccgccg                                             >  2:739484/3‑143 (MQ=255)                                                                                                          tATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccg                                       >  2:881906/1‑143 (MQ=255)                                                                                                           aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAATAATCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                      <  2:997101/143‑1 (MQ=255)                                                                                                             cTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCCCCGTGGATGGTGGCGATGATTtccggcgtatccggcg                                    >  1:765267/1‑143 (MQ=255)                                                                                                             cTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggcg                                    >  1:258738/1‑143 (MQ=255)                                                                                                               ggAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggcgta                                  >  2:850331/1‑143 (MQ=255)                                                                                                                      aGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGa                           <  1:103031/143‑1 (MQ=255)                                                                                                                         ttGAAAAATAATAGCGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTTAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAAt                        <  1:659139/143‑1 (MQ=255)                                                                                                                         ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGCCTAATCTTCTAGATTACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAAt                        >  1:417320/1‑143 (MQ=255)                                                                                                                         ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGACTGGCGGCGATGACTTCCGGCGTATCCGGCGAAGATTCGAGAt                        >  2:787206/1‑143 (MQ=255)                                                                                                                              aaaTAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAACGGTGGCGATGATTTCCGGCGTATCCTGCGTAGATCCGTACTGTCTa                   >  1:26889/1‑143 (MQ=255)                                                                                                                               aaTAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                  <  1:271134/143‑1 (MQ=255)                                                                                                                                 tAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTGGATTCGAAATGTGTACGa                >  1:742246/1‑143 (MQ=255)                                                                                                                                 tAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCCTAGCTTCGAAATGTTCACGa                >  1:487194/1‑143 (MQ=255)                                                                                                                                 aaaaaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGAACACGCGCTTACCTTCACCGAGCAAGATGCCGATAATTTCCGCCGTCTCCGGCGTCGATTCGGAATGTTTACGa                >  1:863040/2‑143 (MQ=255)                                                                                                                                    aaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGACATGTTTACGAATa             >  1:998966/1‑143 (MQ=255)                                                                                                                                    aaaCGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATa             <  1:278333/143‑1 (MQ=255)                                                                                                                                            cTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAt     <  1:881906/143‑1 (MQ=255)                                                                                                                                             tAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTCCGGCTACGGCGGATc    >  2:761135/1‑143 (MQ=255)                                                                                                                                             tAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGGATCTGGTGTAGACACGAAATGTTTACGAATACGGCGGCTc    >  2:829184/1‑143 (MQ=255)                                                                                                                                               aaaaGCGCCG‑TTTTTTTTGTCGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACATTCACCGTTAATGGTGGACATGATTTCCGGCGTATCCGGCGTAGATTCGATATGTTCACGAATCCGGCGGAGCgg  >  2:34763/1‑142 (MQ=255)                                                                                                                                                |                                                                                                                                               GCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATG‑AGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGT  >  NC_000913/4639414‑4639697

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