BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,631,832	(T)_6→5	coding (259/1467 nt)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,631,827	0	T	.	99.7%	1395.5 / ‑3.1	308	coding (264/1467 nt)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base . (152/155); total (153/155)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 337 total aligned reads displayed.

GCAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAA  >  NC_000913/2631735‑2631915

gCAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCgtgt                                                                                          >  1:2700546/1‑93 (MQ=255)

 cAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtt                                                                                         >  1:5104267/1‑93 (MQ=255)

     aCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCCGCTTCGTGtttttc                                                                                     >  2:572514/1‑92 (MQ=255)

        gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcaca                                                                                  <  1:598062/93‑5 (MQ=255)

        gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcaca                                                                                  >  1:2189419/1‑89 (MQ=255)

         ttACGCTCGGTCAGTTCTTTCACTTAGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCacac                                                                                >  1:4930464/1‑93 (MQ=255)

           aCGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACgg                                                                              >  2:2621609/1‑93 (MQ=255)

           aCGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACgg                                                                              >  1:5703947/1‑93 (MQ=255)

              cTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGa                                                                           <  2:651361/93‑1 (MQ=255)

               tCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGaa                                                                          >  1:8606603/1‑93 (MQ=255)

                cGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAAc                                                                         <  1:548031/93‑1 (MQ=255)

                 ggTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAAct                                                                        >  1:1644488/1‑93 (MQ=255)

                 ggTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAAct                                                                        >  2:2561966/1‑93 (MQ=255)

                   tCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActtc                                                                      <  2:7592057/93‑1 (MQ=255)

                     aGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActtctt                                                                    >  2:2984169/1‑93 (MQ=255)

                     aGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActtctt                                                                    <  2:7044232/93‑1 (MQ=255)

                      gTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCTCACGGCGAACttcttc                                                                   <  2:9559756/93‑1 (MQ=255)

                      gTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACttcttc                                                                   <  2:3183758/93‑1 (MQ=255)

                      gTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACttcttc                                                                   >  2:886142/1‑93 (MQ=255)

                       ttCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTtcttct                                                                  >  2:513665/1‑93 (MQ=255)

                       ttCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTtcttct                                                                  >  2:4426103/1‑93 (MQ=255)

                         cTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGc                                                                >  2:3786720/1‑93 (MQ=255)

                          tttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGcc                                                               <  1:3406547/93‑1 (MQ=255)

                          tttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGcc                                                               <  2:4477868/93‑1 (MQ=255)

                           ttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCt                                                              <  1:373221/93‑1 (MQ=255)

                           ttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCt                                                              <  1:8873418/93‑1 (MQ=255)

                             cccTTCGCGCAGCGGCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTgg                                                            <  2:1529670/91‑1 (MQ=255)

                              aCTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGc                                                           <  1:3880800/93‑1 (MQ=255)

                               cTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCATATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCg                                                          >  2:5923969/1‑93 (MQ=255)

                                ttCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGt                                                         <  1:9807351/93‑1 (MQ=255)

                                 tCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGtt                                                        <  2:9541234/93‑1 (MQ=255)

                                  cgcgCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTc                                                       <  1:2147947/93‑1 (MQ=255)

                                    cgcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCaa                                                     <  1:6386984/93‑1 (MQ=255)

                                    cgcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCaa                                                     <  1:6728688/93‑1 (MQ=255)

                                     gcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAAt                                                    >  1:241777/1‑93 (MQ=255)

                                      cAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATg                                                   >  1:8678858/1‑93 (MQ=255)

                                        gCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCACTGGa                                                 >  2:1437925/1‑93 (MQ=255)

                                        gCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGa                                                 >  2:9356032/1‑93 (MQ=255)

                                         cgtcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGAc                                                <  1:812612/93‑1 (MQ=255)

                                          gtcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACa                                               <  2:20625/93‑1 (MQ=255)

                                          gtcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACa                                               <  1:324139/93‑1 (MQ=255)

                                          gtcgtGGGTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACa                                               <  1:6874228/93‑1 (MQ=255)

                                           tcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACAt                                              >  2:3799954/1‑93 (MQ=255)

                                           tcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACAt                                              >  1:2934469/1‑93 (MQ=255)

                                            cgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATg                                             <  2:424537/93‑1 (MQ=255)

                                            cgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATg                                             <  2:1409312/93‑1 (MQ=255)

                                             gtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGt                                            <  2:891538/93‑1 (MQ=255)

                                             gtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGt                                            >  2:7975253/1‑93 (MQ=255)

                                             gtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGt                                            <  2:7214003/93‑1 (MQ=255)

                                             gtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGt                                            >  2:1253726/1‑93 (MQ=255)

                                               ggTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGttt                                          >  2:4579050/1‑93 (MQ=255)

                                                gTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGtttt                                         <  2:9344104/93‑1 (MQ=255)

                                                gTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGtttt                                         >  1:1142307/1‑93 (MQ=255)

                                                gTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGtttt                                         >  1:3635920/1‑93 (MQ=255)

                                                gTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGtttt                                         >  1:3169739/1‑93 (MQ=255)

                                                  tGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTtg                                       <  2:4880422/93‑1 (MQ=255)

                                                   ggCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTtgt                                      <  1:1961438/93‑1 (MQ=255)

                                                      aGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGa                                   <  2:3738847/93‑1 (MQ=255)

                                                       gAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGAt                                  >  1:7980838/1‑93 (MQ=255)

                                                       aaaCAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGAt                                  <  1:3521525/92‑1 (MQ=255)

                                                         aCAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaa                                <  2:6230858/93‑1 (MQ=255)

                                                         aCAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaa                                <  1:4803869/93‑1 (MQ=255)

                                                          cAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaaa                               <  1:10045586/93‑1 (MQ=255)

                                                           aGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaagg                              >  1:7753039/1‑91 (MQ=255)

                                                           aGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAg                              >  2:7843322/1‑93 (MQ=255)

                                                            gTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGc                             >  1:1242728/1‑93 (MQ=255)

                                                             tCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGcc                            >  1:4259248/1‑93 (MQ=255)

                                                               tGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGa                          >  1:5442266/1‑93 (MQ=255)

                                                               tGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGa                          >  2:5412497/1‑93 (MQ=255)

                                                                gCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGAt                         <  1:8270344/93‑1 (MQ=255)

                                                                 cGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATa                        <  1:4364974/93‑1 (MQ=255)

                                                                  ggATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAc                       >  2:4343346/1‑93 (MQ=255)

                                                                   gATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAcc                      >  2:9136013/1‑93 (MQ=255)

                                                                   gATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAcc                      <  1:753953/93‑1 (MQ=255)

                                                                   gATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAcc                      >  1:4451321/1‑93 (MQ=255)

                                                                    aTCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGACTGGCGTTCAATGTACATGTTTTTGCGGATAAAGCCGATAcc                      >  1:9219610/1‑93 (MQ=255)

                                                                     tCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGc                    >  1:4034120/1‑93 (MQ=255)

                                                                      cAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGcc                   <  1:7420594/93‑1 (MQ=255)

                                                                       aGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCt                  <  2:2403862/93‑1 (MQ=255)

                                                                        gTCACCCCACCAGATTCGTG‑TTTTTCACACGTCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCtt                 >  2:3511563/1‑93 (MQ=255)

                                                                         tCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTc                <  2:5104266/93‑1 (MQ=255)

                                                                         tCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTc                >  1:1712900/1‑93 (MQ=255)

                                                                          caccacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTcc               <  1:5660096/93‑1 (MQ=255)

                                                                          caccacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTcc               <  1:6237028/93‑1 (MQ=255)

                                                                             cacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCTCCTTCCTGa            >  2:8871799/1‑93 (MQ=255)

                                                                             cacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGa            <  1:9160956/93‑1 (MQ=255)

                                                                              acaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAg           >  2:6860708/1‑93 (MQ=255)

                                                                               caccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGc          >  2:8450726/1‑93 (MQ=255)

                                                                               caccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGc          >  2:1928963/1‑93 (MQ=255)

                                                                                 ccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCa        >  1:7203369/1‑93 (MQ=255)

                                                                                  caGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCag       <  1:8773873/93‑1 (MQ=255)

                                                                                   cgATTCGTG‑TTTTTCACACGGCGAGCTTCCTCGGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCaga      <  2:2193993/92‑1 (MQ=255)

                                                                                   aGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCaga      <  1:5562540/93‑1 (MQ=255)

                                                                                   aGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCaga      <  1:6612560/93‑1 (MQ=255)

                                                                                     aTTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGc    <  2:6434416/93‑1 (MQ=255)

                                                                                     aTTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGc    <  1:6049871/93‑1 (MQ=255)

                                                                                     aTTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGc    <  1:6483330/93‑1 (MQ=255)

                                                                                       tCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCac  >  2:6081686/1‑92 (MQ=255)

                                                                                       tCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCaa  <  2:5285882/93‑1 (MQ=255)

                                                                                       tCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCaa  >  2:8125457/1‑93 (MQ=255)

GCAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAA  >  NC_000913/2631735‑2631915

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