BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	335,882	T→C	I245T (ATC→ACC)	yahE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	335,882	0	T	C	99.5%	649.6 / NA	199	I245T (ATC→ACC)	yahE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (112/86); total (113/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

CATTATTAAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACATCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTT  >  NC_000913/335792‑335971

caTTATTAAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACAccc                                                                                          >  2:5496335/1‑93 (MQ=255)

 attattAAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACAcccc                                                                                         <  1:1138306/93‑1 (MQ=255)

 attattAAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACAcccc                                                                                         <  2:4721077/93‑1 (MQ=255)

  ttattaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCt                                                                                        >  2:460380/1‑93 (MQ=255)

   tattaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTg                                                                                       <  2:6096128/93‑1 (MQ=255)

   tattaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTg                                                                                       >  2:1185462/1‑93 (MQ=255)

    attaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTgg                                                                                      <  1:3950978/93‑1 (MQ=255)

    attaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTgg                                                                                      >  1:629852/1‑93 (MQ=255)

    attaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTgg                                                                                      >  2:400076/1‑93 (MQ=255)

     ttaAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTggg                                                                                     <  2:6015238/93‑1 (MQ=255)

      taAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGa                                                                                    >  2:5627334/1‑93 (MQ=255)

        aGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAAc                                                                                  >  2:1182674/1‑93 (MQ=255)

         gTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACg                                                                                 >  1:1231789/1‑93 (MQ=255)

           gCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCa                                                                               >  1:1453950/1‑93 (MQ=255)

            ccATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCaa                                                                              >  1:2274263/1‑93 (MQ=255)

             cATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAAc                                                                             >  2:1210250/1‑93 (MQ=255)

             cATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAAc                                                                             <  1:375251/93‑1 (MQ=255)

             cATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAAc                                                                             >  1:4834035/1‑93 (MQ=255)

              aTAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACt                                                                            >  1:292959/1‑93 (MQ=255)

               tAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCGCAACATTATTTATGACACCCCTGGGAACGCAACTc                                                                           <  1:1025501/93‑1 (MQ=255)

                aaCGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACATCATTATTTATGACACCCCTGGGAACGCAACTCa                                                                          <  1:497752/93‑1 (MQ=255)

                aaCGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGGACGCAACTCa                                                                          <  2:3373067/93‑1 (MQ=255)

                aaCGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCa                                                                          >  1:3083926/1‑93 (MQ=255)

                aaCGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCa                                                                          <  1:787877/93‑1 (MQ=255)

                 aCGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATTGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAg                                                                         >  1:4088692/1‑93 (MQ=255)

                 aCGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAg                                                                         <  2:5159128/93‑1 (MQ=255)

                   gCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGc                                                                       >  2:1910990/1‑93 (MQ=255)

                    cTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCa                                                                      >  1:2021779/1‑93 (MQ=255)

                    cTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCa                                                                      >  2:5803415/1‑93 (MQ=255)

                     tGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCaa                                                                     >  2:1008909/1‑93 (MQ=255)

                     tGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCaa                                                                     >  1:331137/1‑93 (MQ=255)

                      ggAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAAt                                                                    <  2:3083194/93‑1 (MQ=255)

                       gAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATa                                                                   <  1:3769389/93‑1 (MQ=255)

                       gAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATa                                                                   <  2:3528424/93‑1 (MQ=255)

                          gCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑Gaa                                                               <  1:5998433/93‑1 (MQ=255)

                           ccGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑Gaaa                                                              <  2:3869845/93‑1 (MQ=255)

                            cGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑Gaaaa                                                             >  1:1388573/1‑93 (MQ=255)

                              cATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑Gaaaaaa                                                           <  1:2908191/93‑1 (MQ=255)

                                ttACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCCATA‑GAAAAAAtt                                                         >  2:339988/1‑93 (MQ=255)

                                 tACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTa                                                        >  2:5953159/1‑93 (MQ=255)

                                  acaacaACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTaa                                                       >  1:3645223/1‑93 (MQ=255)

                                   caacaaCGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTaaa                                                      <  2:4909131/93‑1 (MQ=255)

                                    aacaacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAAc                                                     <  2:524884/93‑1 (MQ=255)

                                     acaacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACa                                                    <  2:288723/93‑1 (MQ=255)

                                     acaacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACa                                                    <  2:1572299/93‑1 (MQ=255)

                                     acaacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACa                                                    <  2:2598687/93‑1 (MQ=255)

                                      caacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACat                                                   >  1:3044948/1‑93 (MQ=255)

                                      caacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACat                                                   <  2:3678811/93‑1 (MQ=255)

                                       aacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACata                                                  <  2:3661865/93‑1 (MQ=255)

                                       aacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACata                                                  >  1:4623576/1‑93 (MQ=255)

                                        acGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACCCCCCTGGGAACGCAACTCAGGCAATAGGAAAAAATTAAACata                                                  <  2:2153110/93‑1 (MQ=255)

                                        acGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACatat                                                 >  1:5234779/1‑93 (MQ=255)

                                        acGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACatat                                                 <  1:1705656/93‑1 (MQ=255)

                                         cGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATAtt                                                <  2:2343310/93‑1 (MQ=255)

                                         cGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATAtt                                                <  2:209880/93‑1 (MQ=255)

                                         cGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATAtt                                                <  2:3896403/93‑1 (MQ=255)

                                          gCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTg                                               <  1:5788256/93‑1 (MQ=255)

                                           cTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTgg                                              >  1:4930809/1‑93 (MQ=255)

                                            tGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGc                                             >  2:3597460/1‑93 (MQ=255)

                                            tGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGc                                             <  1:3156976/93‑1 (MQ=255)

                                             gCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCt                                            >  2:4056815/1‑93 (MQ=255)

                                              ccGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTc                                           >  1:3123613/1‑93 (MQ=255)

                                                  gAAAATATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGt                                       >  2:5394684/1‑93 (MQ=255)

                                                    aaaTATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTc                                     >  1:904460/1‑93 (MQ=255)

                                                     aaTATTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTcc                                    >  1:1641282/1‑93 (MQ=255)

                                                       tatTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCgg                                  >  2:2980863/1‑93 (MQ=255)

                                                        atTCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGc                                 >  2:6183032/1‑93 (MQ=255)

                                                          tCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTg                               >  2:1386618/1‑93 (MQ=255)

                                                          tCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTg                               >  1:2632120/1‑93 (MQ=255)

                                                          tCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTg                               >  2:3460351/1‑93 (MQ=255)

                                                          tCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTg                               >  1:5047277/1‑93 (MQ=255)

                                                          tCATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTg                               <  1:4527695/93‑1 (MQ=255)

                                                           cATCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGc                              >  1:1778828/1‑93 (MQ=255)

                                                            aTCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCg                             <  1:1556979/93‑1 (MQ=255)

                                                             tCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGa                            <  1:1804282/93‑1 (MQ=255)

                                                             tCGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGa                            >  1:5466918/1‑93 (MQ=255)

                                                              cGTTTTATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGAc                           <  2:2731227/93‑1 (MQ=255)

                                                                  ttATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGc                       >  2:4980631/1‑93 (MQ=255)

                                                                   tATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCt                      >  2:5550688/1‑93 (MQ=255)

                                                                   tATTCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCt                      <  1:2478776/93‑1 (MQ=255)

                                                                     ttCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGt                    >  2:793818/1‑93 (MQ=255)

                                                                     ttCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGt                    >  2:4693747/1‑93 (MQ=255)

                                                                      tCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTa                   <  2:5186695/93‑1 (MQ=255)

                                                                      tCACAACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTa                   <  1:4425156/93‑1 (MQ=255)

                                                                       cacaACATTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTat                  <  2:1178258/93‑1 (MQ=255)

                                                                        acaacaTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatg                 >  2:2769595/1‑93 (MQ=255)

                                                                        acaacaTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatg                 >  2:5880838/1‑93 (MQ=255)

                                                                          aacaTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggc               >  2:6235643/1‑93 (MQ=255)

                                                                           acaTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggca              <  2:4176896/93‑1 (MQ=255)

                                                                           acaTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggca              >  1:2989073/1‑93 (MQ=255)

                                                                            caTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggcat             >  1:5842092/1‑93 (MQ=255)

                                                                            caTTATTTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggcat             <  1:553678/93‑1 (MQ=255)

                                                                             attattTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggcatg            <  1:2099329/93‑1 (MQ=255)

                                                                             attattTATGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggcatg            <  2:2106164/93‑1 (MQ=255)

                                                                              ttatttatGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggcatgg           >  1:687560/1‑93 (MQ=255)

                                                                              ttatttatGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTatggcatgg           >  1:47357/1‑93 (MQ=255)

                                                                                 tttatGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCg        <  1:3426583/93‑1 (MQ=255)

                                                                                 tttatGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCg        >  1:434093/1‑93 (MQ=255)

                                                                                    atGACACCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATg     >  1:2853949/1‑93 (MQ=255)

                                                                                       acacCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGtt  <  1:4231286/93‑1 (MQ=255)

CATTATTAAGTGCCATAACGCTGGAAGCCGCATTACAACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACATCCCTGGGAACGCAACTCAGGCAATA‑GAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTT  >  NC_000913/335792‑335971

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