BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)_6→7	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.1%	573.8 / 20.9	164	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base . (4/4); new base C (74/82); total (78/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                     >  1:3079803/1‑100 (MQ=255)

gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                     >  1:3938727/1‑100 (MQ=255)

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

 ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATCTTGAAAT‑‑c                                                                                                  >  2:3627302/1‑100 (MQ=255)

  cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c                                                                                                >  2:1342748/1‑100 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc                                                                                             <  2:792712/100‑1 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc                                                                                             <  2:2043995/100‑1 (MQ=255)

    cGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGGTGAAAT‑‑c‑c‑cc                                                                                             >  2:2738635/1‑100 (MQ=255)

      cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc‑cc                                                                                          >  1:2177720/1‑100 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc‑ccc                                                                                         <  2:2438245/100‑2 (MQ=255)

       ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc‑ccc                                                                                         <  2:511801/100‑2 (MQ=255)

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc‑ccctg                                                                                       >  1:3302744/1‑97 (MQ=255)

         cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑c‑c‑cc‑ccctg                                                                                       >  2:2901783/1‑97 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTc                                                                                     >  1:2486421/1‑100 (MQ=255)

            cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTc                                                                                     <  1:1415748/100‑1 (MQ=255)

            cTTCACCGTCACTTCACATAGCGTCAAATTTCTCCCATATAGTCTTCCTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑C‑CACC‑CCTGTTc                                                                                     <  2:2817398/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTg                                                                                   >  2:2418109/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTg                                                                                   <  2:3155642/100‑1 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTg                                                                                   >  2:618715/1‑100 (MQ=255)

              tCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTg                                                                                   <  1:3723507/100‑1 (MQ=255)

                aCCTTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTg                                                                                   <  2:1239666/98‑1 (MQ=255)

                aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGa                                                                                 >  1:1292295/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGaa                                                                                <  1:3337938/100‑1 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGaa                                                                                >  1:3609480/1‑100 (MQ=255)

                 ccGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGaa                                                                                >  2:4103895/1‑100 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGAc                                                                          <  1:1581055/100‑1 (MQ=255)

                       cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGAc                                                                          <  2:3883396/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGcc                                                                     <  2:1916066/100‑1 (MQ=255)

                            caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGcc                                                                     >  2:2685511/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:552102/1‑100 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGtt                                                                 <  2:3087840/100‑1 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGt                                                                  >  1:420936/1‑99 (MQ=255)

                                gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCCGTTCTGGAAAATGACAAGCCTGt                                                                  >  1:968233/1‑99 (MQ=255)

                                 cTTCAAATTCTTCCCACATACTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGttt                                                                >  2:2553662/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               >  1:2093552/1‑100 (MQ=255)

                                  ttCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTg                                                               <  2:1196737/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:1521946/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  2:1128416/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  1:1618530/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              >  2:2219308/1‑100 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  1:1323118/100‑1 (MQ=255)

                                   tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTgg                                                              <  2:3764822/100‑1 (MQ=255)

                                    cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:732772/1‑100 (MQ=255)

                                      aaTTCTTCCCACATAGTCTCCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                           <  2:397920/100‑1 (MQ=255)

                                       aTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                          >  1:3598359/1‑100 (MQ=255)

                                        ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:2676546/100‑1 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  2:1742887/1‑100 (MQ=255)

                                             cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    >  1:3696507/1‑100 (MQ=255)

                                              tcacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  1:1947252/99‑1 (MQ=255)

                                              ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1089041/100‑1 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGca                                                  >  2:1947238/1‑99 (MQ=255)

                                               cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCt                                                  >  2:3715905/1‑100 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2412343/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  2:2789745/100‑1 (MQ=255)

                                                   tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2885779/100‑1 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  2:295409/100‑1 (MQ=255)

                                                    aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑C‑C‑CCACCTGTTCTGGATTATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  2:1391330/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:3463568/100‑1 (MQ=255)

                                                     gTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  2:2289341/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGGCTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:2292025/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            >  2:3502828/1‑99 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           >  1:1114432/1‑100 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  2:1830652/100‑1 (MQ=255)

                                                      tCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                           <  1:537867/100‑1 (MQ=255)

                                                       cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  2:1564109/1‑100 (MQ=255)

                                                          cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  1:501369/100‑1 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:1360272/1‑100 (MQ=255)

                                                           gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:3079770/100‑1 (MQ=255)

                                                             aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  2:3805719/1‑100 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:3416104/100‑1 (MQ=255)

                                                              tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   <  1:1755220/100‑1 (MQ=255)

                                                               ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  1:1062458/1‑100 (MQ=255)

                                                                  gctgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                               <  1:2991245/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  1:3214723/100‑1 (MQ=255)

                                                                   ctgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:2879150/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              <  2:780504/100‑1 (MQ=255)

                                                                    tgcCATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  1:2926974/100‑1 (MQ=255)

                                                                      ccATTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           <  1:1305389/100‑1 (MQ=255)

                                                                        aTTGCAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                         >  1:1307439/1‑100 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:88471/100‑1 (MQ=255)

                                                                            cAAAGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     >  2:3846477/1‑100 (MQ=255)

                                                                             aaaGGAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    >  1:1221211/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  2:413124/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAAT‑‑CTC‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 >  1:1148972/1‑100 (MQ=255)

                                                                                ggAGAAGACTATGTTGAAAT‑‑CAC‑CC‑CCTGTTCTGGAAAACGACTAGCCCGGTTGGTGCAGCATTGCTGATGTTGATTCCGATAATGATGACTCggaacgt                 >  2:3276304/1‑95 (MQ=255)

                                                                                   gaagaCTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  2:1516253/100‑1 (MQ=255)

                                                                                   gaagaCTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              >  1:4027125/1‑100 (MQ=255)

                                                                                     agaCTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:2629588/1‑100 (MQ=255)

                                                                                      gaCTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  1:2240796/1‑100 (MQ=255)

                                                                                       aCTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  2:2107676/100‑1 (MQ=255)

                                                                                        cTATGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:2974066/100‑1 (MQ=255)

                                                                                           tGTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      <  2:2486396/100‑1 (MQ=255)

                                                                                            gTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:2999965/1‑100 (MQ=255)

                                                                                            gTTGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:79158/1‑100 (MQ=255)

                                                                                             ttGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:3062689/1‑100 (MQ=255)

                                                                                             ttGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:3364926/1‑100 (MQ=255)

                                                                                              tGAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg   >  1:694629/1‑100 (MQ=255)

                                                                                               gAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTgn  <  1:2011623/100‑2 (MQ=255)

                                                                                               gAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  2:3598320/100‑1 (MQ=255)

                                                                                               gAAATC‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  1:207505/1‑100 (MQ=255)

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑‑C‑C‑CC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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