BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,621,746	T→C	I43M (ATA→ATG)	ydeH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,746	0	T	C	98.4%	1005.2 / 0.4	316	I43M (ATA→ATG)	ydeH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/4); new base C (174/137); total (175/141)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 325 total aligned reads displayed.

CCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATT  >  NC_000913/1621656‑1621835

  gaccGTAAGGTAATTCATCGTTATCGAGTGCCCCCAGATGATCAACCCCCCGACCAAACTGGCACAGTCCATAAGAATGTTTATCCGTCAttt                                                                                       <  2:5878798/90‑1 (MQ=255)

   gaCGTAAGGCAATTGATCGTTGTCGAGTGGCCCCAGATCATCAATCCACCGACCAAACGGGCACAGTCCATAAGAATGGTTATCCGTTATTTc                                                                                      <  1:9615463/92‑1 (MQ=255)

   aaCGTAAGGTAATTCATCGATATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGGTATCCGTCATTTc                                                                                      >  2:3651132/1‑93 (MQ=255)

    aCGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCt                                                                                     >  1:9212392/1‑93 (MQ=255)

      gTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTgg                                                                                   >  2:2662303/1‑93 (MQ=255)

      gTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTgg                                                                                   <  1:1824643/93‑1 (MQ=255)

       tAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGc                                                                                  <  2:565201/93‑1 (MQ=255)

         aGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCtt                                                                                >  2:3023939/1‑93 (MQ=255)

         aGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCtt                                                                                <  2:2270892/93‑1 (MQ=255)

         aGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCtt                                                                                >  1:3809464/1‑93 (MQ=255)

          ggTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTa                                                                               <  1:7486343/93‑1 (MQ=255)

          ggTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTa                                                                               <  1:596650/93‑1 (MQ=255)

            tAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACt                                                                             <  2:2592130/93‑1 (MQ=255)

              aTTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTgg                                                                           <  1:6307672/93‑1 (MQ=255)

               ttCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGc                                                                          <  2:1981270/93‑1 (MQ=255)

                tCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATTCGTCATTTCTGGCTTACTGGCa                                                                         >  1:1051930/1‑93 (MQ=255)

                 taTCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAt                                                                        <  1:7197976/92‑1 (MQ=255)

                 cATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAt                                                                        >  2:2935236/1‑93 (MQ=255)

                   tCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtct                                                                      >  1:3128891/1‑93 (MQ=255)

                   tCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtct                                                                      <  1:8696009/93‑1 (MQ=255)

                    cGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtctc                                                                     <  2:3194275/93‑1 (MQ=255)

                     gTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtctct                                                                    <  1:8387605/93‑1 (MQ=255)

                      ttATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtctctc                                                                   >  1:9370254/1‑93 (MQ=255)

                      ttATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtctctc                                                                   <  2:1507807/93‑1 (MQ=255)

                      ttATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtctctc                                                                   <  1:2687400/93‑1 (MQ=255)

                      ttATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCAtctctc                                                                   >  1:6074496/1‑93 (MQ=255)

                       tATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCATCTCTcg                                                                  <  1:5140147/93‑1 (MQ=255)

                       tATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTcg                                                                  <  2:7081321/93‑1 (MQ=255)

                         tCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTcgcg                                                                >  2:9469429/1‑93 (MQ=255)

                         tCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTcgcg                                                                <  1:5615489/93‑1 (MQ=255)

                           gAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGAc                                                              >  1:8851869/1‑93 (MQ=255)

                            aGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGAcc                                                             >  2:6491990/1‑93 (MQ=255)

                            aGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGAcc                                                             <  1:1163062/93‑1 (MQ=255)

                                gCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGc                                                         <  1:1682040/93‑1 (MQ=255)

                                  cccAGATGATCAATCCACCGACCCAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtg                                                       <  2:5496239/93‑1 (MQ=255)

                                   ccAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgt                                                      >  1:5785942/1‑93 (MQ=255)

                                   ccAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgt                                                      >  1:9439319/1‑93 (MQ=255)

                                   ccAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgt                                                      >  2:3750788/1‑93 (MQ=255)

                                    cAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCtgtg                                                     >  1:5080494/1‑93 (MQ=255)

                                         gatCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACa                                                >  1:7829651/1‑93 (MQ=255)

                                         gatCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGAATCTCTCGCGACCACGCTGTGAAACa                                                >  1:2633134/1‑93 (MQ=255)

                                           tCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATa                                              <  1:811588/93‑1 (MQ=255)

                                           tCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATa                                              >  1:8224066/1‑93 (MQ=255)

                                           tCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATa                                              >  1:2900990/1‑93 (MQ=255)

                                            cAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATAc                                             <  1:3452981/93‑1 (MQ=255)

                                            cAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATAc                                             <  1:3263809/93‑1 (MQ=255)

                                             aaTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACt                                            >  1:6934821/1‑93 (MQ=255)

                                             aaTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACt                                            >  2:1428825/1‑93 (MQ=255)

                                             aaTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACt                                            >  1:8342301/1‑93 (MQ=255)

                                             aaTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACt                                            >  2:4633563/1‑93 (MQ=255)

                                              aTCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTc                                           <  2:7846042/93‑1 (MQ=255)

                                               tCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCa                                          >  1:6513230/1‑93 (MQ=255)

                                               tCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCa                                          >  1:4498978/1‑93 (MQ=255)

                                                ccACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCAc                                         <  1:5040924/93‑1 (MQ=255)

                                                 cACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCAcc                                        <  1:6742482/93‑1 (MQ=255)

                                                   ccGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAg                                      >  2:4176846/1‑93 (MQ=255)

                                                    cGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGc                                     >  1:104781/1‑93 (MQ=255)

                                                      aCCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCa                                   >  2:10018147/1‑93 (MQ=255)

                                                      aCCAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCa                                   <  1:112387/93‑1 (MQ=255)

                                                       ccAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCAc                                  <  2:110546/93‑1 (MQ=255)

                                                       ccAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCAc                                  >  1:2040960/1‑93 (MQ=255)

                                                        cAAACTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACt                                 >  2:2219693/1‑93 (MQ=255)

                                                         aaaCTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTg                                <  1:9860996/93‑1 (MQ=255)

                                                          aaCTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTgg                               >  2:3547992/1‑93 (MQ=255)

                                                            cTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTa                             >  1:9032202/1‑93 (MQ=255)

                                                            cTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTa                             <  2:4732027/93‑1 (MQ=255)

                                                            cTGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTa                             <  2:7519269/93‑1 (MQ=255)

                                                            cTGGCACAGTCCATAAGAATGGTTATCCGGCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTa                             <  1:2261462/93‑1 (MQ=255)

                                                             tGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAg                            <  1:465361/93‑1 (MQ=255)

                                                             tGGCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAg                            <  2:4908642/93‑1 (MQ=255)

                                                              ggCACAGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGt                           >  1:6258010/1‑93 (MQ=255)

                                                                  caGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGc                       >  1:8261672/1‑93 (MQ=255)

                                                                   aGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCa                      >  1:7793487/1‑93 (MQ=255)

                                                                   aGTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCa                      <  1:4036264/93‑1 (MQ=255)

                                                                    gTCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCAt                     >  2:5962556/1‑93 (MQ=255)

                                                                     tCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATc                    >  2:2381757/1‑93 (MQ=255)

                                                                     tCCATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATc                    <  1:1722759/93‑1 (MQ=255)

                                                                       cATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGa                  <  2:4138244/93‑1 (MQ=255)

                                                                       cATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGa                  <  1:1359759/93‑1 (MQ=255)

                                                                       cATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGa                  <  1:3862495/93‑1 (MQ=255)

                                                                       cATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGa                  >  2:1305259/1‑93 (MQ=255)

                                                                       cATAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGa                  <  1:3576629/93‑1 (MQ=255)

                                                                        aTAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGAt                 <  1:9947175/93‑1 (MQ=255)

                                                                         tAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATa                >  2:9297762/1‑93 (MQ=255)

                                                                         tAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATa                >  1:3526706/1‑93 (MQ=255)

                                                                         tAAGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATa                >  2:1636167/1‑93 (MQ=255)

                                                                           aGAATGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGc              <  2:5141834/93‑1 (MQ=255)

                                                                             aaTGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCt            <  2:6302760/93‑1 (MQ=255)

                                                                             aaTGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCt            >  1:3062741/1‑93 (MQ=255)

                                                                             aaTGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCt            >  2:337806/1‑93 (MQ=255)

                                                                               tGGTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTa          >  1:4785084/1‑93 (MQ=255)

                                                                                ggTTATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTAt         <  2:6691986/93‑1 (MQ=255)

                                                                                   tATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTga      <  1:242040/93‑1 (MQ=255)

                                                                                   tATCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTga      <  1:4214766/93‑1 (MQ=255)

                                                                                    aTCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTgag     >  2:3256149/1‑93 (MQ=255)

                                                                                    aTCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTgag     >  2:8200193/1‑93 (MQ=255)

                                                                                     tCCGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTgaga    >  2:6155414/1‑93 (MQ=255)

                                                                                      ccGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGAt   >  2:804450/1‑93 (MQ=255)

                                                                                      ccGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGAt   >  1:5343198/1‑93 (MQ=255)

                                                                                       cGTCATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGAtt  >  2:8577848/1‑93 (MQ=255)

CCGAACGTAAGGTAATTCATCGTTATCGAGTGGCCCCAGATGATCAATCCACCGACCAAACTGGCACAGTCCATAAGAATGGTTATCCGTTATTTCTGGCTTACTGGCATCTCTCGCGACCACGCTGTGAAACATACTCACCAGCCACTGGTAGTGGGCATCGATAGCCTTATTGAGATT  >  NC_000913/1621656‑1621835

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