BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	949,354	A→G	Y155C (TAT→TGT)	ycaK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	949,354	0	A	G	79.0%	419.3 / 121.8	219	Y155C (TAT→TGT)	ycaK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (29/17); new base G (82/91); total (111/108)
Fisher's exact test for biased strand distribution p-value = 6.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 221 total aligned reads displayed.

TTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTAT  >  NC_000913/949262‑949446

 tGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTAt                                                                                             <  2:643060/93‑1 (MQ=255)

     tGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTg                                                                                         <  2:892270/93‑1 (MQ=255)

     tGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTg                                                                                         >  2:3721955/1‑93 (MQ=255)

      ggCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTgg                                                                                        <  1:2592370/93‑1 (MQ=255)

       gcgcTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGt                                                                                       <  1:5568113/93‑1 (MQ=255)

          cTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTAtt                                                                                    >  2:7447539/1‑93 (MQ=255)

           tGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTg                                                                                   <  2:4845880/93‑1 (MQ=255)

               tGGAGGAGACAAAGAATCGTTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGa                                                                               <  1:988499/93‑1 (MQ=255)

               tGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGa                                                                               <  1:20916/93‑1 (MQ=255)

               tGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGa                                                                               <  1:3442322/93‑1 (MQ=255)

               tGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGa                                                                               <  1:2169488/93‑1 (MQ=255)

                ggaggaGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGAt                                                                              <  2:7695/93‑1 (MQ=255)

                ggaggaGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGAt                                                                              <  1:3077880/93‑1 (MQ=255)

                  aggagACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGc                                                                            <  1:48538/93‑1 (MQ=255)

                   ggagACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGcc                                                                           <  1:3004634/93‑1 (MQ=255)

                   ggagACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATACGTGCAGTTGTCTTGGTATTGAAGATGcc                                                                           <  2:5960873/93‑1 (MQ=255)

                    gagaCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCg                                                                          <  2:2509194/93‑1 (MQ=255)

                    gagaCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCg                                                                          >  2:2352030/1‑93 (MQ=255)

                    gagaCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCg                                                                          <  2:6915155/93‑1 (MQ=255)

                      gaCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGAt                                                                        >  1:2110421/1‑93 (MQ=255)

                      gaCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGAt                                                                        >  1:2970693/1‑93 (MQ=255)

                       aCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATg                                                                       <  1:3255783/93‑1 (MQ=255)

                       aCAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATg                                                                       >  1:5412048/1‑93 (MQ=255)

                          aaGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCa                                                                    <  1:5953403/93‑1 (MQ=255)

                          aaGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCa                                                                    >  2:1721979/1‑93 (MQ=255)

                             aaTCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACtt                                                                 >  1:3592808/1‑93 (MQ=255)

                             aaTCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACtt                                                                 <  1:5664997/93‑1 (MQ=255)

                             aaTCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACtt                                                                 <  2:420772/93‑1 (MQ=255)

                              aTCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACttt                                                                <  2:3277370/93‑1 (MQ=255)

                                   ttGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTtgt                                                           <  1:5822220/93‑1 (MQ=255)

                                      tCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTa                                                        <  1:2022937/93‑1 (MQ=255)

                                      tCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTa                                                        <  1:1028153/93‑1 (MQ=255)

                                         aGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATa                                                     <  1:4034421/93‑1 (MQ=255)

                                          gATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATAc                                                    <  2:6204541/93‑1 (MQ=255)

                                           aTGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACa                                                   <  2:7480008/93‑1 (MQ=255)

                                           aTGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACa                                                   <  1:4079155/93‑1 (MQ=255)

                                             gggCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGt                                                 <  1:590005/93‑1 (MQ=255)

                                               gCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTgg                                               <  2:2598605/93‑1 (MQ=255)

                                                 tGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTa                                             >  1:7386310/1‑93 (MQ=255)

                                                 tGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTa                                             >  1:416135/1‑93 (MQ=255)

                                                 tGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTa                                             >  2:6542792/1‑93 (MQ=255)

                                                  gggAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTAt                                            <  1:4658147/93‑1 (MQ=255)

                                                  gggAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTAt                                            >  2:5122047/1‑93 (MQ=255)

                                                   ggAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTAtt                                           >  2:1179302/1‑93 (MQ=255)

                                                    gAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTc                                          <  1:1584140/93‑1 (MQ=255)

                                                       aaaaaTATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGAt                                       >  2:6269885/1‑93 (MQ=255)

                                                       aaaaaTATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGAt                                       >  2:2851049/1‑93 (MQ=255)

                                                       aaaaaTATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGAt                                       <  1:6665568/93‑1 (MQ=255)

                                                         aaaTATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATgg                                     <  2:1865132/93‑1 (MQ=255)

                                                          aaTATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATggg                                    <  2:3236584/93‑1 (MQ=255)

                                                          aaTATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATggg                                    >  1:6145015/1‑93 (MQ=255)

                                                          aaTATAAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATggg                                    >  2:2239300/1‑93 (MQ=255)

                                                          aaTATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATggg                                    >  2:6741363/1‑93 (MQ=255)

                                                            tataAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGga                                  >  2:3371461/1‑93 (MQ=255)

                                                            tataAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGga                                  >  2:6640148/1‑93 (MQ=255)

                                                            tataAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGga                                  <  2:3977305/93‑1 (MQ=255)

                                                             ataAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaa                                 <  2:587356/93‑1 (MQ=255)

                                                              taAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaag                                >  2:6801965/1‑93 (MQ=255)

                                                              taAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaag                                >  2:6842256/1‑93 (MQ=255)

                                                              taAGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTCTGTAATACAGTGGTATTCGATGGGgaag                                >  1:3401125/1‑93 (MQ=255)

                                                              taAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaag                                >  2:6024919/1‑93 (MQ=255)

                                                               aaGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaaga                               <  1:1168804/93‑1 (MQ=255)

                                                               aaGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaaga                               >  1:2442341/1‑93 (MQ=255)

                                                               aaGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaaga                               >  1:4272832/1‑93 (MQ=255)

                                                                aGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaagaa                              >  1:1337958/1‑93 (MQ=255)

                                                                aGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaagaa                              <  2:6069076/93‑1 (MQ=255)

                                                                aGCGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGgaagaa                              <  2:4037170/93‑1 (MQ=255)

                                                                  cGATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACt                            >  2:3123708/1‑93 (MQ=255)

                                                                   gATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACtt                           >  1:2548198/1‑93 (MQ=255)

                                                                   gATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACtt                           >  1:3889915/1‑93 (MQ=255)

                                                                   gATTATTTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACtt                           >  1:467432/1‑93 (MQ=255)

                                                                   gATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACtt                           >  1:6866710/1‑93 (MQ=255)

                                                                   gATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACtt                           >  1:6399440/1‑93 (MQ=255)

                                                                    attattTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTc                          >  2:1343783/1‑93 (MQ=255)

                                                                    attattTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTc                          >  1:488982/1‑93 (MQ=255)

                                                                     ttattTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCa                         >  1:3110063/1‑93 (MQ=255)

                                                                      tattTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAc                        <  2:1086355/93‑1 (MQ=255)

                                                                      tattTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAc                        <  2:5277525/93‑1 (MQ=255)

                                                                      tattTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAc                        >  1:7446841/1‑93 (MQ=255)

                                                                       attTAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAcg                       >  2:699120/1‑93 (MQ=255)

                                                                        tttAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAcgc                      <  2:5414165/93‑1 (MQ=255)

                                                                         ttAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAcgcg                     >  1:6111483/1‑93 (MQ=255)

                                                                         ttAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAcgcg                     <  1:2866676/93‑1 (MQ=255)

                                                                         ttAAAAAATATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCAcgcg                     <  1:3689112/93‑1 (MQ=255)

                                                                             aaaaTATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCt                 >  2:406438/1‑93 (MQ=255)

                                                                               aaTATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCtat               >  1:6281652/1‑93 (MQ=255)

                                                                               aaTATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCtat               >  1:1651147/1‑93 (MQ=255)

                                                                               aaTATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCtat               >  1:5232856/1‑93 (MQ=255)

                                                                               aaTATGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCtat               <  1:6001809/93‑1 (MQ=255)

                                                                                 tatGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCtatta             >  1:1754963/1‑93 (MQ=255)

                                                                                  atGTGCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCtattat            <  1:6745583/93‑1 (MQ=255)

                                                                                    gtgCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCa          <  1:2309725/93‑1 (MQ=255)

                                                                                    gtgCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCa          >  1:5493956/1‑93 (MQ=255)

                                                                                    gtgCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCa          >  1:78231/1‑93 (MQ=255)

                                                                                      gCAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGt        >  2:7048420/1‑93 (MQ=255)

                                                                                       cAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTc       <  2:1364219/93‑1 (MQ=255)

                                                                                       cAGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTc       >  1:3801790/1‑93 (MQ=255)

                                                                                       cAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTc       >  2:4537721/1‑93 (MQ=255)

                                                                                        aGTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCg      >  1:7283254/1‑93 (MQ=255)

                                                                                         gTTGTCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGt     >  1:6808837/1‑93 (MQ=255)

TTGGGTGGCGCTGGTTGGAGGAGACAAAGAATCATTTGTCCAGATGGGCTGGGAAAAAAATATAAGCGATTATTTAAAAAATATGTGCAGTTATCTTGGTATTGAAGATGCCGATGTCACTTTCTTGTGTAATACAGTGGTATTCGATGGGGAAGAACTTCACGCGAGCTATTATCAGTCGTTAT  >  NC_000913/949262‑949446

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