BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	932,454	A→G	Q3R (CAG→CGG)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	932,454	0	A	G	98.4%	391.7 / ‑2.4	123	Q3R (CAG→CGG)	ftsK	cell division protein
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (44/77); total (44/79)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

Only 100 of 125 total aligned reads displayed.

ACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGT  >  NC_000913/932365‑932543

 cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGaa                                                                                       <  2:4144771/93‑1 (MQ=255)

 cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGaa                                                                                       <  1:822186/93‑1 (MQ=255)

 cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGaa                                                                                       <  1:667157/93‑1 (MQ=255)

   gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATa                                                                                     <  1:4306995/93‑1 (MQ=255)

   gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATa                                                                                     >  1:1602293/1‑93 (MQ=255)

     aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   <  1:3013255/93‑1 (MQ=255)

     aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   <  1:3874570/93‑1 (MQ=255)

     aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   >  1:690014/1‑93 (MQ=255)

     aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   <  1:3730668/93‑1 (MQ=255)

      aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  <  2:3008754/93‑1 (MQ=255)

      aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  >  1:4097373/1‑93 (MQ=255)

      aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  >  2:4134174/1‑93 (MQ=255)

       cAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAtt                                                                                 <  2:3586049/93‑1 (MQ=255)

       cAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAtt                                                                                 <  1:1893098/93‑1 (MQ=255)

        aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                <  2:2617403/93‑1 (MQ=255)

        aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                <  2:4017928/93‑1 (MQ=255)

          gaCTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGaa                                                                              <  2:4719497/92‑1 (MQ=255)

            cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGa                                                                            >  2:3238749/1‑93 (MQ=255)

             tGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGAc                                                                           >  1:557596/1‑93 (MQ=255)

              ggggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACa                                                                          >  1:4517009/1‑93 (MQ=255)

               gggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGCGCCGGGAATACATTGAAGACaa                                                                         <  2:2990177/93‑1 (MQ=255)

               gggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACaa                                                                         <  1:2377341/93‑1 (MQ=255)

                ggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACaaa                                                                        <  1:2595977/93‑1 (MQ=255)

                ggTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACaaa                                                                        <  1:3607776/93‑1 (MQ=255)

                 gTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaag                                                                       <  2:1240061/93‑1 (MQ=255)

                  tAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaaga                                                                      <  1:4454647/93‑1 (MQ=255)

                   aaCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaa                                                                     >  1:1750212/1‑93 (MQ=255)

                     ccTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGt                                                                   >  2:78491/1‑93 (MQ=255)

                      cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  <  2:2918470/93‑1 (MQ=255)

                      cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  <  1:4679749/93‑1 (MQ=255)

                      cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  <  2:3642209/93‑1 (MQ=255)

                      cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  <  2:3642210/93‑1 (MQ=255)

                      cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGACTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  >  2:4163483/1‑93 (MQ=255)

                         gTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcaca                                                               >  1:2787250/1‑93 (MQ=255)

                           aCTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACAtt                                                             >  2:4129953/1‑93 (MQ=255)

                           aCTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACAtt                                                             >  1:4398950/1‑93 (MQ=255)

                            cTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTg                                                            >  1:318478/1‑93 (MQ=255)

                             tgttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           <  2:2173452/93‑1 (MQ=255)

                             tgttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           <  2:340487/93‑1 (MQ=255)

                             tgttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           <  1:3898336/93‑1 (MQ=255)

                               ttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACa                                                         >  1:4120292/1‑93 (MQ=255)

                               ttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACa                                                         >  2:4604183/1‑93 (MQ=255)

                                  tCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAg                                                      >  1:1717742/1‑93 (MQ=255)

                                   ccGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGt                                                     <  1:4140097/93‑1 (MQ=255)

                                     gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAGTACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   >  1:3233569/1‑93 (MQ=255)

                                     gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   <  1:3145705/93‑1 (MQ=255)

                                     gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   <  2:3506650/93‑1 (MQ=255)

                                      ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  <  2:4452585/93‑1 (MQ=255)

                                      ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  >  2:2369342/1‑93 (MQ=255)

                                       tttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAg                                                 <  2:1022619/93‑1 (MQ=255)

                                        ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                                <  1:2015690/93‑1 (MQ=255)

                                        ttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                                <  2:4242018/93‑1 (MQ=255)

                                         tAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTa                                               <  1:2432379/93‑1 (MQ=255)

                                         tAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTa                                               <  2:2730987/93‑1 (MQ=255)

                                          aGCATCGGGCAGGAAAAGCCTGTAACCTGGAGACCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAg                                              <  2:4605921/93‑1 (MQ=255)

                                           gCATCGGGCAGGAAAAGCCTGAAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGc                                             <  1:3410082/93‑1 (MQ=255)

                                            cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCg                                            <  1:2115726/93‑1 (MQ=255)

                                            cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCg                                            <  2:279030/93‑1 (MQ=255)

                                            cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCg                                            <  1:2293529/93‑1 (MQ=255)

                                            cATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCg                                            <  1:3076885/93‑1 (MQ=255)

                                             aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCgg                                           <  1:1462547/93‑1 (MQ=255)

                                             aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCgg                                           <  1:3164968/93‑1 (MQ=255)

                                             aTCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCgg                                           >  2:4344296/1‑93 (MQ=255)

                                              tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          >  2:2116533/1‑93 (MQ=255)

                                              tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          <  1:4035900/93‑1 (MQ=255)

                                               cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgcc                                         <  2:2875507/93‑1 (MQ=255)

                                                 ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgc                                       <  1:2928501/93‑1 (MQ=255)

                                                    aGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgc                                    <  1:826151/93‑1 (MQ=255)

                                                     ggAAAAGCCTGTAACCTGGAGAGCTTTTGTTGATCCGGGAATACATTGAAGACAAAGAAGTAACATTGACAAAGTTAAGTAGCGgccgccgcc                                   <  2:4585186/93‑1 (MQ=255)

                                                     ggAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                                   <  2:620546/93‑1 (MQ=255)

                                                     ggAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                                   <  2:1973629/93‑1 (MQ=255)

                                                       aaaaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCtt                                 <  1:4505597/93‑1 (MQ=255)

                                                        aaaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                                <  2:1751632/93‑1 (MQ=255)

                                                          aGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                              >  2:3671889/1‑93 (MQ=255)

                                                          aGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                              <  2:4035339/93‑1 (MQ=255)

                                                          aGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                              >  2:3671888/1‑93 (MQ=255)

                                                            ccTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa                            <  2:1944098/93‑1 (MQ=255)

                                                            ccTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa                            <  1:1312989/93‑1 (MQ=255)

                                                             cTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg                           <  2:2128506/93‑1 (MQ=255)

                                                                 aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGt                       <  2:3924703/93‑1 (MQ=255)

                                                                 aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGt                       <  1:2297214/93‑1 (MQ=255)

                                                                  aCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGtt                      <  1:4652057/93‑1 (MQ=255)

                                                                    cTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGc                    >  1:2713457/1‑93 (MQ=255)

                                                                      ggAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTg                  >  2:595701/1‑93 (MQ=255)

                                                                      ggAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTg                  <  2:4331875/93‑1 (MQ=255)

                                                                       gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 <  1:1176118/93‑1 (MQ=255)

                                                                       gagaGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa                 <  1:461828/93‑1 (MQ=255)

                                                                           gCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCt             >  2:1425831/1‑93 (MQ=255)

                                                                            ccTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCtt            <  2:1692975/93‑1 (MQ=255)

                                                                              tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  2:1384627/1‑93 (MQ=255)

                                                                              tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  1:1272035/1‑93 (MQ=255)

                                                                              tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  1:2002628/1‑93 (MQ=255)

                                                                               ttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCTGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAtt         >  1:3863682/1‑93 (MQ=255)

                                                                               ttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAtt         <  2:3515695/93‑1 (MQ=255)

                                                                                tCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTg        <  2:4597570/93‑1 (MQ=255)

                                                                                tCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGATGATCCTTATTg        >  1:2180816/1‑93 (MQ=255)

                                                                                   tGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTcc     <  1:56238/93‑1 (MQ=255)

                                                                                    gAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCt    >  1:2877049/1‑93 (MQ=255)

                                                                                     aGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTg   <  1:2711482/93‑1 (MQ=255)

                                                                                     aGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTg   <  2:448292/93‑1 (MQ=255)

ACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGT  >  NC_000913/932365‑932543

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