Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	55.4%	14.0 / 75.9	56	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (24/1);  new base C (0/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 5.74e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-08

GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCA  >  NC_000913/3556059‑3556350                                                                                                                                                            |                                                                                                                                                gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTc                                                                                                                                                >  1:805618/1‑150 (MQ=255)           ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAgtgt                                                                                                                                     >  2:489773/1‑151 (MQ=255)              cGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAg                                                                                                                                   >  1:739584/1‑150 (MQ=255)                          gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCt                                                                                                                       >  1:429762/1‑150 (MQ=255)                              tccccgagcgaaAAAAGAGGGTCCGTTAATATTTCCACTTGCG‑ACTACAGGCGAGACGTTATT‑‑CCCCG‑CCCCTGGTCCCCGCGGGGGTTTAGAAATTTTCTGTTC‑ATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                                  <  1:580640/143‑1 (MQ=255)                                 accgccaaaaaaaagGCTTTGGGTAATATTTTCATGTGCG‑ATGGAAAGAGAGAGGTTTTT‑CCCCGC‑CCCCGGGC‑CGCGCGGCTGGTTATAAAAATTCTGTAC‑AGGCGGGGAAAA‑T‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                               <  2:200865/142‑1 (MQ=255)                                            aaaaGCGTGGGGTAATTTTTTC‑CTTTCGAACTTCAAGCGAGCCGTTTTC‑CCCCGC‑CCCCTGTC‑CCCGGGGCAGTTTAAAAATATTCTTTTC‑AGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                    <  1:1018845/151‑1 (MQ=255)                                                  tgtggtttttttttcCACTTTCG‑ACTTAAGGCGAG‑CGGTTTTCCCCCGC‑CCCCGGGC‑CGCGCGGCAGGTTAGAATTTTTCTGTCC‑AGGAGGGGAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                              <  1:132735/140‑1 (MQ=255)                                                       gTGATTTTTCCCCTTGCG‑ACTTCAAGCGAGACGTTATC‑ACCCGG‑CCCCTGGC‑CGCGCGGCAGGTTATAAATTTTCTGTTC‑ATGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                          <  1:1031346/150‑1 (MQ=255)                                                       gTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                         >  2:28408/1‑151 (MQ=255)                                                        ttatatTTCCCCTTTCG‑ACTTCAAGCGAGACGTTATT‑CCCCGG‑CCCCTGTC‑CGCGCGGCAGTTTAGAAATTTTCTGTTC‑AGGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                        <  1:696074/149‑1 (MQ=255)                                                        tGATTTTTTCTCTTGCG‑ACTTCAAGCGAGACGTTATT‑CCCCGC‑CCCCTGGT‑CGCGGGGCAGGTTATGAATTTTCTGTTC‑ATGCAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                        <  1:608005/151‑1 (MQ=255)                                                         aattttttCGCTTGCG‑TCTACAAGCAGGACGTTTTT‑CCCGGC‑CCCCGGGT‑CGGGCGGCAGTTTATGAATTTTCTGTTC‑AGGGAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                       <  1:47712/147‑1 (MQ=255)                                                          atatTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                      >  2:866448/1‑151 (MQ=255)                                                           tttttccccTTTCC‑ACTTAAGGCGAGGCGTTTTT‑CCCCGC‑CCCCTGGC‑CGCGGGGCAGTTTATGAATTTTCT‑TTT‑AAGGAGGGAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                     <  1:489568/143‑1 (MQ=255)                                                           tatTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                    >  1:373475/1‑136 (MQ=255)                                                              tttCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCCCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                   >  1:947572/1‑150 (MQ=255)                                                                     tccG‑ACTTAAAGGGAGTCGTTTTT‑CCCCGG‑CCCCTGGG‑CCCGCGGCAGTTTAGAAATTTTTTGTCC‑AAGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                           <  2:186454/149‑1 (MQ=255)                                                                      tCG‑ACTTAAAGAGAGACGGTATT‑CCCCGG‑CCCCGGGG‑CGCGGGGCAGTTTATAAATATTCTGTTC‑ATGCGGGGAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                          <  2:44045/151‑1 (MQ=255)                                                                          ttttaaagagggagggTATG‑CCCGGG‑CGCCGGGG‑CGCGGGGCGGGTTCAGCTTCTTCTGTTC‑ATGGGGGGAAAA‑A‑CCCCCCGCCAGTGTAGCGCTTTCCCGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                          <  2:389532/134‑1 (MQ=255)                                                                            ttCAAGCGAGGCGGTATC‑CCCCGC‑CCCCGGGT‑CGCGGGGCGGTTTTAGATTATTCTGTTC‑AGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                     <  2:629559/151‑1 (MQ=255)                                                                              aaagggagggggTTTT‑CCCCGG‑CCCCGGGC‑CGCGCGGCAGGTTATGAATTTTCTGTCC‑AGGCAGGGAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                   <  2:495450/141‑1 (MQ=255)                                                                                 gCGAGAGGGTTTT‑CCCCGC‑CCCCCGGC‑CGCGCGGCAGTTTTAGAATTTTCTGTTC‑ATCCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACaa                                                                 <  2:514584/150‑1 (MQ=255)                                                                                    agcgggTATT‑CCCCGG‑CCCCTGCC‑CGCGCGGCGGTTTATAAATATTCTGTTC‑ATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                              <  1:407039/146‑1 (MQ=255)                                                                                    agccgttATT‑CCCCGC‑CCCTTGGC‑CGCGCGGCAGTTTATGAATATTCTGTTC‑ATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                              <  1:1000593/144‑1 (MQ=255)                                                                                    agagGGTATC‑CCCCGC‑CCCCGGGC‑CGCGCGGAGGGTTATAAATATTCTGTCC‑AGGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                              <  1:994573/150‑1 (MQ=255)                                                                                    acaCGGTATT‑CCCCGG‑CCCCTGGT‑CGCGCGGCGGTTTATAATTTTTCTGTTC‑ATGCAGGGAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                              <  1:523743/148‑1 (MQ=255)                                                                                         gTATT‑AGCTGG‑CCCCTGGT‑CTCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGCAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                         >  1:352773/1‑150 (MQ=255)                                                                                         gTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACCCCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                         >  1:578864/1‑150 (MQ=255)                                                                                         gTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                         >  2:13543/1‑150 (MQ=255)                                                                                         gTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                         >  1:1019294/1‑150 (MQ=255)                                                                                          tATG‑ACCCGG‑CCACTTGT‑CCCCGGGCTAGTTATTAATTTTCTTTTC‑CAGGGGGGAAAACC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                         <  1:596814/150‑1 (MQ=255)                                                                                           gtacaCCCGGCCCCCGGGG‑CCCGCGGCAGTTAAAAAAAATTCTGTCC‑AGGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCCCGCCAGCTAATAGAGCTTCTCCACGCATCTGCACAATGTTCCcg                                                         <  2:449216/146‑1 (MQ=255)                                                                                               aCCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATCCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                    >  1:489931/1‑150 (MQ=255)                                                                                                     tccccAGGT‑CGCGCGGCAGGTTATGAATATCCTGGTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                              >  2:9382/2‑151 (MQ=255)                                                                                                      ccccTGGC‑CGCGCGGC‑GGTTAATAAATTTTTGGTC‑CAGGAGGGAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                             <  2:455420/151‑1 (MQ=255)                                                                                                                      cggtttttaaaTTTTCTGTTC‑AGGAAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                              <  2:858459/142‑1 (MQ=255)                                                                                                                       ggtttttaaaTATTCTTTTC‑ATG‑AGGAAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                              <  2:435447/142‑1 (MQ=255)                                                                                                                       ggtttATGAATTTTC‑GTTCAAGGGAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                              <  2:1034521/147‑1 (MQ=255)                                                                                                                       gggTTATAAATATTCTTTTC‑ATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                              <  2:70765/149‑1 (MQ=255)                                                                                                                        ggTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                            >  2:198999/1‑151 (MQ=255)                                                                                                                         tttATGAATTTTCTTTTC‑ATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                             <  2:255692/148‑1 (MQ=255)                                                                                                                         tttATAAATTTTCTTTTC‑GTGCAGGGAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                            <  2:935728/149‑1 (MQ=255)                                                                                                                         gTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                           >  1:306774/1‑151 (MQ=255)                                                                                                                         gTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                           >  2:874323/1‑151 (MQ=255)                                                                                                                         gTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                           >  1:613418/1‑151 (MQ=255)                                                                                                                         gTTAAGAATTTTCTGTTC‑AGGGGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                              <  2:646268/148‑1 (MQ=255)                                                                                                                          ttATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                          >  1:995941/1‑151 (MQ=255)                                                                                                                            aTGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                         >  2:285899/1‑150 (MQ=255)                                                                                                                             tGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                       >  2:810597/1‑151 (MQ=255)                                                                                                                                atatTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                    >  2:528287/1‑151 (MQ=255)                                                                                                                                        ttccaggaaGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc             <  1:942303/143‑1 (MQ=255)                                                                                                                                        tcc‑ATGCAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc             <  1:323467/148‑1 (MQ=255)                                                                                                                                              gCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCg        <  1:809516/150‑1 (MQ=255)                                                                                                                                                   gAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGGTGCCACGGGCGAGACTTCCGTTGcca  >  2:863432/1‑151 (MQ=255)                                                                                                                                                   gAAAA‑A‑CTCCACGCCAGTGTAGCGATTTCACGCTAAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTTCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcca  >  2:829118/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                GCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCG‑ACTTCAAGCGAGACGGTATT‑ACCCGG‑CCCCTGGT‑CGCGCGGCAGGTTATGAATATTCTGTTC‑ATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCA  >  NC_000913/3556059‑3556350

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