BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,900,545	A→T	16.8%	A11A (GCT→GCA)	yieL ←	putative xylanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,900,545	0	A	T	16.8%	489.9 / 91.4	203	A11A (GCT→GCA)	yieL	putative xylanase
Reads supporting (aligned to +/- strand): ref base A (84/85); new base T (15/19); total (99/104)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-03

Only 100 of 208 total aligned reads displayed.

ATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGT  >  NC_000913/3900446‑3900644

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

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

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

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

    aCCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCa                                                                                                 <  2:2521648/100‑1 (MQ=255)

      cTGATTCACATCCTGGTTCACCGGGCTAGTGGGTGCCGGGCTGGCCGGCATATCCGCTGCTATCGACCACTGAGCAGAAATACCGCTGGCAATAGc                                                                                                   <  2:3897009/96‑1 (MQ=255)

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

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

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

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

        gAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGt                                                                                             >  1:3797097/1‑100 (MQ=255)

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

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

          gTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCTATTGCCAGCGGTa                                                                                           <  2:3730546/100‑1 (MQ=255)

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

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

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

            cacaTACTGTTTCACCGGGAGAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTaaa                                                                                         <  2:2253791/100‑1 (MQ=255)

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

             acaTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAag                                                                                        >  2:1281271/1‑100 (MQ=255)

              caTACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACTGCTGGCAATAGCCAGCGTTaa                                                                                          >  2:3349001/1‑97 (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)

                  cTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAgcagca                                                                                   <  2:625878/100‑1 (MQ=255)

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

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

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

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

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

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

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

                        cACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAAttttt                                                                             >  2:2352262/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)

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

                          ccGGGATAGTGGGTGCCGGGCTGACTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTtat                                                                           >  2:851514/1‑100 (MQ=255)

                             ggATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTatat                                                                         >  2:3955987/1‑99 (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)

                                      gTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCGCTGGGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAat                                                               >  2:1394481/1‑100 (MQ=255)

                                        gCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAatat                                                             >  2:294173/1‑100 (MQ=255)

                                          cGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTc                                                           <  2:1119585/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)

                                                        tatCCGCTGCTATGGCCCAGTGAGCAGAAATATCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATCTTAATCGTAATATTGCTTTCATTAAGTTc                                             >  2:4106707/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)

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

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

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

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

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

                                                                gctATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTacca                                     >  2:1906599/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)

                                                                           cTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAg                          <  2:1493359/100‑1 (MQ=255)

                                                                           cTGAGCAGAAATACCGCTAGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAg                          <  2:2405355/100‑1 (MQ=255)

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

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

                                                                              aGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCt                       >  2:7641/1‑100 (MQ=255)

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

                                                                                cAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTGATTAAGTTCCATTACCACCAGATTTCAGCCTgg                     <  2:2760171/100‑1 (MQ=255)

                                                                                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)

                                                                                   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)

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

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

                                                                                       aCCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGa              >  2:2847200/1‑100 (MQ=255)

                                                                                        ccGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGAc             <  2:2650205/100‑1 (MQ=255)

                                                                                        acGCTGGCTATTGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTACATTACCACCAGATTTCAGCCTGGGCACCGAc             >  2:4079353/2‑100 (MQ=255)

                                                                                         cGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACa            <  2:3211549/100‑1 (MQ=255)

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

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

                                                                                           cTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGc          >  2:3028897/1‑100 (MQ=255)

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

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

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

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

                                                                                                aaTAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACt     <  2:500906/100‑1 (MQ=255)

ATTGACCTGAGTCACATACTGTTTCACCGGGATAGTGGGTGCCGGGCTGGCTGGCATATCCGCTGCTATGGCCCACTGAGCAGAAATACCGCTGGCAATAGCCAGCGTTAAAGCAGCAATTTTTATATTCATCGTAATATTCCTTTCATTAAGTTCCATTACCACCAGATTTCAGCCTGGGCACCGACAGCAAACTGGT  >  NC_000913/3900446‑3900644

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