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.2%	420.9 / 13.3	126	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (2/4); new base C (60/60); total (62/64)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638089‑4638280

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

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

  tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc                                                                                             <  1:1660088/100‑1 (MQ=255)

    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cc                                                                                          <  2:749747/100‑1 (MQ=255)

    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑c                                                                                           <  2:2140176/99‑1 (MQ=255)

     cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑ccc                                                                                         >  1:2043679/1‑100 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccc                                                                                        >  1:1932170/1‑99 (MQ=255)

      ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCACACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGATGAAA‑T‑ccc‑cccc                                                                                        <  2:3009221/100‑2 (MQ=255)

         gACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCC‑CCCTGtt                                                                                     <  2:1554709/100‑1 (MQ=255)

          aCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑ccc‑cccctgtt                                                                                    >  1:1583399/1‑95 (MQ=255)

           cctCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                    >  2:621900/3‑100 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                    >  2:405822/1‑100 (MQ=255)

           cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTc                                                                                    >  2:1958775/1‑100 (MQ=255)

            ttCACCGTCACTTCGCATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCt                                                                                   >  1:3057764/1‑100 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                  <  2:910074/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                  <  1:2244055/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                  <  1:616998/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                  >  1:2107856/1‑100 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                  <  1:424675/100‑1 (MQ=255)

             tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTg                                                                                  >  2:232867/1‑100 (MQ=255)

              cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTgg                                                                                 <  2:1139230/100‑1 (MQ=255)

                ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGAATATGTTGAAA‑TCCCC‑CCCTGTTCTGGaa                                                                               >  2:686794/1‑100 (MQ=255)

                   tCACTTCATATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGaaaa                                                                             <  2:2652524/99‑1 (MQ=255)

                   tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAAt                                                                            <  2:892520/100‑1 (MQ=255)

                    cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATg                                                                           >  2:3150242/1‑100 (MQ=255)

                      cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGAc                                                                         <  2:448169/100‑1 (MQ=255)

                      cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGAc                                                                         <  1:2854047/100‑1 (MQ=255)

                       ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACt                                                                        <  2:2964262/100‑1 (MQ=255)

                         cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAg                                                                      >  1:1529730/1‑100 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAg                                                                      >  1:84580/1‑98 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGcc                                                                    >  1:60717/1‑100 (MQ=255)

                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGcc                                                                    <  1:2467551/100‑1 (MQ=255)

                               gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGtt                                                                >  1:2478076/1‑100 (MQ=255)

                                cTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGttt                                                               <  1:2305861/100‑1 (MQ=255)

                                 ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTg                                                              >  2:2774418/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             >  1:2689658/1‑100 (MQ=255)

                                  tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                             >  1:1352584/1‑100 (MQ=255)

                                    aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                           >  2:3165534/1‑100 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                          <  1:2608091/100‑1 (MQ=255)

                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                          <  1:453303/100‑1 (MQ=255)

                                       ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                        <  2:416570/100‑1 (MQ=255)

                                        tcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                       <  2:2252615/100‑1 (MQ=255)

                                         cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGACGACTATGTTGAAA‑T‑CCCACCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGt                                                       >  1:2249580/1‑99 (MQ=255)

                                         cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                      >  2:1427255/1‑100 (MQ=255)

                                           tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                     >  2:2533620/1‑99 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                   <  1:1526575/99‑1 (MQ=255)

                                             ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                  <  2:1189740/100‑1 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                   >  2:2694548/1‑98 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  2:2611838/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  2:1626741/1‑100 (MQ=255)

                                              cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACCATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                 >  1:1294242/1‑100 (MQ=255)

                                                  tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                             <  2:586169/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  2:921091/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  2:1526436/100‑1 (MQ=255)

                                                    gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCGGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                           <  2:1042191/100‑1 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                         >  1:564612/1‑100 (MQ=255)

                                                      cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                         >  1:2337231/1‑100 (MQ=255)

                                                       ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        <  1:68028/100‑1 (MQ=255)

                                                        tCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                        >  1:1184707/1‑99 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     <  2:1231473/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     <  2:2751551/100‑1 (MQ=255)

                                                          gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                     >  1:777714/1‑100 (MQ=255)

                                                           tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGGTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                    >  2:929366/1‑100 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                  <  2:48508/99‑1 (MQ=255)

                                                              ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                 >  1:1043269/1‑100 (MQ=255)

                                                              ccTGCCGCCATTGCAAAAGGGAAGACTATGTTGTAC‑TCCCC‑CCCTGTTCTGGAAATTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                  <  2:2779353/99‑1 (MQ=255)

                                                               ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                >  1:2160803/1‑100 (MQ=255)

                                                               ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATaa                                >  2:1562387/1‑100 (MQ=255)

                                                                tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                               >  1:2317152/1‑100 (MQ=255)

                                                                tgctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                               >  1:1259328/1‑100 (MQ=255)

                                                                  cttccgatctaaaGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             <  2:3101965/90‑1 (MQ=255)

                                                                  ctgcCATTGCAAAGGAGAAGACTATGTTGAAATT‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             <  1:2970580/100‑1 (MQ=255)

                                                                  ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             >  2:616401/1‑100 (MQ=255)

                                                                  ctgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             <  1:2989079/100‑1 (MQ=255)

                                                                   tgcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                             <  1:3023041/99‑1 (MQ=255)

                                                                    gcCATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                           <  2:2799969/100‑1 (MQ=255)

                                                                      cATTGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                         >  1:1557651/1‑100 (MQ=255)

                                                                       aTTGCAAAGGAGAAGACTATGTTGAAA‑GCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                        <  1:2718390/100‑1 (MQ=255)

                                                                        ttGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                        >  2:3158447/1‑99 (MQ=255)

                                                                        ttGCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                       <  2:1435422/100‑1 (MQ=255)

                                                                          gCAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCgg                     >  2:1153000/1‑100 (MQ=255)

                                                                           cAAAGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                    >  2:672163/1‑100 (MQ=255)

                                                                            aaaGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                   <  1:108341/100‑1 (MQ=255)

                                                                             aaGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                  <  1:2887776/100‑1 (MQ=255)

                                                                             aaGGAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                  <  2:1184701/100‑1 (MQ=255)

                                                                               ggAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                >  1:2171050/1‑100 (MQ=255)

                                                                               ggAGAAGACTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                >  2:2061658/1‑100 (MQ=255)

                                                                                  gaagaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt             >  2:2251418/1‑100 (MQ=255)

                                                                                    agaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           >  2:1025315/1‑100 (MQ=255)

                                                                                    agaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           >  2:1984233/1‑100 (MQ=255)

                                                                                    agaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg           <  2:2007335/100‑1 (MQ=255)

                                                                                     gaCTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt          <  1:2019012/100‑1 (MQ=255)

                                                                                      aCTATGTTGAAT‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTCTTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         <  1:1871111/100‑1 (MQ=255)

                                                                                       cTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc         <  1:1092962/99‑1 (MQ=255)

                                                                                       cTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg        >  1:78267/1‑100 (MQ=255)

                                                                                       cTATGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg        >  2:2328408/1‑100 (MQ=255)

                                                                                         aTGTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGaa      <  2:2993283/100‑1 (MQ=255)

                                                                                           gTTGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg    <  2:2293856/100‑1 (MQ=255)

                                                                                             tGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCGTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAGGCTGATTCGGCAGGTGATTGTCGAACGTg  >  2:1952256/1‑100 (MQ=255)

                                                                                             tGAAA‑TCCCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg  >  1:294361/1‑100 (MQ=255)

CCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAA‑T‑CCC‑CCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTG  >  NC_000913/4638089‑4638280

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