BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,382,820	T→C	24.0%	Y292H (TAT→CAT)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,820	0	T	C	24.0%	221.3 / 65.0	124	Y292H (TAT→CAT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base T (52/42); new base C (17/13); total (69/55)
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 125 total aligned reads displayed.

TTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATAAGA  >  NC_000913/1382725‑1382919

tttcctgtcaTGCTTTCCAGAACGGACCTTCCGTTTTCCTCAAATACGCGTTTGAGTGTCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGc                                                                                                 <  2:3463902/91‑1 (MQ=255)

tttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGc                                                                                                 <  2:3583977/100‑1 (MQ=255)

tttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGc                                                                                                 <  1:708874/100‑1 (MQ=255)

 ttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCa                                                                                                >  2:3744187/1‑100 (MQ=255)

  tACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAg                                                                                               <  1:2429126/100‑1 (MQ=255)

    cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgt                                                                                             >  2:1130358/1‑100 (MQ=255)

    cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGAAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAgg                                                                                              <  2:515444/99‑1 (MQ=255)

      gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgtgt                                                                                           >  2:1439010/1‑100 (MQ=255)

      gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgtgt                                                                                           >  2:3046833/1‑100 (MQ=255)

      gTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGgtgt                                                                                           >  1:353521/1‑100 (MQ=255)

        tATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCg                                                                                         >  2:449480/1‑100 (MQ=255)

         aTGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCgg                                                                                        >  2:2571521/1‑100 (MQ=255)

          tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGc                                                                                       >  1:2108629/1‑100 (MQ=255)

          tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGc                                                                                       >  1:220551/1‑100 (MQ=255)

          tGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGc                                                                                       >  1:1327030/1‑100 (MQ=255)

           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCg                                                                                      >  2:2462668/1‑100 (MQ=255)

           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTACGCAGGTGTCGGCg                                                                                      <  1:3369445/100‑1 (MQ=255)

            aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGt                                                                                     <  2:3416615/100‑1 (MQ=255)

             tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTa                                                                                    >  1:1558127/1‑100 (MQ=255)

             tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTa                                                                                    >  1:2605605/1‑100 (MQ=255)

             tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                                    <  1:1798955/100‑1 (MQ=255)

               tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTaaa                                                                                  >  2:291350/1‑100 (MQ=255)

               tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGAGAAATTCCATCATGCAGGTGTCGGCGTaa                                                                                   >  1:4017998/1‑99 (MQ=255)

               tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGAAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTaaa                                                                                  >  1:3875705/1‑100 (MQ=255)

                ccAGAACGGACTTACCGTTTCGCTGGCATACGCGTTTGTGTGGCAGGATCACGACGAAGGCGACAGTGCTAACTTCCATTATGCAGGTGTCGGCGTACAt                                                                                 >  2:2553136/1‑100 (MQ=255)

                 cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGCTAAATTCCATTATGCAGGTGTACGAGTAAAtt                                                                                >  1:1787949/1‑100 (MQ=255)

                 cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAAtt                                                                                >  1:2124739/1‑100 (MQ=255)

                 cAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAAtt                                                                                <  1:2399382/100‑1 (MQ=255)

                   gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTAc                                                                              >  1:1430933/1‑100 (MQ=255)

                      cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTc                                                                            >  2:1882825/1‑99 (MQ=255)

                      cGGACTTTCCGTTTCGCTGGAATACGCGTTCGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCg                                                                           >  1:3128421/1‑100 (MQ=255)

                       ggACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGTGTAAATTACTCGt                                                                          <  2:1442946/100‑1 (MQ=255)

                        gACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGtt                                                                         >  1:894807/1‑100 (MQ=255)

                        gACTTTCCGTTTCGCTGGAATACGCGTTTCGGTGGCAGGATCACGACGCAGGCTAAAGTGATAGATTCCAATACGCAGGTTTCGGCGTAAATTACTc                                                                            >  2:1613498/1‑97 (MQ=255)

                          cTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                       >  1:4035929/1‑100 (MQ=255)

                           tttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTg                                                                      <  1:1711637/100‑1 (MQ=255)

                            ntCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGa                                                                     >  2:3142298/2‑100 (MQ=255)

                             tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                    <  2:3535560/100‑1 (MQ=255)

                             tCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGAt                                                                    <  1:3179897/100‑1 (MQ=255)

                               cGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATaa                                                                  >  2:699003/1‑100 (MQ=255)

                               cGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATaa                                                                  >  2:3724034/1‑100 (MQ=255)

                                gTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAAt                                                                 >  2:417174/1‑100 (MQ=255)

                                    cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                             >  1:2447365/1‑100 (MQ=255)

                                        ggAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTaaa                                                         >  1:3825033/1‑100 (MQ=255)

                                            tACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGc                                                     >  1:2916248/1‑100 (MQ=255)

                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                   >  1:2049949/1‑100 (MQ=255)

                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                   >  1:3441973/1‑100 (MQ=255)

                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCTTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                   <  2:1054383/100‑1 (MQ=255)

                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                   <  2:627721/100‑1 (MQ=255)

                                               gcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCgg                                                  <  2:3687538/100‑1 (MQ=255)

                                               gcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCgg                                                  <  2:1314849/100‑1 (MQ=255)

                                               gcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCgg                                                  <  1:3411595/100‑1 (MQ=255)

                                                cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                                 >  1:3548934/1‑100 (MQ=255)

                                                cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAAGTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                                 <  2:3774575/100‑1 (MQ=255)

                                                 gTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGAt                                                >  1:3049974/1‑100 (MQ=255)

                                                     gAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcgg                                            <  1:858317/100‑1 (MQ=255)

                                                      aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggc                                           >  1:151198/1‑100 (MQ=255)

                                                       gTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATACTGGGCTAAATTACCGGATgcggcg                                          >  2:2738725/1‑100 (MQ=255)

                                                         gtcAGGATCACGACGAAGTCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        <  1:3446013/98‑1 (MQ=255)

                                                         ggCAGGATCAGGACGAAGGCGACAGTGATAAATGCCATCATGCAGGTGGCGGCGTAAATTACTCGTTCTGATATTGGGCTAAATTGCCGGATgcggc                                           >  2:3611873/1‑97 (MQ=255)

                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        <  2:89959/100‑1 (MQ=255)

                                                           cacgATAACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                      >  2:2468848/4‑100 (MQ=255)

                                                           cAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                      <  1:668074/100‑1 (MQ=255)

                                                           cAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTCCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                      <  1:3187850/100‑1 (MQ=255)

                                                           cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                      >  1:437805/1‑100 (MQ=255)

                                                             ggATCACGACGAAGGCGACAGTGATTAATTCCATTATGCAGGTGTTTGCGTAAATTACTCGTTCTGTTAATGGGCTAAATTGCCGGATGCGGCGCGAGTa                                    <  2:227035/100‑1 (MQ=255)

                                                             ggATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTa                                    >  1:3883989/1‑100 (MQ=255)

                                                               aTCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGCGTACt                                  >  2:953432/1‑100 (MQ=255)

                                                                tCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                                 >  2:1881605/1‑100 (MQ=255)

                                                                tCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                                 >  2:652700/1‑100 (MQ=255)

                                                                tCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                                 <  1:1334434/100‑1 (MQ=255)

                                                                 cACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACttt                                <  2:735293/100‑1 (MQ=255)

                                                                   cgacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                              >  1:3218957/1‑100 (MQ=255)

                                                                   cgacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                              >  1:768365/1‑100 (MQ=255)

                                                                   cgacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                              >  2:2965300/1‑100 (MQ=255)

                                                                    gacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATc                             <  1:2611308/100‑1 (MQ=255)

                                                                     acgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATcc                            >  2:3218491/1‑100 (MQ=255)

                                                                     acgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATcc                            <  1:3684772/100‑1 (MQ=255)

                                                                        aaGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGAt                         <  1:1284774/100‑1 (MQ=255)

                                                                         aGGCGACAGTGATAAATTCCATTGTGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                        <  2:2363367/100‑1 (MQ=255)

                                                                         aGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                        <  1:1440733/100‑1 (MQ=255)

                                                                         aGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                        <  2:849744/100‑1 (MQ=255)

                                                                         aGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATc                        <  1:2899909/100‑1 (MQ=255)

                                                                          ggCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                       <  1:2819320/100‑1 (MQ=255)

                                                                             gACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCtata                    <  1:1216902/100‑1 (MQ=255)

                                                                              aCAGTGATAAATTCCATTATGCAGGCGTCGGCGTAAATTACTCGTCCTGATAATGGGCGTAATTTCCGGATGCGACGCGAGCACTTTATCCGCTCtatca                   <  2:3397366/100‑3 (MQ=255)

                                                                               cAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa                  <  2:912569/100‑1 (MQ=255)

                                                                               cAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTAGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa                  <  2:2126007/100‑1 (MQ=255)

                                                                                aGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAAt                 <  1:792624/100‑1 (MQ=255)

                                                                                 gTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg                >  1:1612954/1‑100 (MQ=255)

                                                                                 gTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg                <  2:4055311/100‑1 (MQ=255)

                                                                                 gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg                >  1:2640435/1‑100 (MQ=255)

                                                                                   gATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa              >  1:706557/1‑100 (MQ=255)

                                                                                   gATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTa              >  1:3428007/1‑100 (MQ=255)

                                                                                    aTAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg             <  1:790735/100‑1 (MQ=255)

                                                                                       aaTTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGcc          >  1:2538212/1‑100 (MQ=255)

                                                                                       aaTTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGc           >  2:2686248/1‑99 (MQ=255)

                                                                                         ttCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg        >  2:3790611/1‑100 (MQ=255)

                                                                                          tCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGa       <  2:2107760/100‑1 (MQ=255)

                                                                                               tATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATAAGa  <  1:99547/100‑1 (MQ=255)

TTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATAAGA  >  NC_000913/1382725‑1382919

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