BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,382,820	T→C	62.1%	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	62.1%	111.7 / 126.4	130	Y292H (TAT→CAT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand): ref base T (21/28); new base C (45/35); total (67/63)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 131 total aligned reads displayed.

TTTACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATAAGA  >  NC_000913/1382725‑1382919

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

tttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGc                                                                                                 <  2:3327075/100‑1 (MQ=255)

 ttACGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGCCAGTGATAAATTCCATCATGCa                                                                                                <  1:1961135/100‑1 (MQ=255)

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

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

   aCGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAgg                                                                                              <  1:3165686/100‑1 (MQ=255)

    cGGTTATGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGgt                                                                                             >  2:2740485/1‑100 (MQ=255)

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

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

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

         aTGATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCgg                                                                                        <  2:2113382/100‑1 (MQ=255)

           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCg                                                                                      <  1:421497/100‑1 (MQ=255)

           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                      <  2:2945942/100‑1 (MQ=255)

           gATTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCg                                                                                      <  2:3086438/100‑1 (MQ=255)

            aTTTCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGt                                                                                     >  2:2874952/1‑100 (MQ=255)

             tttCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTa                                                                                    >  2:1234011/1‑100 (MQ=255)

               tCCAGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTaaa                                                                                  >  2:2170744/1‑100 (MQ=255)

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

                  aGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTa                                                                               <  2:3330200/100‑1 (MQ=255)

                  aGAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCAATCATGCAGGTGTCGGCGTAAATTa                                                                               >  2:2480657/1‑100 (MQ=255)

                   gAACGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTAc                                                                              >  1:1077948/1‑100 (MQ=255)

                   gAACGGACTTTCCGTTTCGCTGGAATACGAGTTTGAGAGGCAGGATCACGACGAAGGCGATAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTAc                                                                              <  1:2101341/100‑1 (MQ=255)

                    aaCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACt                                                                             >  1:403974/1‑100 (MQ=255)

                     aCGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCGATCATGCAGGTGTCGGCGTAAATTACTc                                                                            >  2:238121/1‑100 (MQ=255)

                      cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                           >  1:3186087/1‑100 (MQ=255)

                      cGGACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCg                                                                           <  2:2738359/100‑1 (MQ=255)

                       ggACTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGt                                                                          <  1:2501867/100‑1 (MQ=255)

                          cTTTCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                       <  2:260128/100‑1 (MQ=255)

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

                          cTTTCCGTTTCGCTGGAATACGCGTTCGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCt                                                                       <  2:500510/100‑1 (MQ=255)

                           tttCCGTTTCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTg                                                                      >  2:3278003/1‑100 (MQ=255)

                           tttCCGTTTCGCTGGAATACGAGTTTTAGTATCAGGATCACGACGAAGGCGACAGTGATATATTCCATCATGCAGGTTTCGCCTTAAATTACTCGTCCTg                                                                      <  1:1416430/100‑1 (MQ=255)

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

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

                                 tttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATg                                                                <  1:1534010/100‑1 (MQ=255)

                                 tttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATg                                                                >  2:1565838/1‑100 (MQ=255)

                                 tttCGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATg                                                                <  2:1744460/100‑1 (MQ=255)

                                    cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATgg                                                               >  1:1633454/1‑98 (MQ=255)

                                    cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTACATTACTCGTTCTGATAATGTGc                                                             >  1:111708/1‑100 (MQ=255)

                                    cGCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                             >  1:3111289/1‑100 (MQ=255)

                                     ggtGGAATACGCGTTTGAGTGGCAGTTTCACGACGAAGGCGACAGTGATAAATTCCATCATTCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGc                                                             <  1:143247/97‑1 (MQ=255)

                                     gCTGGAATACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCt                                                            >  1:908023/1‑100 (MQ=255)

                                          aaTACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAAtt                                                       >  1:1731662/1‑100 (MQ=255)

                                           aTACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTg                                                      <  2:1500011/100‑1 (MQ=255)

                                           aTACGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTg                                                      >  2:2255411/1‑100 (MQ=255)

                                             aCGCGTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGcc                                                    <  2:1060386/100‑1 (MQ=255)

                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTATATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                   >  1:3095898/1‑100 (MQ=255)

                                              cgcgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCg                                                   >  1:3149709/1‑100 (MQ=255)

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

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

                                                cgTTTGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGa                                                 <  1:1820067/100‑1 (MQ=255)

                                                  tttGAGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATg                                               >  2:2321460/1‑100 (MQ=255)

                                                      aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATCACTCGTTCTGATAATGGGCTAAATTGCCGGATgcggc                                           <  1:1811031/100‑1 (MQ=255)

                                                      aGTGGCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACACCCTCTGATAATGGGCTAAATTGCCGGATgcggc                                           >  1:514046/1‑100 (MQ=255)

                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        >  2:1896188/1‑100 (MQ=255)

                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        >  1:1007878/1‑100 (MQ=255)

                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        <  1:566311/100‑1 (MQ=255)

                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        <  2:2994804/100‑1 (MQ=255)

                                                         ggCAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGcgcg                                        <  2:2639219/100‑1 (MQ=255)

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

                                                           cAGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAg                                      <  1:1747021/100‑1 (MQ=255)

                                                            aGGATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGt                                     >  2:892053/1‑100 (MQ=255)

                                                             ggATCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTa                                    >  2:1597837/1‑100 (MQ=255)

                                                                tCACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACtt                                 <  1:414823/100‑1 (MQ=255)

                                                                 cACGACGAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACttt                                <  1:2457486/100‑1 (MQ=255)

                                                                  acgacgAAGGCGACAGTGATAAATTCCATGATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTa                               >  1:3251164/1‑100 (MQ=255)

                                                                  acgacgAAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTa                               >  1:2465779/1‑100 (MQ=255)

                                                                   cgacgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTAt                              <  1:1613365/100‑1 (MQ=255)

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

                                                                      cgaAGGCGACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCg                           <  2:2369571/100‑1 (MQ=255)

                                                                      cgaAGGCGACAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCg                           >  2:3387728/1‑100 (MQ=255)

                                                                          ggCTACAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCt                       <  1:101527/100‑1 (MQ=255)

                                                                               cAGTGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa                  >  1:2214980/1‑100 (MQ=255)

                                                                               cAGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATaaa                  <  2:847503/100‑1 (MQ=255)

                                                                                aGTGGTAAATTCCATTATGCAAGTGTCGGCTTCACGTACTCGTTCCGATAATGGGCTAAATTGCCGGATGCGGCGAGAGTACTTTATCCGATCTATAAAt                 <  1:35196/100‑1 (MQ=255)

                                                                                aGTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAAt                 <  1:296458/100‑1 (MQ=255)

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

                                                                                 gTGATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATg                <  1:148152/100‑1 (MQ=255)

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

                                                                                  tGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGt               <  2:898722/100‑1 (MQ=255)

                                                                                  tGATAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCAGATCTATAAATGt               <  2:1333065/100‑1 (MQ=255)

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

                                                                                   gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAAt                 >  1:1665262/1‑97 (MQ=255)

                                                                                   gATAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTCATCCGATCTATAAATGTa              <  1:2201395/100‑1 (MQ=255)

                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg             >  1:2878094/1‑100 (MQ=255)

                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg             >  2:602046/1‑100 (MQ=255)

                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg             >  2:2795521/1‑100 (MQ=255)

                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg             >  1:1187502/1‑100 (MQ=255)

                                                                                    aTAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAg             <  2:1980510/100‑1 (MQ=255)

                                                                                     tAAATTCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAgg            >  2:1296218/1‑100 (MQ=255)

                                                                                     tAAATTCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAgg            >  2:1905510/1‑100 (MQ=255)

                                                                                       aaTTCCATCATGCAGGTGTCGGAGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGcc          >  1:1962769/1‑100 (MQ=255)

                                                                                         ttCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg        <  1:304292/100‑1 (MQ=255)

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

                                                                                         ttCCATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCgg        <  1:13315/100‑1 (MQ=255)

                                                                                          tCCATTATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGa       >  1:2638032/1‑100 (MQ=255)

                                                                                           ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt      >  2:3133100/1‑100 (MQ=255)

                                                                                           ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt      <  1:1565831/100‑1 (MQ=255)

                                                                                           ccATCATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGAt      <  2:2866789/100‑1 (MQ=255)

                                                                                               tATGCAGGTGTCGGCGTAAATTACTCGTTCTGATAATGGGCTAAATTGCCGGATGCGGCGCGAGTACTTTATCCGATCTATAAATGTAGGCCGGATAAGa  <  2:1226569/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` `< 32 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`