BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	867,645	(C)_8→7	coding (870/1839 nt)	b0829 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	867,638	0	C	.	98.9%	831.9 / ‑1.4	188	coding (863/1839 nt)	b0829	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/2); new base . (95/91); total (95/93)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

TTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGT  >  NC_000913/867547‑867726

ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcc                                                                                         <  2:6738956/93‑1 (MQ=255)

  aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccc                                                                                       >  1:7285133/1‑93 (MQ=255)

   gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccc                                                                                      <  1:6665244/93‑1 (MQ=255)

   gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccc                                                                                      <  2:2888566/93‑1 (MQ=255)

     tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccc                                                                                    >  1:3505361/1‑93 (MQ=255)

     tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccc                                                                                    <  1:2109055/93‑1 (MQ=255)

      gctgTTCCGCAACTTGGTGGGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccccca                                                                                   <  2:1503432/93‑2 (MQ=255)

         gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatcg                                                                                >  1:1994313/1‑89 (MQ=255)

         gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatcg                                                                                >  1:7683048/1‑89 (MQ=255)

         gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGcccccccatcg                                                                                >  1:4973215/1‑89 (MQ=255)

         gTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGccccaccaccg                                                                                >  2:5853982/1‑89 (MQ=255)

          ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGa                                                                              <  2:5510641/93‑1 (MQ=255)

          ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGa                                                                              <  1:5379094/93‑1 (MQ=255)

           tCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAg                                                                             >  1:3254794/1‑93 (MQ=255)

              gCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAg                                                                          >  1:3183174/1‑93 (MQ=255)

               cAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGa                                                                         <  2:350555/93‑1 (MQ=255)

                 aCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaa                                                                       <  1:3103362/93‑1 (MQ=255)

                  cTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                      <  2:4223080/93‑1 (MQ=255)

                  cTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                      <  1:7168385/93‑1 (MQ=255)

                  cTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                      <  2:489059/93‑1 (MQ=255)

                  cTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaaa                                                                      <  2:5159706/93‑1 (MQ=255)

                   ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGaaacc                                                                     >  1:5658194/1‑91 (MQ=255)

                   ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGACAAc                                                                     >  1:7408718/1‑93 (MQ=255)

                   ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAAc                                                                     <  2:6629917/93‑1 (MQ=255)

                   ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAAc                                                                     <  2:4929693/93‑1 (MQ=255)

                   ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAAc                                                                     <  1:6761458/93‑1 (MQ=255)

                     ggTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACgg                                                                   >  1:4020805/1‑93 (MQ=255)

                     ggTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACgg                                                                   >  2:3839913/1‑93 (MQ=255)

                      gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACggt                                                                  <  2:4463441/93‑1 (MQ=255)

                        gCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑ACCCCCCATCGAGCAGAAAACggtgg                                                                <  2:1449286/93‑1 (MQ=255)

                           aTGAAAGGGTTAAATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGtggtgg                                                             <  2:6800902/93‑1 (MQ=255)

                            tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                            <  1:4690958/93‑1 (MQ=255)

                            tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                            >  2:1114227/1‑93 (MQ=255)

                             gAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGAt                                                           <  1:1667459/93‑1 (MQ=255)

                             gAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGAt                                                           <  2:7369410/93‑1 (MQ=255)

                             gAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGAt                                                           <  2:1136684/93‑1 (MQ=255)

                               aaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                         <  1:7823289/93‑1 (MQ=255)

                                  gtttAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGa                                                      <  2:5482978/91‑1 (MQ=255)

                                  ggTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGa                                                      >  1:1903521/1‑93 (MQ=255)

                                  ggTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGa                                                      >  1:5400769/1‑93 (MQ=255)

                                   gTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                     <  2:7246941/93‑1 (MQ=255)

                                      aGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCt                                                  <  2:6281308/93‑1 (MQ=255)

                                       gATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTg                                                 >  1:1354295/1‑93 (MQ=255)

                                        aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                <  1:682165/93‑1 (MQ=255)

                                         ttATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                               >  1:3300288/1‑93 (MQ=255)

                                         tcaTCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                               <  1:2822386/91‑1 (MQ=255)

                                         atATCCCCGACGGTTCCCGTTGATATCGGCTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                               <  1:1352000/92‑1 (MQ=255)

                                          tATTCCCGCCGGTTCACGGTGCTATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCGATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                              <  1:1910878/93‑1 (MQ=255)

                                          tATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                              <  1:2773176/93‑1 (MQ=255)

                                          tATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                              >  1:363799/1‑93 (MQ=255)

                                             ccccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                           <  1:833685/93‑1 (MQ=255)

                                             ccccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                           <  1:2647629/93‑1 (MQ=255)

                                             ccccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                           <  1:1593158/93‑1 (MQ=255)

                                               ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                         >  2:6929975/1‑93 (MQ=255)

                                               ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                         >  2:2144486/1‑93 (MQ=255)

                                                cGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAg                                        >  1:6530154/1‑93 (MQ=255)

                                                 gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAgg                                       >  1:2279283/1‑92 (MQ=255)

                                                 gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAgg                                       >  2:2274860/1‑92 (MQ=255)

                                                 gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGt                                       >  1:2898783/1‑93 (MQ=255)

                                                   cGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                     >  2:6200864/1‑93 (MQ=255)

                                                     tttCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGTACCTGTTTTACGAGTGCGt                                   >  1:7022976/1‑93 (MQ=255)

                                                       tCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTaa                                 <  2:3579065/93‑1 (MQ=255)

                                                         ccGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGTAAACGGTGGTGGATGGCGAACCTGTTTTACGAGCGCGTACTc                               >  2:4690095/1‑93 (MQ=255)

                                                         ccGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATc                               >  2:7178238/1‑93 (MQ=255)

                                                          cGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCt                              <  1:4127847/93‑1 (MQ=255)

                                                          cGTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCt                              <  2:2822440/93‑1 (MQ=255)

                                                           gTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                             <  2:1346707/93‑1 (MQ=255)

                                                           gTTGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                             >  2:5352832/1‑93 (MQ=255)

                                                           gTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTg                            >  1:1892444/1‑93 (MQ=255)

                                                            ttGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTg                            >  1:2191977/1‑93 (MQ=255)

                                                            ttGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTg                            <  1:474604/93‑1 (MQ=255)

                                                             tGATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGt                           >  2:600140/1‑93 (MQ=255)

                                                              gATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGGACCTGTTTTTCGAGGGCGTAATTTTGta                          >  2:2982320/1‑92 (MQ=255)

                                                              gATATCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTc                          >  1:7586988/1‑93 (MQ=255)

                                                               atatCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCa                         >  1:5730613/1‑93 (MQ=255)

                                                                  tCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAccc                      >  2:2492840/1‑93 (MQ=255)

                                                                  tCGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAccc                      >  1:5310930/1‑93 (MQ=255)

                                                                   cGCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCg                     <  1:4389259/93‑1 (MQ=255)

                                                                    gCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGt                    <  1:378322/93‑1 (MQ=255)

                                                                    gCTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGt                    <  2:2718469/93‑1 (MQ=255)

                                                                     cTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                   <  2:3359619/93‑1 (MQ=255)

                                                                     cTTGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                   <  1:7907075/93‑1 (MQ=255)

                                                                      ttGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGttt                  <  1:7133317/93‑1 (MQ=255)

                                                                      ttGAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGttt                  <  2:2636164/93‑1 (MQ=255)

                                                                        gAACATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTcc                >  1:2050444/1‑93 (MQ=255)

                                                                         aaCATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTccc               >  1:555556/1‑93 (MQ=255)

                                                                           cATCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt             <  1:6173317/93‑1 (MQ=255)

                                                                            aTCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt            <  1:7138304/93‑1 (MQ=255)

                                                                            aTCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt            >  1:112405/1‑93 (MQ=255)

                                                                             tCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTg           >  1:6040881/1‑93 (MQ=255)

                                                                             tCCAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTg           >  2:4009757/1‑93 (MQ=255)

                                                                              ccAGCCAAACAGGCCCCCCCC‑TCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgc          <  2:5783095/93‑1 (MQ=255)

                                                                               cAGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcg         >  1:4539244/1‑93 (MQ=255)

                                                                                aGCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcgc        >  1:1142308/1‑93 (MQ=255)

                                                                                 gCGAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTTACCCGGTTCCCTTTGCGCa       >  2:4331066/1‑93 (MQ=255)

                                                                                   gAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGcgccgc     >  2:3119403/1‑93 (MQ=255)

                                                                                   gAAACAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGc     <  2:7740404/93‑1 (MQ=255)

                                                                                    aaaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg    >  2:5253119/1‑93 (MQ=255)

                                                                                    aaaCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCg    >  2:1462025/1‑93 (MQ=255)

                                                                                      aCAGG‑CCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGt  >  2:1945829/1‑93 (MQ=255)

TTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGT  >  NC_000913/867547‑867726

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