BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	974,456	A→G	V305V (GTA→GTG)	mukF →	MukF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	974,456	0	A	G	86.4%	662.6 / 94.7	279	V305V (GTA→GTG)	mukF	MukF
Reads supporting (aligned to +/- strand): ref base A (14/24); new base G (116/125); total (130/149)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 280 total aligned reads displayed.

TTGGCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTGCGCGATGA  >  NC_000913/974365‑974548

 tGGCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCa                                                                                            <  2:3296092/93‑1 (MQ=255)

 tGGCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCa                                                                                            <  1:5749891/93‑1 (MQ=255)

  ggCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCaa                                                                                           >  2:5914017/1‑93 (MQ=255)

   gCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCaaa                                                                                          <  1:1146660/93‑1 (MQ=255)

   gCTACGACCGCCACGTACACAAATTTATTCGTACCACGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCaaa                                                                                          >  2:4526513/1‑93 (MQ=255)

     tACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAAcc                                                                                        >  2:284449/1‑93 (MQ=255)

       cGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTa                                                                                      <  2:4683846/93‑1 (MQ=255)

           cGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTAtttt                                                                                  <  2:1971316/93‑1 (MQ=255)

             ccACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTtga                                                                                <  2:6792560/93‑1 (MQ=255)

               aCGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTtgatg                                                                              <  1:2176221/93‑1 (MQ=255)

               aCGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTtgatg                                                                              <  2:608566/93‑1 (MQ=255)

               aCGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTtgatg                                                                              >  1:337313/1‑93 (MQ=255)

                cGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTtgatga                                                                             <  2:1411311/93‑1 (MQ=255)

                  tACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGc                                                                           >  2:2388484/1‑93 (MQ=255)

                     acaAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGt                                                                        <  1:7422076/93‑1 (MQ=255)

                     acaAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGt                                                                        <  2:6728940/93‑1 (MQ=255)

                      caAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTg                                                                       >  1:6733342/1‑93 (MQ=255)

                       aaaTTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTgg                                                                      >  2:5213283/1‑93 (MQ=255)

                         aTTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGgc                                                                    <  2:4385505/93‑1 (MQ=255)

                           ttATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGgcgc                                                                  >  1:3920534/1‑93 (MQ=255)

                              ttCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTaa                                                               >  1:1790933/1‑93 (MQ=255)

                               tCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAAc                                                              <  1:6118133/93‑1 (MQ=255)

                                cGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACt                                                             <  1:3527013/93‑1 (MQ=255)

                                  tACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTa                                                           <  1:727322/93‑1 (MQ=255)

                                  tACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTa                                                           <  2:760767/93‑1 (MQ=255)

                                  tACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTa                                                           <  1:3313406/93‑1 (MQ=255)

                                   aCCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTaa                                                          >  1:5499428/1‑92 (MQ=255)

                                   aCCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTAt                                                          >  2:5969525/1‑93 (MQ=255)

                                   aCCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTAt                                                          <  2:7162264/93‑1 (MQ=255)

                                    ccGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATg                                                         >  2:3066055/1‑93 (MQ=255)

                                    ccGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATg                                                         <  1:6240490/93‑1 (MQ=255)

                                    ccGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATg                                                         >  2:1496284/1‑93 (MQ=255)

                                     cgcgATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGc                                                        <  2:3275291/93‑1 (MQ=255)

                                      gcgATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGcc                                                       <  1:4196777/93‑1 (MQ=255)

                                        gatcgatATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCaa                                                     <  1:3103567/93‑1 (MQ=255)

                                        gatcgatATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCaa                                                     <  2:180905/93‑1 (MQ=255)

                                         atcgatATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAAt                                                    <  2:5759712/93‑1 (MQ=255)

                                         atcgatATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAAt                                                    <  2:1073923/93‑1 (MQ=255)

                                          tcgatATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATg                                                   >  1:5343513/1‑93 (MQ=255)

                                           cgatATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGc                                                  <  2:3591816/93‑1 (MQ=255)

                                              tatGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGa                                               >  1:1059844/1‑93 (MQ=255)

                                              tatGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGa                                               <  1:2050464/93‑1 (MQ=255)

                                               atGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGAt                                              <  1:1607378/93‑1 (MQ=255)

                                               atGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGAt                                              <  1:923974/93‑1 (MQ=255)

                                               atGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGAt                                              >  2:6950529/1‑93 (MQ=255)

                                               atGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGAt                                              <  1:5629912/93‑1 (MQ=255)

                                                tGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTAATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATc                                             >  1:2857848/1‑93 (MQ=255)

                                                tGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATc                                             >  1:3259603/1‑93 (MQ=255)

                                                 ggATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCg                                            >  2:1673457/1‑93 (MQ=255)

                                                   aTAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTc                                          <  1:401518/93‑1 (MQ=255)

                                                      aaaaCCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTctgc                                       >  1:6957706/1‑93 (MQ=255)

                                                      aaaaCCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTctgc                                       >  1:3833281/1‑93 (MQ=255)

                                                      aaaaCCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTctgc                                       <  2:1177941/93‑1 (MQ=255)

                                                          ccGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGTGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGa                                   >  2:4145441/1‑93 (MQ=255)

                                                           cgcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGat                                  <  1:2719025/93‑1 (MQ=255)

                                                           cgcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGat                                  >  1:6202730/1‑93 (MQ=255)

                                                            gcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGata                                 >  2:4024279/1‑93 (MQ=255)

                                                            gcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGata                                 <  2:4965352/93‑1 (MQ=255)

                                                            gcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGata                                 <  1:6346746/93‑1 (MQ=255)

                                                            gcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGata                                 <  1:1005529/93‑1 (MQ=255)

                                                            gcgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGata                                 <  1:2427489/93‑1 (MQ=255)

                                                             cgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGatat                                <  1:6417585/93‑1 (MQ=255)

                                                             cgTCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGatat                                <  1:1941424/93‑1 (MQ=255)

                                                               tCTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTAGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGc                              <  2:4610568/93‑1 (MQ=255)

                                                               tCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGc                              >  2:6144371/1‑93 (MQ=255)

                                                                cTTTGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCg                             >  2:3729878/1‑93 (MQ=255)

                                                                 tttGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGt                            <  2:1311070/93‑1 (MQ=255)

                                                                 tttGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGt                            <  1:6509865/93‑1 (MQ=255)

                                                                   tGCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGa                          <  2:7453555/93‑1 (MQ=255)

                                                                    gCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGAc                         <  2:4706364/93‑1 (MQ=255)

                                                                    gCTCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGAc                         <  2:3435744/93‑1 (MQ=255)

                                                                      tCAGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGa                       >  1:155623/1‑93 (MQ=255)

                                                                       cAGCGGTTAGGTCAGTCGGGGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGGCTGCTGGATATGCGTGACGaa                      >  1:792528/1‑93 (MQ=255)

                                                                        aGCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAag                     <  1:664771/93‑1 (MQ=255)

                                                                         gCGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAaga                    <  1:2959663/93‑1 (MQ=255)

                                                                         gCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAaga                    >  2:557171/1‑93 (MQ=255)

                                                                          cGGTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAagag                   >  2:5063751/1‑93 (MQ=255)

                                                                           ggTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAgaga                  <  1:4106982/93‑1 (MQ=255)

                                                                            gTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGAt                 >  2:6386984/1‑93 (MQ=255)

                                                                            gTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGAt                 <  1:7494948/93‑1 (MQ=255)

                                                                            gTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGAt                 <  1:2481499/93‑1 (MQ=255)

                                                                            gTTACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGAt                 >  1:5110365/1‑93 (MQ=255)

                                                                            gTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGAt                 >  2:3669869/1‑93 (MQ=255)

                                                                            gTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGAt                 <  1:4724684/93‑1 (MQ=255)

                                                                             ttACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATg                <  1:1330870/93‑1 (MQ=255)

                                                                             ttACGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATg                <  1:6832520/93‑1 (MQ=255)

                                                                               aCGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGc              >  1:2625912/1‑93 (MQ=255)

                                                                               aCGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGc              >  1:4478084/1‑93 (MQ=255)

                                                                               aCGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGc              <  1:2000230/93‑1 (MQ=255)

                                                                               aCGTCAGTCGCTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGc              <  1:3702414/93‑1 (MQ=255)

                                                                                cGTCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCa             <  1:4916584/93‑1 (MQ=255)

                                                                                 gTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCAc            <  2:915370/93‑1 (MQ=255)

                                                                                  tCAGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACt           >  2:3504473/1‑93 (MQ=255)

                                                                                    aGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTgc         >  1:7124900/1‑93 (MQ=255)

                                                                                    aGTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTgc         >  1:1231774/1‑93 (MQ=255)

                                                                                     gTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTgcg        <  2:6312718/93‑1 (MQ=255)

                                                                                     gTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTgcg        >  2:7003606/1‑93 (MQ=255)

                                                                                     gTCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTgcg        >  2:3383794/1‑93 (MQ=255)

                                                                                      tCGGTGCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTgcgc       >  2:2329791/1‑93 (MQ=255)

                                                                                           aCAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTGCGCGATga  >  2:6617421/1‑93 (MQ=255)

TTGGCTACGACCGCCACGTACACAAATTTATTCGTACCGCGATCGATATGGATAAAAACCGCGTCTTTGCTCAGCGGTTACGTCAGTCGGTACAAACCTATTTTGATGAGCCGTGGGCGCTAACTTATGCCAATGCCGATCGTCTGCTGGATATGCGTGACGAAGAGATGGCACTGCGCGATGA  >  NC_000913/974365‑974548

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