BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,592,156	(C)_8→7	coding (3955/3978 nt)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,592,149	0	C	.	91.3%	637.5 / 47.0	172	coding (3962/3978 nt)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (11/4); new base . (71/86); total (82/90)
Fisher's exact test for biased strand distribution p-value = 5.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.81e-01

Only 100 of 197 total aligned reads displayed.

TACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCC  >  NC_000913/1592057‑1592238

tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGc                                                                                           >  1:1595984/1‑93 (MQ=255)

  cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccc                                                                                         >  2:4771792/1‑93 (MQ=255)

   aTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccc                                                                                        >  2:5075236/1‑93 (MQ=255)

      ccccGGTTCCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccccccc                                                                                     <  1:2853153/93‑1 (MQ=255)

       cccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccccccca                                                                                    <  1:5678054/93‑2 (MQ=255)

          ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                 <  2:4182610/93‑5 (MQ=255)

          ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                 >  2:1712494/1‑89 (MQ=255)

          ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                 <  2:2798229/93‑5 (MQ=255)

           gTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCa                                                                               >  1:4978129/1‑93 (MQ=255)

           gTTCCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCa                                                                               <  1:5376020/93‑1 (MQ=255)

              tCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATaa                                                                            <  1:385719/93‑1 (MQ=255)

               ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAAc                                                                           >  2:5734673/1‑93 (MQ=255)

                  cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGtt                                                                        >  2:2590044/1‑93 (MQ=255)

                   ccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTc                                                                       <  2:5280913/93‑1 (MQ=255)

                   ccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTc                                                                       >  1:6062998/1‑93 (MQ=255)

                    cTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCa                                                                      <  2:6419537/93‑1 (MQ=255)

                     ttCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCat                                                                     <  1:7914562/93‑1 (MQ=255)

                       cGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCatat                                                                   >  1:507078/1‑93 (MQ=255)

                        gttATTCATTTGAACCGGGCCGCCATTCGCCGGGTTGCCGCCTATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATc                                                                  <  2:6244869/91‑1 (MQ=255)

                        gCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATc                                                                  <  2:4012906/93‑1 (MQ=255)

                         cTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑GCTCCCCAGGCCTAACGTTCATATCa                                                                 >  2:3398102/1‑93 (MQ=255)

                         cTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCa                                                                 >  2:1657852/1‑93 (MQ=255)

                          tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATctc                                                                <  2:4638373/93‑3 (MQ=255)

                          tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAc                                                                <  2:2733501/93‑1 (MQ=255)

                          tATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAc                                                                >  2:2163462/1‑93 (MQ=255)

                           aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAcc                                                               <  1:849697/93‑1 (MQ=255)

                           aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAcc                                                               <  2:175982/93‑1 (MQ=255)

                           aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAcc                                                               <  1:229838/93‑1 (MQ=255)

                           aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAcc                                                               >  2:1860907/1‑93 (MQ=255)

                            ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCt                                                              >  1:128085/1‑93 (MQ=255)

                            ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCt                                                              >  1:6452549/1‑93 (MQ=255)

                            ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCt                                                              >  1:3227158/1‑93 (MQ=255)

                               aTTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGc                                                            >  1:1801611/1‑93 (MQ=255)

                                tttGACCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                          <  1:3731875/93‑1 (MQ=255)

                                tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                          >  2:6263132/1‑93 (MQ=255)

                                tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                          >  2:5026800/1‑93 (MQ=255)

                                tttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAg                                                          >  1:4779082/1‑93 (MQ=255)

                                   gAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                       >  2:1159059/1‑93 (MQ=255)

                                    aaCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAgtggt                                                      <  2:6868790/93‑1 (MQ=255)

                                     aCCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAgtggtg                                                     <  2:824148/93‑1 (MQ=255)

                                      ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGa                                                    <  2:5109878/93‑1 (MQ=255)

                                      ccGTGCCGCCATTCTCCAGCTTGCCGCCGATAGTGGGTTTAGCGACCCGCAGT‑ACCCCCCCAAGCATAAAGTTCATATCACCTGCAGCGGTGa                                                    >  2:601960/1‑93 (MQ=255)

                                        gTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACa                                                  >  2:7701446/1‑93 (MQ=255)

                                        gTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACa                                                  >  1:6789573/1‑93 (MQ=255)

                                          gccgcaATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                <  2:6813472/93‑1 (MQ=255)

                                           ccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                               <  2:3175275/93‑1 (MQ=255)

                                           ccgccATTCCCCAGGTGCCCGCCGTTAGTGGTTTTAGCGACCCGCAGT‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                               <  2:6149570/93‑1 (MQ=255)

                                            cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAttt                                               <  1:7931597/93‑1 (MQ=255)

                                            cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTc                                              <  1:6177891/93‑1 (MQ=255)

                                             gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                             <  2:3444004/93‑1 (MQ=255)

                                             gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                             >  2:2912858/1‑93 (MQ=255)

                                              ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                            <  2:1357415/93‑1 (MQ=255)

                                              ccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                            >  1:5206835/1‑93 (MQ=255)

                                                aTTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGa                                          >  2:7156923/1‑93 (MQ=255)

                                                 ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                         <  1:4426201/93‑1 (MQ=255)

                                                 ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                         >  2:4386397/1‑93 (MQ=255)

                                                 ttCTCCAGGTTGCCGCCGATAGTGGGTTTATCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                         >  1:7063444/1‑93 (MQ=255)

                                                    tcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                      <  1:5956478/93‑1 (MQ=255)

                                                    tcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                      >  1:1284684/1‑93 (MQ=255)

                                                     ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                     >  1:430193/1‑93 (MQ=255)

                                                      cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                    <  1:2385892/93‑1 (MQ=255)

                                                      cAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                    >  1:7773953/1‑93 (MQ=255)

                                                       aGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcaca                                   <  1:7244832/93‑1 (MQ=255)

                                                       aGGTCGCCGCCGAGAGTGGTTTTGGCGACAAGCAGTG‑CCCCCCCAGGCATATCGTTCATATCGCCTGCAGTGGTGACATTTCCCGACATcaca                                   <  2:4077618/93‑1 (MQ=255)

                                                        ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                  >  1:7280425/1‑93 (MQ=255)

                                                         gTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACAcc                                 >  1:6822643/1‑93 (MQ=255)

                                                          ttGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                <  1:7598985/93‑1 (MQ=255)

                                                          ttGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                >  2:1441416/1‑93 (MQ=255)

                                                           tGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                               <  2:6454147/93‑1 (MQ=255)

                                                           tGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                               <  1:2788032/93‑1 (MQ=255)

                                                           tGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGc                               <  2:387682/93‑1 (MQ=255)

                                                             ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                             >  1:3097528/1‑93 (MQ=255)

                                                             ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                             >  2:3534287/1‑93 (MQ=255)

                                                             ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                             <  1:7852975/93‑1 (MQ=255)

                                                             ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                             <  2:7365640/93‑1 (MQ=255)

                                                              cgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCCCACCGCTGt                            >  1:5621194/1‑93 (MQ=255)

                                                                ccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                           >  2:1543686/1‑93 (MQ=255)

                                                                ccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTa                          <  2:174812/93‑1 (MQ=255)

                                                                 cgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                         <  1:6562694/93‑1 (MQ=255)

                                                                 cgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                         >  2:3020176/1‑93 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                         <  1:1538896/93‑1 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        >  1:1737497/1‑93 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        >  1:1024419/1‑93 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        <  1:7226392/93‑1 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        >  2:4597569/1‑93 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        >  2:2598609/1‑93 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        <  1:5497910/93‑1 (MQ=255)

                                                                  gATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                        <  2:1914749/93‑1 (MQ=255)

                                                                       tGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGatat                   <  1:2996029/93‑1 (MQ=255)

                                                                        ggTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATAtt                  >  1:249382/1‑93 (MQ=255)

                                                                           tttAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt               <  1:4363467/93‑1 (MQ=255)

                                                                            ttAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg              <  2:1067104/93‑1 (MQ=255)

                                                                             tAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt             <  1:6226856/93‑1 (MQ=255)

                                                                               gCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATgtgt           <  2:5648914/93‑1 (MQ=255)

                                                                                cGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATgtgt           <  1:6643834/93‑1 (MQ=255)

                                                                                   cacGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc       <  1:2254768/93‑1 (MQ=255)

                                                                                    acGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACt      <  2:4904178/93‑1 (MQ=255)

                                                                                    acGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAAATCACACCGCTGTTATGGATATTAATGTGTTGACt      <  1:1376773/93‑1 (MQ=255)

                                                                                      gCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg    >  2:5396714/1‑93 (MQ=255)

                                                                                         gTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc  >  1:473031/1‑93 (MQ=255)

TACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCC  >  NC_000913/1592057‑1592238

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