BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,542	A→T	16.7%	I12I (ATT→ATA)	yieL ←	putative xylanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,542	0	A	T	16.7%	495.4 / 87.3	203	I12I (ATT→ATA)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base A (87/82); new base T (15/19); total (102/101)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 208 total aligned reads displayed.

GGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACT  >  NC_000913/3900443‑3900641

 gCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCGGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACAGCCGGCaa                                                                                                    >  2:3360458/1‑100 (MQ=255)

  cATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAAt                                                                                                   >  2:1516445/1‑100 (MQ=255)

   aTTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATa                                                                                                  >  1:104893/1‑100 (MQ=255)

   aTTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATa                                                                                                  >  1:3148976/1‑100 (MQ=255)

   aTTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATa                                                                                                  <  2:3370549/100‑1 (MQ=255)

    ttGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAg                                                                                                 <  2:430857/100‑1 (MQ=255)

      gACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGcc                                                                                               >  2:3519361/1‑100 (MQ=255)

         cTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGc                                                                                            <  1:3597713/100‑1 (MQ=255)

           gAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCgg                                                                                          <  1:3410341/100‑2 (MQ=255)

           gAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCgg                                                                                          <  1:1797537/100‑2 (MQ=255)

           gAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCgg                                                                                          <  2:3410683/100‑2 (MQ=255)

           gAGTCACATACTGTTGCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGACATACCGCTGGCTATTGCCAGCgg                                                                                          <  2:256609/100‑2 (MQ=255)

            aGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCggt                                                                                         <  1:319031/100‑3 (MQ=255)

            aGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCggt                                                                                         <  1:617435/100‑3 (MQ=255)

            aGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCggt                                                                                         <  2:4056714/100‑3 (MQ=255)

             gTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTa                                                                                        <  1:3012629/100‑1 (MQ=255)

             gTCACATACGGTTTCACCGGGATAGTCGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCGGTa                                                                                        <  2:3941122/100‑1 (MQ=255)

              tCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaa                                                                                       <  1:3640276/100‑1 (MQ=255)

               cacaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa                                                                                      <  1:886444/100‑1 (MQ=255)

                acaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAag                                                                                     <  1:1471994/100‑1 (MQ=255)

                 caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc                                                                                    <  1:1647973/100‑1 (MQ=255)

                 caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagc                                                                                    <  1:3338734/100‑1 (MQ=255)

                   tACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagcag                                                                                  <  1:1824698/100‑1 (MQ=255)

                   tACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAagcag                                                                                  <  1:4056122/100‑1 (MQ=255)

                     cTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAgcagca                                                                                >  1:805415/1‑100 (MQ=255)

                       gTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAt                                                                              >  2:326713/1‑100 (MQ=255)

                       gTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAt                                                                              >  1:357930/1‑100 (MQ=255)

                       gTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAt                                                                              >  2:1989869/1‑100 (MQ=255)

                        tttCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGATATGGACCACTGAGAAGAAATATCGATGGAAATAGCCAGCGTTAAAGAAGAAAtt                                                                             >  2:1001152/1‑100 (MQ=255)

                         ttCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                            >  2:2177276/1‑100 (MQ=255)

                         ttCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttt                                                                            >  1:1315953/1‑100 (MQ=255)

                          tCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAtttt                                                                           >  2:414267/1‑100 (MQ=255)

                           cACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                          <  1:849499/100‑1 (MQ=255)

                           cACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                          >  2:1602293/1‑100 (MQ=255)

                           cACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                          <  2:2938363/100‑1 (MQ=255)

                           cACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                          >  1:185388/1‑100 (MQ=255)

                           cACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                          >  2:369283/1‑100 (MQ=255)

                            aCCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTta                                                                         <  1:3179756/100‑1 (MQ=255)

                            aCCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTta                                                                         <  2:1794323/100‑1 (MQ=255)

                            aCCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTta                                                                         <  1:2396449/100‑1 (MQ=255)

                             ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTtat                                                                        >  1:1252404/1‑100 (MQ=255)

                             ccGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTtat                                                                        >  1:2097195/1‑100 (MQ=255)

                               gggATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTatat                                                                      >  2:2796218/1‑100 (MQ=255)

                                ggATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATAtt                                                                     >  1:409508/1‑100 (MQ=255)

                                  tctgTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCa                                                                   <  2:1929471/97‑1 (MQ=255)

                                     gTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCaaga                                                                >  1:1929474/1‑97 (MQ=255)

                                     gTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCg                                                                <  1:736451/100‑1 (MQ=255)

                                                gctggctgGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCtt                                                     <  2:769773/100‑1 (MQ=255)

                                                  tggctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTc                                                   <  1:659551/100‑1 (MQ=255)

                                                     ctggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCAtt                                                >  1:112141/1‑100 (MQ=255)

                                                      tggcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTa                                               <  1:2937539/100‑1 (MQ=255)

                                                        gcATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAg                                             <  2:409507/100‑1 (MQ=255)

                                                         cATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGt                                            <  2:300240/100‑1 (MQ=255)

                                                          atatCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGtt                                           >  2:2352521/1‑100 (MQ=255)

                                                           tatCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTc                                          <  1:4029919/100‑1 (MQ=255)

                                                             tCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCa                                        <  1:4172463/100‑1 (MQ=255)

                                                             tCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCa                                        <  1:2766242/100‑1 (MQ=255)

                                                              ccGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCAt                                       >  2:1246617/1‑100 (MQ=255)

                                                               cGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCAtt                                      <  1:2533231/100‑1 (MQ=255)

                                                                 ctgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTac                                    <  1:3233000/100‑1 (MQ=255)

                                                                  tgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacc                                   <  2:3093990/100‑1 (MQ=255)

                                                                  tgctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacc                                   <  2:2048166/100‑1 (MQ=255)

                                                                   gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacca                                  >  2:1906599/1‑100 (MQ=255)

                                                                   gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacca                                  >  1:3781086/1‑100 (MQ=255)

                                                                   gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacca                                  >  2:209720/1‑100 (MQ=255)

                                                                   gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacca                                  >  1:3488940/1‑100 (MQ=255)

                                                                   gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCAGTAAGTTCCATTacca                                  >  1:4129250/1‑100 (MQ=255)

                                                                     tATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTTAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTaccacc                                >  2:2068255/1‑100 (MQ=255)

                                                                      aTGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTAccacca                               >  1:2461958/1‑100 (MQ=255)

                                                                        ggCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGa                             >  1:3354779/1‑100 (MQ=255)

                                                                        ggCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGa                             >  2:1651168/1‑100 (MQ=255)

                                                                          cccACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGAtt                           >  1:3932175/1‑100 (MQ=255)

                                                                            cACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTc                         <  1:267485/100‑1 (MQ=255)

                                                                                gAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTTAGc                      >  2:762192/1‑99 (MQ=255)

                                                                                gAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGcc                     <  2:724740/100‑1 (MQ=255)

                                                                                gAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGcc                     <  1:2502668/100‑1 (MQ=255)

                                                                                  gctgctttaacGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTg                   >  1:223434/11‑100 (MQ=38)

                                                                                   cAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTgg                  >  1:2779933/1‑100 (MQ=255)

                                                                                    aGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTggg                 <  1:1956203/100‑1 (MQ=255)

                                                                                     gctttaacGCTGGCTATTGCCAGCGGTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGc                >  1:2250421/8‑100 (MQ=38)

                                                                                     gAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGc                >  1:2330837/1‑100 (MQ=255)

                                                                                      ctttaacGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATTGTAATATTCCTTTCGTTAAGTTCCATTAGCACCAGATTTCAGCCTGGGCa               >  2:1446444/7‑100 (MQ=255)

                                                                                      ctttaacGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCa               >  1:2697099/7‑100 (MQ=255)

                                                                                       tttaacGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCAc              >  1:2348709/6‑100 (MQ=255)

                                                                                       tttaacGCTGGCTATTGCCAGCGGTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCAc              >  1:12829/6‑100 (MQ=38)

                                                                                       aaTACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCAc              <  1:2801847/100‑1 (MQ=255)

                                                                                        aTACCGCTGGCAATAGCCAGCGTTAACGCTGTAATTCTTATCTTCATCGTCATATTCCTCTCATTAAGTCCCATCACCACCCAAATTCAGCCTAGGCAcc             <  2:3979835/100‑1 (MQ=255)

                                                                                         tACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCg            >  2:601809/1‑100 (MQ=255)

                                                                                         tACCGCTGGCAATAGCCAGCGCTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATAAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCg            <  2:2484732/100‑1 (MQ=255)

                                                                                          aacGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCAc              >  2:3840652/3‑97 (MQ=255)

                                                                                          aacGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGa           >  1:3115753/3‑100 (MQ=255)

                                                                                            cGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACa         >  2:2323283/1‑100 (MQ=255)

                                                                                             gCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATCTTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAg        >  1:3832736/1‑100 (MQ=255)

                                                                                             gCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAg        <  2:2673695/100‑1 (MQ=255)

                                                                                             gCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAg        >  1:1700658/1‑100 (MQ=255)

                                                                                               tGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCa      <  2:1291979/100‑1 (MQ=255)

                                                                                                ggCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCaa     >  1:957796/1‑100 (MQ=255)

                                                                                                  cAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAAc   <  1:1107246/100‑1 (MQ=255)

                                                                                                  cAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAAc   >  1:911533/1‑100 (MQ=255)

                                                                                                  cAATAGCCAGCGTTAAAGCAGAAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACGAGATTTCAGCCTGGGCACCGACAGCAAAc   <  2:3734517/100‑1 (MQ=255)

GGCATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACT  >  NC_000913/3900443‑3900641

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