BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,200,006:1	+C	5.8%	coding (829/1278 nt)	PP_1944 →	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,200,006	1	.	C	5.8%	430.5 / 15.0	104	coding (829/1278 nt)	PP_1944	aminomethyltransferase
Reads supporting (aligned to +/- strand): ref base . (49/49); new base C (3/3); total (52/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGT  >  NC_002947/2199863‑2200153

cTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgc                                                                                                                                               <  1:1407530/150‑1 (MQ=255)

  gACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctg                                                                                                                                             <  2:1454425/150‑1 (MQ=255)

   aCTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                            <  1:39678/150‑1 (MQ=255)

    cTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctggg                                                                                                                                           >  1:2497097/1‑150 (MQ=255)

      aCATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAGGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggct                                                                                                                                         <  1:1471488/150‑1 (MQ=255)

        aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgggctgg                                                                                                                                       <  2:1367147/150‑1 (MQ=255)

          gCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCtgggctgggc                                                                                                                                     <  2:1756226/150‑1 (MQ=255)

           cccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGcc                                                                                                                                    >  2:2317061/1‑150 (MQ=255)

             cgcgATTTTCGACGCAGAGAGGGCGGGTTATCTAAAGGAATTCCGCGCGGGTATGCCAAGATTTACCTCGCCTTTAGATGTTGCGGGGAGCTTTGAGTAGAACAAAGTAAGAGATTGATATGGCTGGCCAG‑TCGTGctgagctgggataa                                                                                                                                  >  2:1032856/1‑146 (MQ=255)

             cgcgATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCaa                                                                                                                                  >  1:1151610/1‑150 (MQ=255)

             cgcgATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCaa                                                                                                                                  >  1:396318/1‑150 (MQ=255)

                  ttttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCa                                                                                                                             >  2:929787/1‑150 (MQ=255)

                        aCCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGTAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCa                                                                                                                             <  1:1581006/144‑1 (MQ=255)

                        aCCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGTAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCa                                                                                                                             >  2:1581006/1‑144 (MQ=255)

                          ccAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGATCTTTGAGTCGAATAAGGTAAGCGATTGGGATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATGAAGTTc                                                                                                                     >  2:633315/1‑150 (MQ=255)

                          ccAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                     >  1:1686151/1‑150 (MQ=255)

                          ccAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                     >  1:1802007/1‑150 (MQ=255)

                           cAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCg                                                                                                                    >  1:2554573/1‑150 (MQ=255)

                           cAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCg                                                                                                                    >  1:344660/1‑150 (MQ=255)

                                aTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCAc                                                                                                               >  2:2804224/1‑150 (MQ=255)

                                             tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGAc                                                                                                            >  1:1192309/1‑140 (MQ=255)

                                             tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGAc                                                                                                            <  2:1192309/140‑1 (MQ=255)

                                             tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                  >  2:1302796/1‑150 (MQ=255)

                                             tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                  >  1:2300799/1‑150 (MQ=255)

                                              cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCg                                                                                                 >  1:277954/1‑150 (MQ=255)

                                              cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCCCAGCCCAGCTCGa‑ctgggctgcgata                                                                                                                                   <  1:1966/116‑6 (MQ=255)

                                              cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCCCAGCCCAGCTCGa‑ctgggctgcgata                                                                                                                                   >  2:1966/1‑111 (MQ=255)

                                               aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGt                                                                                                >  1:2063913/1‑150 (MQ=255)

                                               aaGGAATTCCGCGCGGGGATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGACAAGACAAGGAAGTGCGATCACGACTTTATCGGTCGt                                                                                                >  1:1735714/1‑150 (MQ=255)

                                                aGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTg                                                                                               <  2:2063913/150‑1 (MQ=255)

                                                 ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGTAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCAt                                                                                                                            <  1:89295/120‑1 (MQ=255)

                                                 ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGTAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCAt                                                                                                                            >  2:89295/1‑120 (MQ=255)

                                                 ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑Tc                                                                                                                                                   >  2:561054/1‑97 (MQ=255)

                                                    aTTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCcgc                                                                                           <  1:2950013/150‑1 (MQ=255)

                                                       ccGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                        <  2:1802007/150‑1 (MQ=255)

                                                             gggTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaaaa                                                                                  >  1:1861907/1‑150 (MQ=255)

                                                              ggTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaaaaa                                                                                 <  2:1151610/150‑1 (MQ=255)

                                                               gTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGGCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAag                                                                                <  1:2507733/150‑1 (MQ=255)

                                                                 aTGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagag                                                                              <  1:602349/150‑1 (MQ=255)

                                                                    ccAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                     >  2:1346743/1‑108 (MQ=255)

                                                                    ccAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                     <  1:1346743/108‑1 (MQ=255)

                                                                     cAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                   <  1:2703623/109‑1 (MQ=255)

                                                                     cAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                   >  2:2703623/1‑109 (MQ=255)

                                                                     cAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTg                                                                          >  1:656204/1‑150 (MQ=255)

                                                                           ttGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGc                                                                    <  1:3045385/150‑1 (MQ=255)

                                                                            tGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGcc                                                                   <  2:2300799/150‑1 (MQ=255)

                                                                               cTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCg                                                                >  2:2544083/1‑150 (MQ=255)

                                                                                tCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGa                                                               <  1:2317061/150‑1 (MQ=255)

                                                                                  gACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGCGCTGGGCCAAAACCATCAAGTTCGATCAGCACTTTATCGGTCGTGCAGCGTTGTAAAAACAGGTTGCAGAGCCTCGACg                                                             <  2:396318/150‑1 (MQ=255)

                                                                                    cTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACggg                                                           <  1:929787/150‑1 (MQ=255)

                                                                                     tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATc                                                                                                                 >  2:2291592/1‑95 (MQ=255)

                                                                                     tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGa                                                                                                             <  1:2323054/99‑1 (MQ=255)

                                                                                     tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGa                                                                                                             >  2:2323054/1‑99 (MQ=255)

                                                                                     tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                          >  1:252079/1‑150 (MQ=255)

                                                                                      ttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAGAGAGGTTGCAGAGCCCCGACGGGTg                                                         >  1:1006715/1‑150 (MQ=255)

                                                                                      ttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTg                                                         >  1:1007996/1‑150 (MQ=255)

                                                                                         aaTGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTTCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATg                                                      >  1:1096354/1‑150 (MQ=255)

                                                                                           tGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATgcg                                                    >  2:1162621/1‑150 (MQ=255)

                                                                                              tGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCAc                                                 <  2:277954/150‑1 (MQ=255)

                                                                                                  gggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                             <  2:2497097/150‑1 (MQ=255)

                                                                                                  ggcaGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGCTCGATCACGAATTTATCGGTCGTGCCNCGTTGCAAANAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                             <  1:1032856/147‑1 (MQ=255)

                                                                                                     aGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGtt                                          <  1:1794485/150‑1 (MQ=255)

                                                                                                     aGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGtt                                          <  1:317697/150‑1 (MQ=255)

                                                                                                       cTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                        <  2:1007996/150‑1 (MQ=255)

                                                                                                       cTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCCCGACGGGTGATGCGCACCTTGGTTTg                                        <  2:1006715/150‑1 (MQ=255)

                                                                                                        tttGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCTAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagag                                                                              >  1:2242452/1‑111 (MQ=255)

                                                                                                        tttGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagag                                                                              <  2:2242452/111‑1 (MQ=255)

                                                                                                        tttGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTgg                                       >  1:2788629/1‑150 (MQ=255)

                                                                                                            ggTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATg                                   >  1:2432392/2‑150 (MQ=255)

                                                                                                            aGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCaa                                                                                                                         >  2:254405/1‑64 (MQ=255)

                                                                                                            aGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCaa                                                                                                                         <  1:254405/64‑1 (MQ=255)

                                                                                                            aGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATg                                   <  1:1162621/150‑1 (MQ=255)

                                                                                                            aGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCAGCGTTACACTAACTAGCAGCATTGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATg                                   <  1:1302796/150‑1 (MQ=255)

                                                                                                             gTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATgg                                  >  1:423267/1‑150 (MQ=255)

                                                                                                               cGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGGCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGGCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATAGtg                                <  2:250639/150‑1 (MQ=255)

                                                                                                               cGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtg                                <  1:2544083/150‑1 (MQ=255)

                                                                                                                 aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAACACCATCAAGTTCGATCACGACttt                                                                                                         >  1:1754252/1‑75 (MQ=255)

                                                                                                                 aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACttt                                                                                                         <  2:1754252/75‑1 (MQ=255)

                                                                                                                 aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGttt                                         >  1:951546/1‑139 (MQ=255)

                                                                                                                  aTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgatg                             <  1:1939028/150‑1 (MQ=255)

                                                                                                                   tAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGaa                                     >  1:3038186/1‑141 (MQ=255)

                                                                                                                   tAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGaa                                     <  2:3038186/141‑1 (MQ=255)

                                                                                                                         aaGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGc                      <  2:423267/150‑1 (MQ=255)

                                                                                                                          aGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTTCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCa                     <  2:1096354/150‑1 (MQ=255)

                                                                                                                          aGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCCCCTTGGTTTGGAATGGTGATGACGTTGCa                     <  2:1735714/150‑1 (MQ=255)

                                                                                                                             gATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGAt                  >  2:1816142/1‑150 (MQ=255)

                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGAt                                                                                                                  >  1:2650151/1‑49 (MQ=255)

                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGAt                                                                                                                  <  2:2650151/49‑1 (MQ=255)

                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCa             >  1:1603051/1‑150 (MQ=255)

                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCa             >  1:420928/1‑150 (MQ=255)

                                                                                                                                  gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCa             >  2:2304367/1‑150 (MQ=255)

                                                                                                                                   tATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGGGATGACGTTGCAGATGTTCAt            <  2:2788629/150‑1 (MQ=255)

                                                                                                                                    aTCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATg           <  1:2804224/150‑1 (MQ=255)

                                                                                                                                       gCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCt        >  1:2866725/1‑150 (MQ=255)

                                                                                                                                       ccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc          >  1:1164707/2‑148 (MQ=255)

                                                                                                                                       ccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc          <  2:1164707/147‑1 (MQ=255)

                                                                                                                                        cAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACCTTAGAGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTc       >  1:561339/1‑150 (MQ=255)

                                                                                                                                            ccAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTg   >  1:618778/1‑150 (MQ=255)

                                                                                                                                             cAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGt  >  1:148962/1‑150 (MQ=255)

                                                                                                                                             cAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGt  <  2:252079/150‑1 (MQ=255)

CTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGT  >  NC_002947/2199863‑2200153

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