BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,224,926	+C	coding (1052/1434 nt)	ygjI →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,224,924	1	.	C	99.0%	638.1 / ‑1.7	198	coding (1050/1434 nt)	ygjI	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (87/109); total (88/110)
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.98e-01

Only 100 of 227 total aligned reads displayed.

TTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑CC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTC  >  NC_000913/3224832‑3225011

tttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑c                                                                                         >  1:924291/1‑93 (MQ=255)

tttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑c                                                                                         <  2:1440151/93‑1 (MQ=255)

 ttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cc                                                                                        <  2:4502177/93‑1 (MQ=255)

 ttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cc                                                                                        <  2:1732530/93‑1 (MQ=255)

  tGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑ccc                                                                                       <  1:6055153/93‑2 (MQ=255)

  tGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑ccc                                                                                       >  2:2601202/1‑92 (MQ=255)

   ggTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccg                                                                                      >  2:4745854/1‑91 (MQ=255)

    gTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccga                                                                                     <  1:1938610/93‑4 (MQ=255)

     tAAGAAAACCGTCGAACTGAACGAAACCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccgat                                                                                    <  2:1073913/93‑5 (MQ=255)

     tAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccgat                                                                                    <  1:4714881/93‑5 (MQ=255)

     tAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccgat                                                                                    >  2:464669/1‑89 (MQ=255)

     tAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑cccgat                                                                                    <  2:3489748/93‑5 (MQ=255)

       aGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGc                                                                                  <  2:6014554/93‑1 (MQ=255)

        gAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCt                                                                                 >  2:4064744/1‑93 (MQ=255)

         aaaaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGTTCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTc                                                                                <  2:430612/93‑1 (MQ=255)

           aaCCGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCgg                                                                              >  1:4572658/1‑93 (MQ=255)

             ccGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGtt                                                                            >  1:3186401/1‑93 (MQ=255)

             ccGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGGTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCgggt                                                                            >  2:916190/1‑91 (MQ=255)

              cGTCGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATT‑CA‑GATGCTCGGTTc                                                                           >  2:5295580/1‑93 (MQ=255)

                 cGAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCaa                                                                        <  2:2160610/93‑1 (MQ=255)

                  gAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAAt                                                                       >  2:233272/1‑93 (MQ=255)

                  gAACTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAAt                                                                       <  2:5079475/93‑1 (MQ=255)

                   aaCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATa                                                                      <  1:2512957/93‑1 (MQ=255)

                   aaCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATa                                                                      <  2:4164203/93‑1 (MQ=255)

                    aCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATAc                                                                     <  1:1098910/93‑1 (MQ=255)

                    aCTGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATAc                                                                     <  1:2223980/93‑1 (MQ=255)

                     cTGAACGAAAACGGCGTCCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑CCAGATGCTCGGTACCAATAcc                                                                    <  2:6149191/93‑1 (MQ=255)

                      tGAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCg                                                                   >  1:5632539/1‑93 (MQ=255)

                       gAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGt                                                                  <  1:1922138/93‑1 (MQ=255)

                       gAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGt                                                                  >  1:5538756/1‑93 (MQ=255)

                       gAACGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGt                                                                  <  1:4783431/93‑1 (MQ=255)

                         aCGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGc                                                                <  1:2882021/93‑1 (MQ=255)

                          cGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCa                                                               <  1:1761662/93‑1 (MQ=255)

                          cGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCa                                                               >  2:2421296/1‑93 (MQ=255)

                          cGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCa                                                               <  1:1227127/93‑1 (MQ=255)

                          cGAAAACGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCa                                                               <  2:1911722/93‑1 (MQ=255)

                            aaaaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAgg                                                             <  1:5051997/93‑1 (MQ=255)

                            aaaaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAgg                                                             <  2:1586725/93‑1 (MQ=255)

                            aaaaCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAgg                                                             >  2:6174666/1‑93 (MQ=255)

                               aCGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATc                                                          <  2:4411528/93‑1 (MQ=255)

                                cGGCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCt                                                         >  1:4150492/1‑93 (MQ=255)

                                  gCGTTCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCtga                                                       >  2:2547090/1‑93 (MQ=255)

                                     ttCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCtgatga                                                    >  2:5737698/1‑93 (MQ=255)

                                      tCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGaa                                                   <  1:851250/93‑1 (MQ=255)

                                      tCCGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGaa                                                   <  1:6138432/93‑1 (MQ=255)

                                       ccGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAAt                                                  <  2:2774914/93‑1 (MQ=255)

                                       ccGGCGC‑GCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAAt                                                  <  2:130802/93‑1 (MQ=255)

                                         ggcgc‑gcgcaGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATAc                                                >  1:531980/1‑93 (MQ=255)

                                           gcgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACt                                               >  1:5419194/5‑93 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  1:5235577/90‑1 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  2:3794975/90‑1 (MQ=255)

                                            cgcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACta                                              <  2:3439407/90‑1 (MQ=255)

                                             gcgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACtat                                             <  2:3589208/91‑1 (MQ=255)

                                              cgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACtatt                                            >  2:2692058/2‑93 (MQ=255)

                                               gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGCTGAATACtatta                                           >  2:3006868/1‑93 (MQ=255)

                                               gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACtatta                                           >  1:3257006/1‑93 (MQ=255)

                                               gcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACtatta                                           <  1:6009255/93‑1 (MQ=255)

                                                cgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACtattat                                          <  1:5169307/93‑1 (MQ=255)

                                                 gcAGCGTTGATCCAGTGCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTattatt                                         <  2:5147636/93‑1 (MQ=255)

                                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTattatt                                         >  1:4179884/1‑93 (MQ=255)

                                                 gcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTattatt                                         >  1:5928017/1‑93 (MQ=255)

                                                   cgCGTGGATCCAGTTCCTGGTCGTCATCCCGCTGCTGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTaa                                       <  1:2005629/92‑1 (MQ=255)

                                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTaa                                       <  1:5261948/93‑1 (MQ=255)

                                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTaa                                       >  1:3488892/1‑93 (MQ=255)

                                                   aGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTaa                                       <  2:1347603/93‑1 (MQ=255)

                                                    gCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAat                                      >  2:5315343/1‑93 (MQ=255)

                                                      gTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAatat                                    <  1:2443602/93‑1 (MQ=255)

                                                       tGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATg                                   <  1:2301589/93‑1 (MQ=255)

                                                       tGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATg                                   >  2:4966349/1‑93 (MQ=255)

                                                        ggATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGa                                  <  2:906548/93‑1 (MQ=255)

                                                         gATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAc                                 >  1:169518/1‑93 (MQ=255)

                                                         gATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAc                                 >  1:1527850/1‑93 (MQ=255)

                                                         gATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAc                                 <  1:3261628/93‑1 (MQ=255)

                                                          aTCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAcc                                <  1:6228120/93‑1 (MQ=255)

                                                           tCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTCTTATTAATATGAccg                               >  1:373584/1‑93 (MQ=255)

                                                           tCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccg                               <  2:4815658/93‑1 (MQ=255)

                                                           tCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGAccg                               >  2:636070/1‑93 (MQ=255)

                                                                ttCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCa                          >  1:1016643/1‑93 (MQ=255)

                                                                ttCCTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCa                          >  1:5049698/1‑93 (MQ=255)

                                                                  ccTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGc                        <  1:5964210/93‑1 (MQ=255)

                                                                   cTGATCGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCg                       >  2:6041553/1‑93 (MQ=255)

                                                                        cGTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCAt                  <  2:3103877/93‑1 (MQ=255)

                                                                         gTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATg                 >  1:3914363/1‑93 (MQ=255)

                                                                         gTCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATg                 <  2:5454871/93‑1 (MQ=255)

                                                                          tCATCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGGCCGCCGCAGCGTCCATGc                >  1:1583090/1‑93 (MQ=255)

                                                                            aTCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCtt              <  1:312977/93‑1 (MQ=255)

                                                                             tCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTc             <  1:5822525/93‑1 (MQ=255)

                                                                             tCCCGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTc             >  1:5502423/1‑93 (MQ=255)

                                                                               ccGCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTccg           >  2:2836404/1‑93 (MQ=255)

                                                                                 gCTGATGATTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTccgcc         >  2:1410774/1‑93 (MQ=255)

                                                                                   tgatgaTTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       >  1:3161822/1‑93 (MQ=255)

                                                                                   tgatgaTTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       >  1:1308604/1‑93 (MQ=255)

                                                                                   tgatgaTTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGt       <  1:2570581/93‑1 (MQ=255)

                                                                                    gatgatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGtt      >  2:1035527/1‑93 (MQ=255)

                                                                                    gatgatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGtt      >  1:5102578/1‑93 (MQ=255)

                                                                                    gatgatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGtt      >  2:3874714/1‑93 (MQ=255)

                                                                                     atgatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTa     >  1:1863718/1‑93 (MQ=255)

                                                                                     atgatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTa     >  1:4867626/1‑93 (MQ=255)

                                                                                       gatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTAtt   >  1:283558/1‑93 (MQ=255)

                                                                                       gatTA‑TCCC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTAtt   <  1:4371663/93‑1 (MQ=255)

TTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTA‑T‑CC‑GATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTC  >  NC_000913/3224832‑3225011

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