BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,061,041	T→G	6.7%	R64R (CGT→CGG)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,061,041	0	T	G	6.7%	454.3 / 22.3	134	R64R (CGT→CGG)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base T (60/65); new base G (6/3); total (66/68)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 134 total aligned reads displayed.

CGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGC  >  NC_000913/3060898‑3061188

cGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCa                                                                                                                                                 <  4:723703/148‑1 (MQ=255)

 gTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCAc                                                                                                                                                <  3:964174/148‑1 (MQ=255)

   cGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACTGCGGGACGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGCTAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGTGTGCCACTa                                                                                                                                              <  4:688581/148‑1 (MQ=255)

   cGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                                                                                                                              <  4:839773/148‑1 (MQ=255)

   cGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAt                                                                                                                                             <  2:288358/149‑1 (MQ=255)

   cGGGAAAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAt                                                                                                                                             <  1:817219/149‑1 (MQ=255)

         aaaaCTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGGTAATCTGGAGGTGACAGGTACCCTTCCGGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAc                                                                                                                                       <  1:816374/149‑1 (MQ=255)

          aaaCTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAcc                                                                                                                                      >  1:689441/1‑149 (MQ=255)

          aaaCTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAc                                                                                                                                       >  4:368282/1‑148 (MQ=255)

               gTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGccc                                                                                                                                  >  4:515947/1‑148 (MQ=255)

                     tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGt                                                                                                                            <  3:633582/148‑1 (MQ=255)

                     tCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTg                                                                                                                           <  1:243464/149‑1 (MQ=255)

                      cGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTgg                                                                                                                          <  1:430978/149‑1 (MQ=255)

                      cGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTg                                                                                                                           <  4:927571/148‑1 (MQ=255)

                       gCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                                                                                                                         >  2:871426/1‑149 (MQ=255)

                        cTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGAc                                                                                                                        <  2:566294/149‑1 (MQ=255)

                        cTGGTTCATCAAACCGCGGAAGGGATCGCCATAAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGa                                                                                                                         >  3:670505/1‑148 (MQ=255)

                               atcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCACGCGGGCCAcgaa                                                                                                                                             >  1:705288/1‑118 (MQ=255)

                               atcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCACGCGGGCCAcgaa                                                                                                                                             <  2:705288/121‑4 (MQ=255)

                                tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGta                                                                                                                                          >  2:283222/1‑123 (MQ=255)

                                tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                                >  1:1056382/1‑149 (MQ=255)

                                tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                                >  1:602466/1‑149 (MQ=255)

                                tcaAACCGCGGAAGGGAGCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGGCAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGt                                                                                                                 >  4:304324/1‑148 (MQ=255)

                                 caAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                                >  4:798199/1‑148 (MQ=255)

                                  aaaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAg                                                                                                              >  1:651258/1‑149 (MQ=255)

                                   aaCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAg                                                                                                              >  3:698038/1‑148 (MQ=255)

                                    aCCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                                                                                                                              >  3:64697/1‑115 (MQ=255)

                                     ccGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTa                                                                                                            >  3:234067/1‑148 (MQ=255)

                                     ccGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCACGCGGGCCAcgaa                                                                                                                                             <  3:518294/115‑4 (MQ=255)

                                       gcgGAAGGGATCGACATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATg                                                                                                          >  3:806876/1‑148 (MQ=255)

                                           aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGt                                                                                                     <  1:753038/149‑1 (MQ=255)

                                             gggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGttt                                                                                                   <  2:1064656/149‑1 (MQ=255)

                                             gggATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtt                                                                                                    <  4:355455/148‑1 (MQ=255)

                                               gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCTTGGCCCGCGTGACACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGttttt                                                                                                 >  2:1037684/1‑149 (MQ=255)

                                               gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATg                                                                                                          >  3:946404/1‑140 (MQ=255)

                                               gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATg                                                                                                          <  4:946404/140‑1 (MQ=255)

                                                aTCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGttttt                                                                                                 <  4:454701/148‑1 (MQ=255)

                                                aTCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTc                                                                                                >  2:657553/1‑149 (MQ=255)

                                                aTCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTc                                                                                                <  1:606163/149‑1 (MQ=255)

                                                  cGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                                <  3:239045/131‑1 (MQ=255)

                                                  cGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                                >  4:239045/1‑131 (MQ=255)

                                                  cGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCa                                                                                               <  3:232553/148‑1 (MQ=255)

                                                    ccATCCAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGc                                                                                                                                                   >  1:12163/1‑94 (MQ=255)

                                                    ccATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAAc                                                                                             <  4:460883/148‑1 (MQ=255)

                                                    ccATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAAc                                                                                             <  3:349962/148‑1 (MQ=255)

                                                    ccATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACa                                                                                            <  2:1028301/149‑1 (MQ=255)

                                                    ccATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGACCAACCGTGGACCATCCGTCAGTATGCTCGTTTTTCAACa                                                                                            <  2:898924/149‑1 (MQ=255)

                                                      aTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGc                                                                                          >  2:714845/1‑149 (MQ=255)

                                                       tCAAGCGGCTGCATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGACCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGc                                                                                          >  4:665288/1‑148 (MQ=255)

                                                       tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCt                                                                                                        <  3:920944/134‑1 (MQ=255)

                                                       tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCCTCAGTATGCTGGTTTTTCAACAGCa                                                                                         >  2:531333/1‑149 (MQ=255)

                                                           gCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAag                                                                                     <  2:421779/149‑1 (MQ=255)

                                                              gCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAagag                                                                                   <  3:267538/148‑1 (MQ=255)

                                                                   ataCCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGTACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCaa                                                                              <  4:698038/148‑1 (MQ=255)

                                                                       cGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTc                                                                                                >  1:462463/1‑126 (MQ=255)

                                                                       cGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCtt                                                                         >  2:97254/1‑149 (MQ=255)

                                                                       cGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCtt                                                                         >  2:116359/1‑149 (MQ=255)

                                                                         aaGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCtttt                                                                       <  2:241983/149‑1 (MQ=255)

                                                                             cGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATc                                                                   <  2:170017/149‑1 (MQ=255)

                                                                              gATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATc                                                                   >  4:415574/1‑148 (MQ=255)

                                                                              gATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCg                                                                  >  2:434648/1‑149 (MQ=255)

                                                                                   cGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTa                                                             >  2:444183/1‑149 (MQ=255)

                                                                                     aTAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTaa                                                            >  3:413903/1‑148 (MQ=255)

                                                                                     aTAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTaa                                                            >  3:103004/1‑148 (MQ=255)

                                                                                     aTAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGACACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAc                                                           <  1:1037684/149‑1 (MQ=255)

                                                                                      tAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                          >  1:277416/1‑149 (MQ=255)

                                                                                       taTCTGGATGGCACAGGTAACTTTCCTGGTCTGCTGCCGTACGTTCGTGGACCTCGTGCCACAATGTATACCGCCCAACCGTGGACCATCCGTCAGTATTCTTGCTTTTCAGAACGAAAAGAGTCCAACGATTTTTATCGCCGTAAcc                                                          <  4:650936/147‑1 (MQ=255)

                                                                                       aaTCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                          <  4:413903/148‑1 (MQ=255)

                                                                                        aTCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCt                                                         <  3:368282/148‑1 (MQ=255)

                                                                                            ggAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGgccg                                                    >  2:121516/1‑149 (MQ=255)

                                                                                                  gACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCggg                                               <  3:404805/148‑1 (MQ=255)

                                                                                                   aCAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                  >  3:43895/1‑112 (MQ=255)

                                                                                                       gTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCaaaaa                                         <  1:500017/149‑1 (MQ=255)

                                                                                                       gTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCaaaa                                          <  4:667199/148‑1 (MQ=255)

                                                                                                        tACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAg                                        <  2:240229/149‑1 (MQ=255)

                                                                                                         aCCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAgg                                       <  2:33457/149‑1 (MQ=255)

                                                                                                          cccTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAgg                                       >  4:264805/1‑148 (MQ=255)

                                                                                                             ttCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCt                                    >  3:773415/1‑148 (MQ=255)

                                                                                                             ttCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGGATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCtt                                   >  1:279248/1‑149 (MQ=255)

                                                                                                              tCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATg                                                                                                          <  3:67387/77‑1 (MQ=255)

                                                                                                               ccTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCttt                                  <  3:858436/148‑1 (MQ=255)

                                                                                                                   gTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGtt                             <  1:936858/149‑1 (MQ=255)

                                                                                                                    tttGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGtt                             >  3:715615/1‑148 (MQ=255)

                                                                                                                       gccgccCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCg                          <  3:425012/148‑1 (MQ=255)

                                                                                                                       gccgccCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCg                          <  4:628260/148‑1 (MQ=255)

                                                                                                                            ccTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGa                     <  4:826452/148‑1 (MQ=255)

                                                                                                                            ccTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGAc                    <  2:365292/149‑1 (MQ=255)

                                                                                                                             cTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGAcc                   <  2:279248/149‑1 (MQ=255)

                                                                                                                                 gTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGc               >  1:845412/1‑149 (MQ=255)

                                                                                                                                 gTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGc               <  2:845412/149‑1 (MQ=255)

                                                                                                                                 gTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGc               >  2:747646/1‑149 (MQ=255)

                                                                                                                                   tagTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGcc              >  3:135391/3‑148 (MQ=255)

                                                                                                                                   tagTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGcc              >  3:62925/3‑148 (MQ=255)

                                                                                                                                   tCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGcc              >  3:157016/1‑148 (MQ=255)

                                                                                                                                    agTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGTCCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGcca             >  4:612996/2‑148 (MQ=255)

                                                                                                                                    agTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGcca             >  4:387578/2‑148 (MQ=255)

                                                                                                                                       ggCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGccaccca         >  2:574625/1‑149 (MQ=255)

                                                                                                                                         cacgcgGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGccacccac        <  4:62925/146‑1 (MQ=255)

                                                                                                                                           cgcgTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCg      >  3:434821/1‑148 (MQ=255)

                                                                                                                                              gTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGc  >  2:333636/1‑149 (MQ=255)

CGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGC  >  NC_000913/3060898‑3061188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`