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	94.5%	297.3 / 11.6	91	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (1/4); new base C (52/33); total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c                                                                                                >  1:2287681/1‑100 (MQ=255)

  tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cc                                                                                             <  1:2369549/100‑1 (MQ=255)

  tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cc                                                                                             <  2:1521933/100‑1 (MQ=255)

   cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑ccc                                                                                            <  2:2316855/100‑1 (MQ=255)

    gCTTCGACTTCACCGTCACTTCACATAGTTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cccc                                                                                           <  1:1018004/100‑1 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑ccccc                                                                                          <  1:1446044/100‑1 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cccc                                                                                           >  1:2608246/1‑99 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cccccc                                                                                         <  2:1752903/100‑2 (MQ=255)

       tCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cccccct                                                                                        <  2:378369/100‑3 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑c‑cccccctgt                                                                                      <  2:1014576/100‑5 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATCGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAATACTATGTTGAAA‑T‑‑c‑cccccctgtt                                                                                     >  2:35739/1‑95 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTc                                                                                     >  2:2654347/1‑100 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTc                                                                                     >  2:321912/1‑100 (MQ=255)

            ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCt                                                                                    <  2:880382/100‑1 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTgg                                                                                  <  1:380235/100‑1 (MQ=255)

               aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGa                                                                                 >  2:2076192/1‑100 (MQ=255)

                   tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAAt                                                                             <  2:2110209/100‑1 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGa                                                                           >  1:2066890/1‑100 (MQ=255)

                     aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGa                                                                           <  2:477227/100‑1 (MQ=255)

                      cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGAc                                                                          >  1:1610928/1‑100 (MQ=255)

                       ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACt                                                                         >  1:1932260/1‑100 (MQ=255)

                        tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGAc                                                                          <  1:2470996/98‑1 (MQ=255)

                        tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTa                                                                        <  1:151918/100‑1 (MQ=255)

                         cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAg                                                                       >  2:2344148/1‑100 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGcc                                                                     >  2:1425619/1‑100 (MQ=255)

                            aTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCt                                                                    >  2:1451532/1‑100 (MQ=255)

                              aGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGTAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGt                                                                  <  2:81341/100‑1 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGGAGACTATGTAGAAA‑TG‑C‑CCCCCTGATCTGGAAAATGACTAGCCTGtt                                                                 >  2:1787312/1‑100 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:1772223/1‑100 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  2:641526/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:1483534/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                <  2:1827520/100‑1 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:1313528/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  1:2058270/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTacccaccaca                                                               >  1:799927/1‑91 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:2223375/1‑100 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:1109555/100‑1 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  2:909208/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  2:400313/100‑1 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:335419/1‑100 (MQ=255)

                                   naaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  1:2635795/2‑100 (MQ=255)

                                   cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:1344603/1‑100 (MQ=255)

                                    agaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            >  1:1797165/3‑99 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:2225152/100‑1 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑C‑CCCCCTGTTCTGGAAAATGACCACCCTGTTTGGTGc                                                           <  1:2566583/100‑1 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  1:1658755/100‑1 (MQ=255)

                                       ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:1483352/100‑1 (MQ=255)

                                            cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTAGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  2:1414994/1‑100 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1804835/100‑1 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:2025027/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:402922/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  1:190266/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:70747/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:1540634/1‑100 (MQ=255)

                                               acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 <  2:2642452/100‑1 (MQ=255)

                                               acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  2:541874/1‑100 (MQ=255)

                                                caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                <  1:2287308/100‑1 (MQ=255)

                                                  tAGTCTTCGTATCCTGCTGCTATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:1810720/100‑1 (MQ=255)

                                                  tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2472521/100‑1 (MQ=255)

                                                  tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              >  2:1271456/1‑100 (MQ=255)

                                                   aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  1:293679/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  1:1271468/100‑1 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttgn                                          <  2:2635764/100‑2 (MQ=255)

                                                       ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         >  1:623893/1‑100 (MQ=255)

                                                             tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  1:1963381/1‑100 (MQ=255)

                                                             tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATAACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  1:1441892/1‑100 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGACA‑TC‑C‑CCCCCTGTTCTGGAAAAGGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:632805/1‑99 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:939597/1‑100 (MQ=255)

                                                                tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                <  2:1262790/100‑1 (MQ=255)

                                                                tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                <  2:2137165/100‑1 (MQ=255)

                                                                 gctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               >  2:428185/1‑100 (MQ=255)

                                                                  ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  1:2009273/100‑1 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             >  2:2264858/1‑100 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:1346339/1‑100 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑CTCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  1:979801/100‑1 (MQ=255)

                                                                       aTTGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         <  2:1205473/100‑1 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAACG‑‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTCGATTCCGATAATGCTGATTc                        <  2:1613816/99‑1 (MQ=255)

                                                                         tGCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        >  1:1101163/1‑99 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAA‑TCCC‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       >  1:477327/1‑100 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:131669/1‑100 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                      >  1:218904/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:1502501/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:1687776/1‑100 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  2:2066868/100‑1 (MQ=255)

                                                                             aaGGAGAAGACTATGTTGAAA‑T‑‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAgg                  <  1:2145163/100‑1 (MQ=255)

                                                                              aGGAGAAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTTCAGTATTGCAGTTATTGATTCAGATAATGCTAATTCGGCa                    >  1:1268545/1‑98 (MQ=255)

                                                                                gagaAGACTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTg                >  2:2160415/1‑100 (MQ=255)

                                                                                 agaagaCTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGt                 >  2:395705/1‑98 (MQ=255)

                                                                                    agaCTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAtt             >  2:1319712/1‑99 (MQ=255)

                                                                                      aCTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:641530/100‑1 (MQ=255)

                                                                                      aCTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:454116/100‑1 (MQ=255)

                                                                                       cTATGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:1344616/100‑1 (MQ=255)

                                                                                         aTGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa       <  1:604778/100‑1 (MQ=255)

                                                                                          tGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:1277003/1‑100 (MQ=255)

                                                                                          tGTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  1:1277015/100‑1 (MQ=255)

                                                                                           gTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:329896/1‑100 (MQ=255)

                                                                                           gTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:682333/1‑100 (MQ=255)

                                                                                           gTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:1990480/1‑100 (MQ=255)

                                                                                           gTTGAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  1:1772240/100‑1 (MQ=255)

                                                                                              gAAA‑TC‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  1:640546/100‑1 (MQ=255)

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑‑C‑CCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638089‑4638281

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: `‑`