BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,853,574	T→C	R8R (AGA→AGG)	yidG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,853,574	0	T	C	99.8%	1325.0 / ‑5.7	422	R8R (AGA→AGG)	yidG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (159/262); total (159/263)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 446 total aligned reads displayed.

AGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGT  >  NC_000913/3853484‑3853663

aGCCGACGCCATCAACGCGCCGTAGCCCAGCATGGTCCGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCt                                                                                         <  2:1965254/93‑1 (MQ=255)

  ccAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTgg                                                                                       >  2:2949110/1‑93 (MQ=255)

   cAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCCATGCGCCTGGc                                                                                      >  2:4540060/1‑93 (MQ=255)

   cAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGc                                                                                      <  1:1750828/93‑1 (MQ=255)

    aaagCCATCAACGCGACGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCt                                                                                     >  2:5526727/4‑93 (MQ=255)

    aaCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCt                                                                                     >  1:1662353/1‑93 (MQ=255)

     aCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGGCGGCAATGCGCCTGGCtt                                                                                    >  2:6543640/1‑93 (MQ=255)

      cGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCttt                                                                                   <  2:4767177/93‑1 (MQ=255)

      cGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCttt                                                                                   <  1:7830670/93‑1 (MQ=255)

      cGCCATCAACGCGCCCTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCttt                                                                                   <  2:7358995/93‑1 (MQ=255)

       gCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTg                                                                                  <  1:6366359/93‑1 (MQ=255)

        ccATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGc                                                                                 >  1:796033/1‑93 (MQ=255)

           tCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGc                                                                              >  2:581373/1‑93 (MQ=255)

            cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCt                                                                             <  1:7035634/93‑1 (MQ=255)

             aaCGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTa                                                                            <  1:5129347/93‑1 (MQ=255)

             aaCGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTa                                                                            <  1:575369/93‑1 (MQ=255)

               cgcgCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATc                                                                          <  1:7180362/93‑1 (MQ=255)

               cgcgCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATc                                                                          <  2:876144/93‑1 (MQ=255)

                gcgcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATcc                                                                         >  1:2087054/1‑93 (MQ=255)

                   ccGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGc                                                                      >  2:7552268/1‑93 (MQ=255)

                    cGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCa                                                                     >  2:5702369/1‑93 (MQ=255)

                      tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATa                                                                   >  2:5309/1‑93 (MQ=255)

                        gCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATaca                                                                 >  2:6988430/1‑93 (MQ=255)

                        gCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATaca                                                                 >  1:529481/1‑93 (MQ=255)

                          ccAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATacaac                                                               >  1:437974/1‑93 (MQ=255)

                            aGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAacac                                                             >  1:3272380/1‑93 (MQ=255)

                            aGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAacac                                                             <  1:756270/93‑1 (MQ=255)

                             gCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACAcc                                                            >  1:6915403/1‑93 (MQ=255)

                              cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCa                                                           <  2:4271442/93‑1 (MQ=255)

                              cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCa                                                           <  1:2666146/93‑1 (MQ=255)

                              cATGGGACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGGCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCa                                                           <  1:2007988/93‑1 (MQ=255)

                               aTGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAg                                                          <  2:4036579/93‑1 (MQ=255)

                               aTGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAg                                                          <  1:2980223/93‑1 (MQ=255)

                               aTGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAg                                                          >  2:3623611/1‑93 (MQ=255)

                                tGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCTCCTGGCTTTGGGGCTATCTGTCATCCAACACCAgc                                                         >  1:5841869/1‑92 (MQ=255)

                                 ggTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTc                                                        <  2:2365235/93‑1 (MQ=255)

                                  gTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTcc                                                       >  1:4120140/1‑93 (MQ=255)

                                   tACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTccc                                                      <  2:251572/93‑1 (MQ=255)

                                      gAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATa                                                   <  1:3528793/93‑1 (MQ=255)

                                      gAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATa                                                   >  1:5510365/1‑93 (MQ=255)

                                       aaaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                                                  >  2:3487278/1‑93 (MQ=255)

                                       aaaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                                                  <  2:7115005/93‑1 (MQ=255)

                                       aaaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                                                  <  2:1926381/93‑1 (MQ=255)

                                       aaaCCACGCCAGTGATGTGCGCTCCGGCTCCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                                                  <  1:590702/93‑1 (MQ=255)

                                        aaCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAAc                                                 >  2:5947113/1‑93 (MQ=255)

                                         aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                                                <  2:4987509/93‑1 (MQ=255)

                                         aCCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                                                >  1:4774463/1‑93 (MQ=255)

                                          ccACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGa                                               <  2:1635994/93‑1 (MQ=255)

                                           cACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                                              <  2:6140406/93‑1 (MQ=255)

                                           cACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                                              <  2:262977/93‑1 (MQ=255)

                                           cACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                                              <  2:1254522/93‑1 (MQ=255)

                                           cACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                                              <  2:3899080/93‑1 (MQ=255)

                                           cACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                                              <  2:1443597/93‑1 (MQ=255)

                                             cGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCa                                            <  2:330948/93‑1 (MQ=255)

                                               ccAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCAcc                                          >  2:5602394/1‑93 (MQ=255)

                                                 aGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGc                                        <  2:7235245/93‑1 (MQ=255)

                                                  gTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGcc                                       >  2:129166/1‑93 (MQ=255)

                                                  gTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGCCATACAACACCAGTCCCATAACGATCACCGcc                                       >  2:5376826/1‑93 (MQ=255)

                                                    gATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc                                     <  2:1680048/93‑1 (MQ=255)

                                                         gcgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCa                                <  1:6249013/93‑1 (MQ=255)

                                                         gcgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGCCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCa                                <  2:3739350/93‑1 (MQ=255)

                                                          cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt                               <  1:1189284/93‑1 (MQ=255)

                                                           gcTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc                              <  1:7291046/93‑1 (MQ=255)

                                                           gcTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc                              >  1:4616183/1‑93 (MQ=255)

                                                            cTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCa                             >  1:669682/1‑93 (MQ=255)

                                                            cTCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCa                             >  1:3688143/1‑93 (MQ=255)

                                                             tCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa                            <  2:5280216/93‑1 (MQ=255)

                                                             tCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa                            <  2:6201964/93‑1 (MQ=255)

                                                             tCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa                            <  1:4340974/93‑1 (MQ=255)

                                                             tCCGGCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa                            >  2:1086754/1‑93 (MQ=255)

                                                                cgCTGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAatat                         <  2:2763878/92‑1 (MQ=255)

                                                                  gtGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTa                       <  1:1266633/92‑1 (MQ=255)

                                                                   tGCCGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTag                      >  2:7839771/1‑92 (MQ=255)

                                                                   tGCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTaa                      <  2:2270194/93‑1 (MQ=255)

                                                                    gCAGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAg                     <  2:7583864/93‑1 (MQ=255)

                                                                      aGTCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCt                   >  2:3432297/1‑93 (MQ=255)

                                                                        tCCTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTaa                 >  1:6487450/1‑93 (MQ=255)

                                                                         ccTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                <  2:63437/93‑1 (MQ=255)

                                                                         ccTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAAt                <  1:3979782/93‑1 (MQ=255)

                                                                          cTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATg               >  2:6737844/1‑93 (MQ=255)

                                                                          cTGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATg               >  1:3597966/1‑93 (MQ=255)

                                                                           tGGTTCGGCAATGCGCCTGGCGTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGa              <  2:2158840/93‑1 (MQ=255)

                                                                           tGGGTCGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGa              <  1:7748013/93‑1 (MQ=255)

                                                                            gggTCGGCAATGCGCCTGGCTTTGCGGCTATCCGTCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt             <  2:4404632/93‑1 (MQ=255)

                                                                              gTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATta           <  1:2290210/93‑1 (MQ=255)

                                                                                cGGCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaat         <  1:1484064/93‑1 (MQ=255)

                                                                                 ggCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata        <  1:3144750/93‑1 (MQ=255)

                                                                                 ggCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaata        >  1:2136213/1‑93 (MQ=255)

                                                                                  gCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaataa       <  2:5159114/93‑1 (MQ=255)

                                                                                  gCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaataa       <  1:6131577/93‑1 (MQ=255)

                                                                                  gCAATGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaataa       <  1:2434216/93‑1 (MQ=255)

                                                                                    aaTGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGc     >  1:1099633/1‑93 (MQ=255)

                                                                                    aaTGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGc     <  1:392922/93‑1 (MQ=255)

                                                                                    aaTGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGc     >  1:1014243/1‑93 (MQ=255)

                                                                                     aTGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCt    <  1:4147325/93‑1 (MQ=255)

                                                                                     aTGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCt    <  1:4888199/93‑1 (MQ=255)

                                                                                     aTGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCt    <  1:5378070/93‑1 (MQ=255)

                                                                                      tGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTg   >  1:3057049/1‑93 (MQ=255)

                                                                                      tGCGCCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTg   <  2:419242/93‑1 (MQ=255)

                                                                                       gcgcCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGt  <  2:3118971/93‑1 (MQ=255)

AGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGT  >  NC_000913/3853484‑3853663

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