BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,193:1	(C)_6→7	100%	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	95.6%	393.5 / 10.9	114	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (4/1); new base C (52/57); total (56/58)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 128 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                   <  2:894362/100‑1 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCTACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                 >  1:955907/1‑100 (MQ=255)

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                 >  1:1132192/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCGCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c                                                                                               >  1:1865309/1‑100 (MQ=255)

     gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccc                                                                                           <  1:259211/100‑1 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccccc                                                                                         <  1:1050733/100‑2 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccccc                                                                                         <  2:1736586/100‑2 (MQ=255)

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTTCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccccctg                                                                                       <  1:2214966/100‑4 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccccctgt                                                                                      >  2:938282/1‑96 (MQ=255)

          gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccccctgt                                                                                      >  1:1796069/1‑96 (MQ=255)

           aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c‑c‑ccccctgtt                                                                                     >  1:563280/1‑95 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTc                                                                                     >  1:752538/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTc                                                                                     >  1:442938/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTc                                                                                     <  2:2438723/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTg                                                                                   <  2:1411012/100‑1 (MQ=255)

               cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTgg                                                                                  >  1:292751/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGaa                                                                                >  2:2803861/1‑100 (MQ=255)

                 ccGCCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGaa                                                                                >  1:37508/1‑100 (MQ=255)

                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGaaa                                                                               <  2:2009500/100‑1 (MQ=255)

                  cGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGaaa                                                                               <  2:828854/100‑1 (MQ=255)

                    tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAAt                                                                             <  1:2994315/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGAc                                                                          <  2:889295/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGAc                                                                          <  2:2266603/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CGC‑CCCCTGTTCTGGAAAATGAc                                                                          >  1:416235/1‑100 (MQ=255)

                         tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTa                                                                        <  1:1713402/100‑1 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAg                                                                       >  2:2666204/1‑100 (MQ=255)

                          cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑C‑CCCCTGTTCTGGAAAATGACTAGc                                                                      >  2:816923/1‑100 (MQ=255)

                           acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGc                                                                      <  2:2163575/100‑1 (MQ=255)

                            cgtAGCTTCAAATACTTCCCACATAGTCTTCATATCCTGCTGGCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     >  2:2064357/3‑100 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  2:1711845/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGcc                                                                     <  1:1836044/100‑1 (MQ=255)

                              tAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTg                                                                   <  2:809792/100‑1 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                <  2:1157942/100‑1 (MQ=255)

                                 cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑CACCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2393204/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:709178/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:2682742/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:1478200/1‑100 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             <  2:336377/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:185077/1‑100 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:2893414/100‑1 (MQ=255)

                                     aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:1371693/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:1836351/100‑1 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           >  1:417849/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:1335666/100‑1 (MQ=255)

                                         tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                        <  1:1620458/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  2:2407778/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  2:2697209/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       <  2:955912/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  1:1222950/1‑100 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTATTTGGTGCAGTa                                                       <  2:2008588/100‑1 (MQ=255)

                                          cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGGTCTTGAAAAAGAATTGCCTGTTTgg                                                              >  1:1380395/1‑93 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAGAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGCTGCAGTATTg                                                    >  1:2014241/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:678487/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  1:31329/100‑1 (MQ=255)

                                                acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  2:1181253/1‑100 (MQ=255)

                                                acaAAGTCCGCGTATCCTCCTGCCATTGCGAAGGAGAAGACTATGGTGAAATCC‑C‑CCCCCGTTCTGGAAAATGACTAGCCTGTTTGGTCCAGTACTGc                                                   <  1:168867/98‑1 (MQ=255)

                                                 caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGG‑AAATGACTAGCCTGTTTGGTGCAGTATTGCtgtt                                               >  1:450255/1‑100 (MQ=255)

                                                  atagccggggTATCCTGCTGCCATTCCAATGGAGAAGACGATGTTGAAGTCC‑C‑CCCCAGTTCTGGAAAATAACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  1:672418/10‑98 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:2247785/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:2245469/1‑100 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:331427/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  2:2119361/1‑100 (MQ=255)

                                                       cTTCGTATCCTGCTGCCGTTGCAAAGCAGAAGACTATGTTGACATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:1388048/1‑100 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:2788304/100‑1 (MQ=255)

                                                        ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:461609/100‑1 (MQ=255)

                                                           gtctCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTCTTATGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:549442/4‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  1:731685/100‑1 (MQ=255)

                                                            tATCCTGCTGCCATTGCAAAGGAGAAGACTGTGATGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTCGTTGATTCCg                                     >  1:2016111/1‑100 (MQ=255)

                                                                ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                 >  1:1800917/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  2:1076237/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  1:2955518/1‑100 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  2:597163/1‑100 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:2884874/100‑1 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  1:3044986/100‑1 (MQ=255)

                                                                          tGCAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  1:2253256/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:1151426/100‑1 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:952416/1‑100 (MQ=255)

                                                                              aaGGAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   <  1:1046308/100‑1 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:1234972/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:2487786/1‑100 (MQ=255)

                                                                                  agaagaCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGa               >  2:2965445/1‑100 (MQ=255)

                                                                                   gaagaCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                   >  1:2412658/1‑95 (MQ=255)

                                                                                     agaCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTATCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  1:257001/100‑1 (MQ=255)

                                                                                     agaCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:2280910/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:893104/100‑1 (MQ=255)

                                                                                       aCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:1447250/100‑1 (MQ=255)

                                                                                       aCTATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:2682758/100‑1 (MQ=255)

                                                                                         tATGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        >  1:197352/1‑100 (MQ=255)

                                                                                          aTGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       >  2:2023879/1‑100 (MQ=255)

                                                                                           tGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  1:236398/100‑1 (MQ=255)

                                                                                           tGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  1:2803844/100‑1 (MQ=255)

                                                                                           tGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  1:1414763/100‑1 (MQ=255)

                                                                                           tGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:133900/1‑100 (MQ=255)

                                                                                           tGTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:133102/1‑100 (MQ=255)

                                                                                            gTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:2463007/1‑100 (MQ=255)

                                                                                            gTTGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2160351/1‑100 (MQ=255)

                                                                                             ttGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:1857259/1‑100 (MQ=255)

                                                                                              tGAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   <  2:2404332/100‑1 (MQ=255)

                                                                                               gAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:253970/100‑1 (MQ=255)

                                                                                               gAAATCC‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:259319/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑C‑C‑CCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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