BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,234,921	T→C	T285A (ACG→GCG)	mglC ←	MglC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,234,921	0	T	C	89.5%	656.8 / 54.4	247	T285A (ACG→GCG)	mglC	MglC
Reads supporting (aligned to +/- strand): ref base T (9/17); new base C (113/108); total (122/125)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

GATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGC  >  NC_000913/2234829‑2235013

gatgTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGt                                                                                              >  1:1657771/1‑93 (MQ=255)

  tgTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGccc                                                                                            >  1:1421431/1‑93 (MQ=255)

  tgTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGccc                                                                                            >  1:2884756/1‑93 (MQ=255)

   gTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGcccc                                                                                           >  1:825233/1‑93 (MQ=255)

   gTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGcccc                                                                                           <  2:6898885/93‑1 (MQ=255)

   gTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGcccc                                                                                           >  1:1936613/1‑93 (MQ=255)

      cTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCcac                                                                                        >  1:2728432/1‑93 (MQ=255)

         ccAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAcc                                                                                     <  1:549334/93‑1 (MQ=255)

         ccAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAcc                                                                                     >  2:4177570/1‑93 (MQ=255)

          cAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAccg                                                                                    <  1:3985619/93‑1 (MQ=255)

          cAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAccg                                                                                    >  1:7067445/1‑93 (MQ=255)

            gTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACAccgcc                                                                                  >  1:373275/1‑93 (MQ=255)

            gTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAccgcc                                                                                  >  2:246849/1‑93 (MQ=255)

             tATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAccgccg                                                                                 >  2:6316480/1‑93 (MQ=255)

             tATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACAccgccg                                                                                 >  1:5396140/1‑93 (MQ=255)

              aTGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACcgccgc                                                                                <  2:6561583/93‑1 (MQ=255)

               tGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCt                                                                               >  1:904165/1‑93 (MQ=255)

               tGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCt                                                                               <  2:3049525/93‑1 (MQ=255)

                 ggTTTACGCCGATATACGTCAGGCCGTAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGa                                                                             <  2:455861/93‑1 (MQ=255)

                  gTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGaa                                                                            >  2:4668959/1‑93 (MQ=255)

                   tttACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAAc                                                                           <  1:1783140/93‑1 (MQ=255)

                   tttACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAAc                                                                           <  2:1282843/93‑1 (MQ=255)

                       cGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATa                                                                       <  2:2032064/93‑1 (MQ=255)

                       cGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATa                                                                       <  2:3481980/93‑1 (MQ=255)

                        gCCGATATACATCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATac                                                                      >  1:1702337/1‑93 (MQ=255)

                         ccGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATaca                                                                     <  2:2397630/93‑1 (MQ=255)

                         ccGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATaca                                                                     <  1:6470313/93‑1 (MQ=255)

                         ccGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATaca                                                                     >  1:3298836/1‑93 (MQ=255)

                          cGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATacac                                                                    >  2:2063870/1‑93 (MQ=255)

                           gATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACAcc                                                                   <  1:5424525/93‑1 (MQ=255)

                             tataCGTCAGGCCAAAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCNc                                                                 >  2:5182338/1‑93 (MQ=255)

                              ataCGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGcc                                                                <  1:4419442/93‑1 (MQ=255)

                               taCGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCt                                                               >  1:4947984/1‑93 (MQ=255)

                                 cGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTAc                                                             >  2:933015/1‑93 (MQ=255)

                                  gTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTAcc                                                            >  2:4250358/1‑93 (MQ=255)

                                   tCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACca                                                           >  2:1469822/1‑93 (MQ=255)

                                   tCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACca                                                           >  1:7355880/1‑93 (MQ=255)

                                    cAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACcac                                                          <  1:5917676/93‑1 (MQ=255)

                                     aGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACcacg                                                         >  1:3636275/1‑93 (MQ=255)

                                      ggCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACcacgc                                                        >  2:1585885/1‑93 (MQ=255)

                                      ggCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACcacgc                                                        >  2:5755174/1‑93 (MQ=255)

                                        ccATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACcacgcac                                                      >  1:2999897/1‑93 (MQ=255)

                                         cATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACcacgcacg                                                     >  2:52748/1‑93 (MQ=255)

                                          aTAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCacgcacgc                                                    <  1:5602065/93‑1 (MQ=255)

                                          aTAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCacgcacgc                                                    >  1:3512270/1‑93 (MQ=255)

                                          aTAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCacgcacgc                                                    <  2:5615901/93‑1 (MQ=255)

                                           tAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgcc                                                   >  2:1207379/1‑93 (MQ=255)

                                           tAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgcc                                                   >  2:4751272/1‑93 (MQ=255)

                                           tAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgcc                                                   >  1:6939980/1‑93 (MQ=255)

                                            aGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgccg                                                  >  1:4860612/1‑93 (MQ=255)

                                             gTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCAcgccgc                                                 >  1:7480512/1‑93 (MQ=255)

                                                gatgaCGGTAAAAATAATTCCCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGAt                                              <  1:202186/93‑1 (MQ=255)

                                                gatgaCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGAt                                              <  2:4971655/93‑1 (MQ=255)

                                                gatgaCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGAt                                              <  2:2224778/93‑1 (MQ=255)

                                                gatgaCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGAt                                              <  2:4207914/93‑1 (MQ=255)

                                                 atgaCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATa                                             <  2:4486312/93‑1 (MQ=255)

                                                  tgaCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAg                                            >  1:3820509/1‑93 (MQ=255)

                                                   gaCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGc                                           >  2:5408664/1‑93 (MQ=255)

                                                    aCGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCa                                          <  2:2810690/93‑1 (MQ=255)

                                                     cGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCAt                                         >  2:5531117/1‑93 (MQ=255)

                                                        tAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                      <  1:4820996/93‑1 (MQ=255)

                                                        tAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCa                                      <  2:334671/93‑1 (MQ=255)

                                                        tAAAAATAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGACTACCACGCACGCCGCGATAGCATCCa                                      >  2:3563557/1‑93 (MQ=255)

                                                         aaaaaTAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAg                                     >  1:3126317/1‑93 (MQ=255)

                                                          aaaaTAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGc                                    <  2:7454467/93‑1 (MQ=255)

                                                           aaaTAATTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCt                                   >  2:1800117/1‑93 (MQ=255)

                                                             ataatTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTca                                 >  1:1169153/1‑93 (MQ=255)

                                                              taatTCCCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcat                                <  1:2552194/93‑1 (MQ=255)

                                                               aatTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcata                               <  1:7155641/93‑1 (MQ=255)

                                                                atTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatac                              <  1:1014890/93‑1 (MQ=255)

                                                                atTACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatac                              <  1:6357429/93‑1 (MQ=255)

                                                                 ttACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcataca                             <  1:5080714/93‑1 (MQ=255)

                                                                 ttACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcataca                             >  2:2581892/1‑93 (MQ=255)

                                                                 ttACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcataca                             >  2:1289965/1‑93 (MQ=255)

                                                                 ttACCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcataca                             <  1:6046324/93‑1 (MQ=255)

                                                                   aCCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacata                           >  2:6958405/1‑93 (MQ=255)

                                                                   aCCCCGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTcatacata                           <  2:831437/93‑1 (MQ=255)

                                                                    ccccGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATTCAGCTCATACATaa                          <  2:6440733/93‑1 (MQ=255)

                                                                    ccccGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaa                          <  1:223408/93‑1 (MQ=255)

                                                                    ccccGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaa                          <  2:5910912/93‑1 (MQ=255)

                                                                     cccGGTCACCACGCCAATCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATaaa                         <  1:3784803/93‑1 (MQ=255)

                                                                      ccGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAg                        >  2:1665989/1‑93 (MQ=255)

                                                                       cGGTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGc                       >  1:358339/1‑93 (MQ=255)

                                                                         gTCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCg                     >  1:4702049/1‑93 (MQ=255)

                                                                          tCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGa                    <  2:4796926/93‑1 (MQ=255)

                                                                          tCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGa                    <  2:7321818/93‑1 (MQ=255)

                                                                          tCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGa                    <  2:4078677/93‑1 (MQ=255)

                                                                          tCACCACGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGa                    <  2:2710722/93‑1 (MQ=255)

                                                                              cacGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgtt                >  1:929188/1‑93 (MQ=255)

                                                                              cacGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgtt                >  2:7210512/1‑93 (MQ=255)

                                                                              cacGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgtt                >  2:669625/1‑93 (MQ=255)

                                                                               acGCCAATCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGgttg               <  2:1051842/93‑1 (MQ=255)

                                                                                    aaTCCCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggt          <  1:5305270/93‑1 (MQ=255)

                                                                                    aaTCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggt          >  2:6173975/1‑93 (MQ=255)

                                                                                    aaTCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggt          >  1:6723599/1‑93 (MQ=255)

                                                                                     aTCACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggtg         <  2:3978136/93‑1 (MQ=255)

                                                                                      tCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTtggtgg        >  1:4606199/1‑93 (MQ=255)

                                                                                       cACCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGc       <  1:1336052/93‑1 (MQ=255)

                                                                                        aCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCa      <  1:1744756/93‑1 (MQ=255)

                                                                                        aCCGCCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCa      >  2:3722317/1‑93 (MQ=255)

GATGTACTGCCAGTATGGGTTTACGCCGATATACGTCAGGCCATAGTTGATGACGGTAAAAATAATTACCCCGGTCACCACGCCAATCACCGTCCCCACACCGCCGCTGAACGATACACCGCCTACCACGCACGCCGCGATAGCATCCAGCTCATACATAAAGCCGAGGTTGTTGGTGGCAGAGC  >  NC_000913/2234829‑2235013

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