BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,842,368	A→C	28.3%	V12G (GTG→GGG)	yicN ←	DUF1198 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,842,368	0	A	C	28.3%	60.8 / 4.3	48	V12G (GTG→GGG)	yicN	DUF1198 family protein
Reads supporting (aligned to +/- strand): ref base A (17/15); new base C (2/14); total (19/29)
Fisher's exact test for biased strand distribution p-value = 1.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGT  >  NC_000913/3842143‑3842586

cTACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGAACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAACCGCTTTTCTGGCCCCGGATGTCGGCACCCCAAAACCAACCACAAacacc                                                                                                                                                                                                                          >  2:19653/1‑227 (MQ=255)

 tACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTTATCGACCATCGAATCCCCGGGTACAACATCGATGTTCGACCGATCGCTGCGACGGCGAATCGCTTTTCTGGCCCCGGATGTCCGCACCCACAAACCACCCACAAacaaaa                                                                                                                                                                                                                         >  2:130150/1‑225 (MQ=255)

      tCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGCTGTTCAGCCTATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGTAGTCACCACCCGAAAAACAACcacacacaaacccccc                                                                                                                                                                                                                    >  2:404223/1‑220 (MQ=255)

      tCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGAACGACATCGATGTTCAGCCGATCGCTGAGCGGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAACCCTACCACAAACACTACCTcc                                                                                                                                                                                                                    >  1:645375/1‑228 (MQ=255)

      tCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAAAGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGAACGACATCGATGTTCAGCCGATCGCTGCGACGGCGAATCGCTTTTCTGGCCCCGGACGTCAGCACCCGAAAAACAACCACAAACACCACCGcc                                                                                                                                                                                                                    >  2:402797/1‑228 (MQ=255)

             tgacaatTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTg                                                                                                                                                                                                              >  1:494833/1‑227 (MQ=255)

             tgacaatTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGAGCTCAGCACCCGAAAACCAACCACAACCACCACcccccacgcgc                                                                                                                                                                                                             >  2:420226/1‑221 (MQ=255)

              tccaatTTGCCAAATG‑GCAACCCCCTGCGGGGGTTTTTGCGGGTGGGCTCCCACCGGACGATGAAAATAACGTAAAAATCCGTCACCAGCTGATTTCCCCATTGGAGCGACCCTCGACCCCACGGGTACGCCATCGATGTTCGGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCAAAACCCAACCACAACCCCTACCGCCAGCGTGGc                                                                                                                                                                                                            <  2:284827/226‑1 (MQ=255)

                        cAGATGAGCAA‑CACCTGCGTGGGGTTTTGCAGGGGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATc                                                                                                                                                                                                  <  1:391647/228‑1 (MQ=255)

                           aTGAGCAA‑CACCTGCGCGGCTTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAg                                                                                                                                                                                               <  1:386398/228‑1 (MQ=255)

                                cAA‑CCCCTGCGCGGCTTTTTGCAGGTGCAACTCATCCGGCCGATGCAAAAAACGAAAAAATTCGTCCCCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCGGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCAt                                                                                                                                                                                          <  1:456010/228‑1 (MQ=255)

                                 aa‑CACCTGACGGGCTTTTTTCAGGTGCGACTCATCCGGACGATGCAAAAACCGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGACCGACCATCGACTCCACGGGTACGCCATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATg                                                                                                                                                                                         <  1:219789/228‑1 (MQ=255)

                                      ccggggcggcgTTTTGCAGGTTCAGCTCACCCGTGCGATGCAAATAACGTAAAAATCTGTCCCCAGCCGACTTTCCCATTTGTTCTCCCATCGACTCCACGGTGACGACATCGATGTTCAGCCGATCGCGGAGACGCCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGtt                                                                                                                                                                                     <  2:30279/223‑1 (MQ=255)

                                          cgcgGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCCACCGAAAACCAACCATAAACACAACCGCCAGCGGGCGGGGCATTAACCAGATCATGCATTCTcc                                                                                                                                                                                 >  2:521399/1‑228 (MQ=255)

                                          cgcgGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGa                                                                                                                                                                                                          >  1:12106/1‑203 (MQ=255)

                                           gcggcgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTCGCCCCGGATGTCAGCACCCGACAACCAACCACAAACACTCCCGCCAGCGTGGCGAGCATTTTCCAGCTCATGACtttccct                                                                                                                                                                                >  1:163099/1‑223 (MQ=255)

                                               cgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGCATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCCACCACCAAACCTACCGCCAGCGGCGCGAGCATTATCACGATCATGGGTTCCCCTGCgg                                                                                                                                                                            >  2:345647/1‑227 (MQ=255)

                                               cgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGt                                                                                                                                                                            >  1:177478/1‑228 (MQ=255)

                                                  tttGCGGGTGCGCCTCTCCCGGATGTTGCAAATAAAGTAAAAATTCGTCACCCGCTGACTTTCCCATTTGACCGACCATCGACTCCGTGGGTACGAACTCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAACCCAACCACAACCCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAAt                                                                                                                                                                         <  2:177478/228‑1 (MQ=255)

                                                    ggaagggtcgccTCAACCGGAACAGGCAAATAACGTAAAAATTGGGCACCAGTTGGCTTTCCCATTTGACCGCCCATCGCTCCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCa                                                                                                                                                                       <  1:586425/217‑1 (MQ=255)

                                                       aGGTGCGACTCACCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTCCCCATTTGATCGACCATCGACGCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGaa                                                                                                                                                                    <  1:463382/228‑1 (MQ=255)

                                                                   tCCGGACGATGCCAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGCCTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGAACCCGAAAAACAACCAAAAACCATACCGCCAGCGTGGCGAGCATTATCACGATCATGAATTCACCTCCGTAATCACAACtactagcatc                                                                                                                                                         >  2:564362/1‑218 (MQ=255)

                                                                            ggaaaaTAACGTAAAAATTCGTCACCAGCTGTCTTTCCAATTTGAGCGACCAGCGACTCCAGGGGTACGCAATCGATGTTCGCCAGATCGCGGAGCCGGGAAATCGCTTTCCGGGCCCCGGAGGTCAGCACCCAAAAACCAACCCAAACCCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCt                                                                                                                                               <  2:645375/225‑1 (MQ=255)

                                                                               aaaCAACGTAAAAATTCGTCACCAGCTGCTTTTCCCATTTGAGCGACCATCGACTCCCGGGTTACCACATCGATGTGCAGCCGATCGCTGAGACGGCGAATCGCTGTTCTGGCCCGGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGc                                                                                                                                            <  1:586913/228‑1 (MQ=255)

                                                                                      gTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCAAAACCATAAGTTCTCCTGCGGAATCAGAACTCAAGATATTGTGAAACCCCCGCCaagtac                                                                                                                                     >  2:589166/1‑223 (MQ=255)

                                                                                       tAAAAAGTCGTCCCCAGCTGATTTTCCCATTTGATCGACCATCGACCCCACGGGTACGACATCGATGTTCAGCCGACCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCACCCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTAc                                                                                                                                     <  1:19653/227‑1 (MQ=255)

                                                                                               cGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGCCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTATTccac                                                                                                                             >  1:492123/1‑227 (MQ=255)

                                                                                                  tccccgCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGa                                                                                                                         <  1:125038/223‑1 (MQ=255)

                                                                                                           gcTTTCCCATTTGATCGCCCTTCGACCCCACGGGTACGACCTCGATGGTCCGCCGATCGTTGAGTCGGCGAATCGGTTTTCTGGCCCCGGGTGCCAGCACCCGAAAACCACCCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTc                                                                                                                <  2:7343/227‑1 (MQ=255)

                                                                                                               tCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCa                                                                                                               >  1:355055/1‑225 (MQ=255)

                                                                                                                     ttGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACa                                                                                                      >  1:11110/1‑228 (MQ=255)

                                                                                                                         tCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACagag                                                                                                   >  1:643612/1‑227 (MQ=255)

                                                                                                                               agcggcTCCACGGGTACGACATGGATGTCCAGCCGATCAGTGAAACGGCGAATCGCTTTTCTGGCCCCGGATGTACGCACCCAAAAACCAACCACAATCACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGac                                                                                            <  1:556242/223‑1 (MQ=255)

                                                                                                                                  caCTCCAGGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCATAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccac                                                                                         <  1:521399/227‑1 (MQ=255)

                                                                                                                                   aCTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCCGAATAATCTTCAGTGCAAACAGACCTGCGaccacc                                                                                        >  1:50464/1‑228 (MQ=255)

                                                                                                                                           ggtaacgactCGCTGTCCCGCCGATCGCTGAGGCGGGGAATCGTTTTGCTGGCCCGGGATGTCAGCCCCCGAAAACCAGCCACAAACCCTCCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCCCCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                <  2:435033/219‑1 (MQ=255)

                                                                                                                                           ggatcGACATCGATGTTCAGCCGTTCGCTGAACCGGCTAATCGCTTTTCTGGCCCCGGATGGCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                <  1:589166/224‑1 (MQ=255)

                                                                                                                                              aCGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCCCcacaacaa                                                                              >  2:528922/1‑223 (MQ=255)

                                                                                                                                                                   tcccTGAGACGGCATATCGGTTTTCCTGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                         <  1:404223/224‑1 (MQ=255)

                                                                                                                                                                   ccGCTGAGCCGGCAAATCGCTTTTCTCGCCCCGGATGTCCGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                         <  1:402797/226‑1 (MQ=255)

                                                                                                                                                                         agccGGCGAATCGCTTTTCGGGCCCCGGATTTCAGCACCCGAAACCCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                  <  1:653044/225‑1 (MQ=255)

                                                                                                                                                                                 acccggttttCGGGCCCCGGTTGGCAGCCCCCAAAACCCAAACAAAAACCCCACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAACCCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCa                                          <  2:50464/222‑1 (MQ=255)

                                                                                                                                                                                       ttttCTGCCACCGGATGTGACCCCCCGAAAACCAACCCCAAACCATACCGCAAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatgatg                                     <  1:326518/227‑1 (MQ=255)

                                                                                                                                                                                       ttttCGGGCCCCGGATGTCAGCACCCGAAAACAAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAGCCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatgatg                                     <  1:432490/227‑1 (MQ=255)

                                                                                                                                                                                           cgggCCCCGGGAGGTAGCCCACGAAAACCAACCACAAACCCTACCGCCAGCGTGGCGAGAATTAACCAGATCATGAGTTCTCCTGCATAATCAGAACTCATGATAGCTGGAAACCCGCGCCAGTTTCTAAGCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCTCCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTa                                <  2:643198/226‑1 (MQ=255)

                                                                                                                                                                                           cTGGCCCCGGATGTCAGCACCCGAAAACCAACCCCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTa                                <  1:420226/228‑1 (MQ=255)

                                                                                                                                                                                            aggCCCGGGATGTCAGCACCCGAAAACCAACCCCAAACCCCACCCCCAGCGTGGCGAGCATTAACCAAATCACTAGTTCTCCTGCGTTATCAGAACCCATTATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCCTCAGTGCAAACAGAACTGGGACCACCACCACACACAGGTTCAAATCGCGCCAGCGCCCGGTGCATGCTTTCATGATGCAGTa                                <  1:530062/226‑1 (MQ=255)

                                                                                                                                                                                                       ttcaaccaccGAAAACCACCCAAAAACCCTCCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTTCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCTCACTCAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGcc                    <  1:102445/220‑1 (MQ=255)

                                                                                                                                                                                                             cccaaaaaaCCACCCCAAAACCCTCCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc              <  2:11110/223‑1 (MQ=255)

                                                                                                                                                                                                               ccGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCaa            >  1:29893/1‑228 (MQ=255)

                                                                                                                                                                                                                         accaCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATCATTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGt  >  1:356091/1‑228 (MQ=255)

CTACCATCGACAATGACAATTTGCCAGATGAGCAA‑CACCTGCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGT  >  NC_000913/3842143‑3842586

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