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	.	98.7%	672.2 / ‑0.1	150	coding (264/1467 nt)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (69/79); total (70/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 175 total aligned reads displayed.

GCAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAA  >  NC_000913/2631735‑2631915

gCAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCgtgt                                                                                          <  1:3073927/93‑1 (MQ=255)

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

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

  aaaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGttt                                                                                        <  2:2205888/93‑1 (MQ=255)

  aaaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGttt                                                                                        >  1:2743004/1‑93 (MQ=255)

   aaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttt                                                                                       >  2:3737972/1‑93 (MQ=255)

   aaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttt                                                                                       <  2:2996624/93‑1 (MQ=255)

   aaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttt                                                                                       >  1:2339855/1‑93 (MQ=255)

   aaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttt                                                                                       <  1:1057854/93‑1 (MQ=255)

    aaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCCTACCAGATTCGTGttttt                                                                                      >  2:35363/1‑93 (MQ=255)

    aaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGttttt                                                                                      >  1:3407713/1‑93 (MQ=255)

           aCGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACgg                                                                              <  1:2501446/93‑1 (MQ=255)

            cGCTCGTGCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGc                                                                             <  1:687245/93‑1 (MQ=255)

              cTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGa                                                                           >  1:932375/1‑93 (MQ=255)

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

               tCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGaa                                                                          <  2:4659348/93‑1 (MQ=255)

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

                 ggTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAAct                                                                        <  2:813211/93‑1 (MQ=255)

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

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

                  gTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActt                                                                       <  1:3077047/93‑1 (MQ=255)

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

                    cAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActtct                                                                     <  1:4495349/93‑1 (MQ=255)

                    cAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActtct                                                                     <  1:4347780/93‑1 (MQ=255)

                    cAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAActtct                                                                     <  2:4220486/93‑1 (MQ=255)

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

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

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

                      gTTCTTTCACTTCGCGCAGCGTCGTGGTGGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACttcttc                                                                   <  1:551416/93‑1 (MQ=255)

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

                       ttCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTtcttct                                                                  <  1:2637059/93‑1 (MQ=255)

                       ttCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTtcttct                                                                  <  2:1640098/93‑1 (MQ=255)

                        tCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTg                                                                 <  1:1586214/93‑1 (MQ=255)

                        tCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTg                                                                 <  1:4637904/93‑1 (MQ=255)

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

                          tttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGcc                                                               >  1:1569178/1‑93 (MQ=255)

                          tttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGcc                                                               >  1:4657456/1‑93 (MQ=255)

                          tttCACTTCGCGCAGCGTCGTGGTTCGCAGGACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTAcc                                                               >  2:4284767/1‑93 (MQ=255)

                           ttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCt                                                              >  1:1461872/1‑93 (MQ=255)

                           ttCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCt                                                              >  2:2053097/1‑93 (MQ=255)

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

                                 tCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGtt                                                        >  1:289293/1‑93 (MQ=255)

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

                                     gcAGCGTCGTGGTTGGCAGAACCGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCact                                                    >  2:3968441/1‑91 (MQ=255)

                                     gcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAAt                                                    <  1:3296430/93‑1 (MQ=255)

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

                                      cAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCGCCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATg                                                   >  2:4067696/1‑93 (MQ=255)

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

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

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

                                              tGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGtt                                           <  2:624655/93‑1 (MQ=255)

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

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

                                                     cAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTgg                                    <  1:2910922/93‑1 (MQ=255)

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

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

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

                                                      aGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGa                                   >  2:1977455/1‑93 (MQ=255)

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

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

                                                        aaCAGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATa                                 <  2:2506205/93‑1 (MQ=255)

                                                          cAGTCTGCGGATCAGTGACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaaa                               <  1:473998/93‑1 (MQ=255)

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

                                                           aGTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGGTCAATGGACATGTTTTTGTGGATAAAg                              >  2:2983424/1‑93 (MQ=255)

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

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

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

                                                            gTCTGCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGc                             <  1:386609/93‑1 (MQ=255)

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

                                                                gCGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGAt                         >  1:2287313/1‑93 (MQ=255)

                                                                 cGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATa                        >  2:2485244/1‑93 (MQ=255)

                                                                 cGGATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATa                        >  2:1759617/1‑93 (MQ=255)

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

                                                                  ggATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAc                       <  1:287531/93‑1 (MQ=255)

                                                                  ggATCAGTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAc                       <  1:276075/93‑1 (MQ=255)

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

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

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

                                                                   gATCAGTCACCACACCAGATTCGTG‑TTTTTCAAACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAcc                      <  1:1134326/93‑1 (MQ=255)

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

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

                                                                        gTCACCACACCAGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCtt                 >  2:4485314/1‑93 (MQ=255)

                                                                          caccacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGACGATACCGCCTTcc               >  1:1996433/1‑93 (MQ=255)

                                                                           accacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCt              >  2:25591/1‑93 (MQ=255)

                                                                            ccacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTg             >  1:1879307/1‑93 (MQ=255)

                                                                             cacaccaGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGCTAAAGCCGATAACGCCTTCCTGa            >  2:1418494/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                                                   aGATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCaga      >  2:3358106/1‑93 (MQ=255)

                                                                                    gATTCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCagag     >  2:3238687/1‑93 (MQ=255)

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

                                                                                      ttCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCa   >  2:1146660/1‑93 (MQ=255)

                                                                                      ttCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCa   <  2:1870320/93‑1 (MQ=255)

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

                                                                                       tCGTG‑TTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCaa  <  1:2878191/93‑1 (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: `‑`